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1 RESUMEN 

El articulo presenta un estudio comparativo del diseno de edificios regulares, resueltos 
con marcos de acero ductiles con contravientos concentricos, a los que se les aplicaron las 
metodologias de analisis y diseno conforme al metodo de diseno convencional del cuerpo 
principal y el planteamiento alterno del Apendice A establecido en las Normas Tecnicas para 
Diseno por Sismo (NTCS-2004) del Reglamento del Distrito Federal. Los edificios estan 
ubicados en la zona lllb de las NTCS-2004 y se diseharon con un factor de comportamiento 
sismico Q=4 y asociados a distorsiones que satisfacen las restricciones por desplazamiento 
cercanas al limite maximo establecidos por las NTCS-2004. Se hicieron comparaciones de 
ambos disenos mediante analisis estaticos no lineales ante carga lateral monotona creciente 
(pushover) y analisis dinamicos no lineales ante acelerogramas asociados al sismo de diseno. 
Los resultados del estudio sugieren la imperiosa necesidad de realizar una revision integral de 
los criterios de diseno de las NTCS-2004 para este tipo de estructuras, dadas las grandes 
reservas de resistencia lateral que lideran a las estructuras a comportarse de manera casi 
elastica, y que los analisis pushover sugieren mecanismos ultimos de falla diferentes e 
incompatibles a los que se suponen en su diseno. 

2 SUMMARY 

The results of a study devoted to compare the design of regular 15-story steel buildings 
structured with special-moment resisting concentric braced frames (SMRCBFs) is presented. 
The SMRCBFs were designed with the traditional method advocated in the main body of the 
seismic norms of Mexico’s Federal District Code (NTCS-2004), as well as a more modern, 
alternative method advocated in its Appendix A. Subject building are located in zone lllb of 
NTCS-2004 and were designed assuming a seismic response modification factor Q=4. Design 
drifts for subject building were design close to the limiting drift established by the building code. 
Pushover analyses and nonlinear time-history dynamic analyses were used to evaluate the 
seismic performance of the resulting design. The results of such analyses show the urgent need 
of performing and in-depth review of the design criteria of NTCS-2004 for medium-rise and 
lowrise SMRCBFs in the soft soils of Mexico City, given the large strength reserves that this 
system develop for designs that fulfill the requirements for the code, as their simulated nonlinear 
dynamic response suggest a likely elastic behavior under the action of acceleration records 
related to the maximum credible earthquake. Besides, collapse mechanisms obtained from 
pushover analyses are very different from the strong-column, weak-beam weaker-brace 
assumed in a design that fulfills the requirements of NTCS-2004. 
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1. INTRODUCCION 

Con la aprobacion y publication de las nuevas Normas Tecnicas Complementarias 
para Diseno por Sismo (NTCS-2004 2004) del Reglamento para Construcciones del Distrito 
Federal de Mexico se introdujo oficialmente un nuevo procedimiento de diseno sismico en su 
Apendice Normativo A cuyo proposito es hacer mas transparente el proceso de diseno y liderar 
hacia una nueva generacion de reglamentos de diseno sismico en Mexico donde los conceptos 
de diseno por desempeno sean claros y transparentes. Entre otros aspectos, este apendice 
aporta requerimientos de diseno que contemplan: un manejo transparente y explicito de las 
demandas y suministros de deformacion en la estructura; un diseno dual que considera 
criterios de desempeno de Operacion Completa y Seguridad de Vida; y una microzonacion que 
permite una mejor caracterizacion numerica de las excitaciones slsmicas de diseno (Ordaz et 
at. 2003). 

Ordaz et at. (2003) senalan los fundamentos de la propuesta del Apendice Normativo 
A, entre los que destacan los nuevos espectros de diseno elastico, con los cuales se pretende 
representar de una manera mas realista los verdaderos niveles de demanda que se presentan 
ante los sismos de diseno y reflejar correctamente los niveles de amplificacion que se producen 
en la realidad para diferentes tipos de suelo (Figura 1b) con respecto al metodo convencional 
que se viene utilizando a partir del RCDF-87, donde los espectros elasticos de diseno no 
reflejan las verdaderas demandas elasticas (Figura la), sino que inherentemente estan 
reducidos por sobrerresistencia por las razones que se indican en Rosenblueth et al. (1989). 
Ordaz et al. (2003) explican detalladamente el procedimiento que se siguio para su obtencion, 
asi como las consideraciones para efectuar las reducciones por concepto de ductilidad (Q y Q') 
y sobrerresistencia (R). 



Aceleracion 





b) Espectros de diseno del Apendice A 



a) Espectros de diseno del cuerpo principal 
Figura 1 Representacion esquematica de los espectros de diseno del cuerpo principal y del 
Apendice A de las NTCS-2004 del Reglamento del Distrito Federal mexicano 



El nuevo procedimiento de diseno sismico establecido en el Apendice Normativo A 
tiene solidas bases conceptuales y empiricas para la definicion de los espectros de diseno y la 
interaccion suelo-estructura. Sin embargo, tambien es cierto que varias propuestas, como el 
factor de reduccion por sobrerresistencia, R, se basa en un muy reducido numero de estudios 
de marcos de concreto reforzado bidimensionales, por lo que, como se ha comprobado 
analitica (Tapia 2005) y experimentalmente (Vazquez 2005), la propuesta de valores de este 
factor no es representativa de otras estructuraciones, como marcos contraventeados de acero 
estructural (Tapia 2005) y de estructuras con base en muros de mamposteria confinada 
(Vazquez 2005). Ademas, la calibracion de las demandas de respuesta no lineal esperadas 
para disenos hechos con el Apendice A se hizo con base en sistemas de un grado de libertad 
considerando comportamiento elastoplastico perfecto (Reyes 2004, Ordaz 2004). 

Por lo tanto, la evaluacion del desempeno estructural de edificios con distintas 
estructuraciones disehados conforme al Apendice A y su comparacion con disenos alternatives 
hechos conforme al metodo tradicional de diseno del cuerpo principal de las NTCS-2004 es del 
mayor interes, para valorar los pros y contras que esta novedosa propuesta ofrece. 
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El presente trabajo presenta solamente una parte de un estudio muy detallado que se 
presenta en Tapia (2005) y donde se comparan distintos disenos de edificios de acero 
estructural de 12 y 15 pisos ubicados en la subzona Nib de las NTCS-2004 y estructurados con 
base en marcos ductiles con contraventeo concentrico para Q=4 utilizando: (a) el metodo 
convencional del cuerpo principal y (b) la metodologia establecida en el Apendice Normativo A 
(ANA). 



Dado lo extenso del estudio presentado en Tapia (2005), en las siguientes secciones 
se presentaran solo algunos de los aspectos mas relevantes de la comparacion de los disenos 
obtenidos para Q=4 para el edificio de 15 niveles, que son los mas interesantes para fines 
ilustrativos. 



2. EDIFICIO EN ESTUDIO Y METODOLOGIA DE DISENO 
2.1 Generalidades 

El edificio en estudio se supuso situado en la zona 1 1 lb de la nueva subzonificacion de 
las NTCS-2004, para asociarlo a las mayores demandas sismicas, con un factor de 
comportamiento sismico Q= 4 (figura 2). Adicionalmente, se busco que su diseno final 
cumpliera al limite los requisitos por deformacion lateral impuestos en cada uno de los 
planteamientos propuestos en el codigo, considerando que los elementos no estructurales se 
encuentran desligados de la estructura. Asimismo, se respeto lo sugerido en las Normas 
Tecnicas sobre Criterios y Acciones para el Diseno Estructural de Edificaciones (NTCA-2004 
2004) para la determinacion de las cargas actuantes sobre los modelos. 





a) Cuerpo principal b) Apendice A 

Figura 2 Conceptualizacion del: a) factor de comportamiento sismico Q y, b) factor de 
reduccion de las fuerzas sismicas con fines de diseno Q\ de los reglamentos de 
diseno sismico de Mexico, considerando que, por simplicidad, se supone que el 
sistema inelastico equivalente tiene un comportamiento elastoplastico perfecto. 



Se supuso que los elementos estructurales de los edificios serian con perfiles 
laminados de seccion I y para observarque no se presentara una direccion global mas debil, se 
alterno la direccion principal de las columnas como se ilustra en la figura 3. Ademas, se vario la 
orientacion de las trabes secundarias para lograr una mejor distribucion de la carga sobre las 
columnas. El sistema de piso se supuso con Losacero, cuyo diseno se realizo, conforme a las 
necesidades especificas del proyecto, con ayuda de los catalogos del proveedor. 
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a) Ubicacion de contravientos en planta tipo b) Isometrico del modelo de 15 niveles 

Figura 3 Edificio en estudio. 



El acero estructural se considero tipo A-36, en trabes y contravientos. Las columnas se 
modelaron con acero tipo A-50, variando la seccion transversal de las mismas cada tres 
niveles. Los elementos se disenaron con la herramienta de diseno del paquete de analisis 
estructural empleado (ETABS) considerando los lineamientos del manual Load and Resistance 
Factor Design (LRFD-96 1996), pero se reviso cumplieran las NTCM-2004 (2004) de Mexico. 



2.2 Excentricidad accidental 

Conforme a las nuevas NTCS-2004, en el proceso de diseno fue necesario realizar 
analisis adicionales para considerar el requisito de excentricidad accidental que se plantea por 
primera vez en el codigo, en el cual se propone trasladar transversalmente en 0.1b las 
fuerzas sismicas resultantes en cada direccion de analisis, el mismo signo en todos los niveles. 
Sin embargo, la redaccion de este analisis adicional no aclara si la revision solicita mover la 
fuerza sismica 0.1b en un eje, mientras que en otro permanezca sobre su origen (figura 4) o si 
es necesario mover la carga en 0.1b en ambos ejes (figura 5). Por otra parte, el parrafo del 
requisito citado no aclara si es posible que dado que se esta realizando un analisis dinamico, 
se considere el efecto de la torsion accidental desplazando la masa en 0.1b, con lo que se 
conseguiria la obtencion de los elementos mecanicos, aunque conllevaria el cambio de periodo 
fundamental del modelo T. 

Por ello, se estudiaron los modelos sometiendolos a un primer analisis dinamico en el 
que se desplazo la masa en 0.1b en ambos ejes, situandola en la esquina 1 de la figura 5. En 
el segundo analisis se aplicaron los cortantes dinamicos del analisis con el metodo 
convencional sin excentricidad como cargas en un analisis estatico en la misma esquina 1. En 
la figura 6 se presentan las distorsiones de entrepiso obtenidas para cada nivel, donde se 
aprecia que la direccion Y, la distorsion obtenida con el analisis estatico es mayor. Como se 
aprecia, el diseno de los edificios se hizo de manera que esta distorsion fuera muy cercana al 
limite en el nivel mas solicitado. El analisis dinamico origino menores desplazamientos, 
probablemente debido al acoplamiento de los modos. 
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Figura 6 Comparacion de las distorsiones de diseno obtenidas de los analisis con 

excentricidad accidental 



2.3 Diseno conforme al metodo convencional del cuerpo principal 

Con base en lo expuesto anteriormente, en la tabla 1 se resumen algunos de los 
resultados obtenidos en el analisis dinamico con el metodo convencional, el estatico con 
excentricidad accidental en la esquina 1 de la figura 5 y el analisis estatico sin excentricidad, el 
cual es necesario para verificar la regularidad de la estructura. 

Como se comento, los modelos se supusieron ubicados en la subzona 1 1 lb. En la figura 
7 se presentan los espectros elasticos de todas las subzonas propuestas en la nueva version 
de las NTC-2004 para su cuerpo principal, donde se aprecia que la zona lllb es la que asocia 
las mayores demandas sismicas y que el primer periodo fundamental de la estructura (T) 
analizada se encuentra en la meseta del espectro. 
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Tabla 1 Resultados de los analisis de modelos considerando el metodo convencional 



Tipo de analisis 


ler. Periodo 

(s) 


2o. Periodo 

(s) 


3er. Periodo 

(s) 


Distorsidn 

maxima 


Dinamieo sin excentricidad 


1.272 


1.216 


1.159 


0.01034 


Translacion X 


Translacion Y 


Translacion X 


Estatico sin excentricidad 


— 


— 


— 


0.01189 


Estatico con excentricidad 


— 


— 


— 


0.01195 




— Zonal -e- Zona II Pferiodo(T) 

-^ZonalHa -^ZbnaHIb 
^Zbnalllc Zona Did 




ni i t — ^ 

Distortion del analisis dinamieo 

Distorsidn del analisis estatico 



Figura 7 Espectros elasticos, zona III de las NTCS- Figura 8 Distorsiones de diseno 
2004 



En la figura 8 se comparan graficamente las distorsiones obtenidas en el analisis 
estatico con excentricidad accidental y en el analisis dinamieo. La distorsion maxima se 
presenta en el treceavo nivel. En la tabla 2 se resumen los perfiles de diseno obtenidos para 
los principales elementos estructurales. Adicionalmente, en la figura 9 se resumen los 
resultados de diseno obtenidos de la aplicacion de la formula de interaccion para la revision de 
las columnas (fig 9a), la revision de las trabes principales y secundarias (fig 9b) y, por ultimo, la 
revision de los contravientos (fig 9c). 

2.4 Diseno conforme ai Apendice A 

El Apendice Normativo A (ANA) propone un planteamiento que considera el periodo 
fundamental del suelo T s directamente en el diseno sismico del modelo estudiado. Conforme a 
este criterio, el modelo debe cumplir dos revisiones por deformacion: para limitacion de danos a 
elementos no estructurales y la otra por seguridad contra colapso. 

2.4.1 Revision de limitacion de danos a elementos no estructurales 



Esta revision debe cumplir la ecuacion 1, por considerar que el modelo propuesto no 
tendra elementos incapaces de soportar deformaciones. En ella, Sanai es la distorsion elastica 
obtenida del analisis realizado, Q’ es el factor de reduccion de fuerzas sismicas con fines de 
diseno y R es el factor por sobrerresistencia (figura 10). 



Q'R 



< 0.004 
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( 1 ) 
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Tabla 2 Secciones que satisface la revision de seguridad contra colapso, diseno conforme al 



metodo convencional del cuerpo principal de las NTCS-2004 



Nivel 


Columna 


Trabe principal 


Trabe secundaria 


Contraviento 


1 


W24"x370 


W 18"x86 


W 18"x65 


OR 10"x10"x5/8" 


2 


W24"x370 


W 18"x86 


W 18"x65 


OR 10"x10"x5/8" 


3 


W24"x370 


W 18"x86 


W 18"x65 


OR 10"x10"x5/8" 


4 


W24"x207 


W 18"x86 


W 18"x65 


OR 10"x10"x5/8" 


5 


W24"x207 


W 18"x86 


W 18"x65 


OR 10"x10"x5/8" 


6 


W24"x207 


W 18"x86 


W 18"x65 


OR 10"x10"x1/2" 


7 


W24"x176 


W 18"x86 


W 18"x65 


OR 10"x10"x1/2" 


8 


W24"x176 


W 18"x86 


W 18"x65 


OR 10"x10"x1/2" 


9 


W24"x176 


W 18"x86 


W 18"x65 


OR 10"x10"x1/2" 


10 


W24"x131 


W 18"x86 


W 18"x65 


OR 10"x10"x1/2" 


11 


W24"x131 


W 18"x86 


W 18"x65 


OR 10"x10"x3/8" 


12 


W24"x131 


W 18"x86 


W 18"x65 


OR 10"x10"x3/8" 


13 


W24"x94 


W 18"x86 


W 18"x65 


OR 10"x10"x3/8" 


14 


W24"x94 


W 18"x86 


W 18"x65 


OR 10"x10"x3/8" 


15 


W24"x94 


W 18"x86 


W 18"x65 


OR 10"x10"x3/8" 




0 1400 

C 

<b 

1 1200 
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relation 



a) Columnas b) Vigas c) Contravientos 

Figura 9 Resumen de las relaciones de esfuerzo de diseno de los principales elementos 
estructurales, diseno conforme al metodo convencional del cuerpo principal de las 
NTCS-2004 



Estos factores tienen cierta dependencia del periodo fundamental del suelo T s , ya que 
su comportamiento involucra a los periodos T a y T b , que se determinan a partir de este (NTCS- 
2004, 2004). El factor de reduccion por sobrerresistencia ( R ) es constante a partir del periodo 
T a , mientras que el factor de comportamiento sismico Q permanece constante. En la figura 10 
se ilustra la forma en que varian estos factores con el periodo (T), considerando T s - 2 s y Q= 4. 
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— ' — Q' Penodo (T) 

R 

Q’R/7 

Figura 10 Comportamiento de los factores Q’y R para T s - 2.0 s y Q- 4 

Despejando de la ecuacion 1, es posible obtener una expresion que debe cumplir la 
distorsion elastica (ec. 2) mediante la cual se satisface la revision de este primer requisito del 
ANA. 



& ANA\ ~ 



0.0280 

Q'R 



( 2 ) 



Ahora, es posible demostrar que la ecuacion 3 describe la distorsion elastica maxima 
que debe tener una estructura en el intervalo T a <T<T b , que resulta el caso de la estructura 
estudiada. En ella, se ha sustituido el valor del factor de comportamiento sismico Q y se ha 
supuesto que el factor de reduccion por sobrerresistencia es igual a R= 2: 
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ANA\ 



0.0140 k 
k +Q — 1 



( 3 ) 



Se puede concluir entonces que para que un modelo satisfaga esta revision del ANA, la 
distorsion elastica debe cumplir con la ecuacion 3, que depende unicamente del periodo 
fundamental del suelo T s y del factor de comportamiento sismico Q elegido. El coeficiente k 
toma valores distintos a k- 0.35 solo para terrenos cuyo periodo dominante sea menor a 
7“ s <1.65 s (NTCS-2004 2004). En caso que el periodo del terreno sea mayor a 7“ s >1.65 s, la 
ecuacion anterior dependera solo del factor de comportamiento sismico Q. 

Por otra parte, la revision del metodo convencional del cuerpo principal de las NTCS- 
2004 propone que la rigidez lateral debe cumplir la ecuacion 4: 

8 mc Q< 0.012 (4) 



Para el intervalo estudiado (T a <T<T b ), el factor de reduccion de fuerzas sismicas con 
fines de diseno Q’ es igual al factor de comportamiento sismico Q propuesto inicialmente. Asi, 
buscando una proporcionalidad entre las distorsiones elasticas maximas en ambos 
razonamientos, para el modelado de una misma estructura que tenga asociado el mismo factor 
de comportamiento sismico Q, se propone sustituir la ec. 4, en la ec. 3 como se muestra en la 
ecuacion 5. 
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8 ANAl — 



0.0140 k 

, 0.012 

k + 

8 mc 



(5) 



Con base en lo anterior, se puede demostrar que siempre que la distorsion elastica 
maxima obtenida en el metodo conventional S M c, de una estructura cuyo primer periodo 
fundamental T este en el intervalo T a <T<T b , satisfaga la ecuacion 5, la revision por rigidez 
lateral para la limitacion de dano a elementos no estructurales propuesta en el ANA sera 
satisfecha por la estructura para ese mismo factor de comportamiento sismico Q, sin importar 
su valor e independientemente del periodo dominante del terreno siempre que este sea mayor 
a T s >1.65 s. 

Para fundamentar el razonamiento anterior, se propuso analizar la estructura aplicando 
el espectro de diseno generado conforme al procedimiento planteado en el Apendice A, para 
un periodo dominante del suelo igual a T s - 2. s, con un factor de comportamiento sismico igual 
a Q= 4. 



La distorsion elastica resulto ser igual a S ANA1 - 0.002259, conforme a los razonamientos 
aplicados para obtener la figura 10 (ver NTCS-2004 2004), es posible demostrar que el factor 
de reduccion de fuerzas sismicas con fines de diseno es Q - 6.07 y que el factor de reduccion 
por sobrerresistencia es R= 2.0, por lo que la distorsion elastica para el estado limite de 
servicio es, aplicando la ecuacion 1, S’ ANM - 0.003918. En otras palabras, en este caso la 
misma configuracion que satisface la revision para el estado limite de prevencion de colapso 
del metodo conventional del cuerpo principal de las NTCS-2004 tambien satisface esta revision 
del estado limite de servicio del Apendice A, resultado que ya se esperaba fundamentado en 
los razonamientos analiticos anteriormente planteados. 

2.4.2 Revision de seguridad contra colapso 

La revision descrita anteriormente, no es la unica que se estipula en el Apendice A, ya 
que se solicita se realice una revision de seguridad contra colapso, que representa el primer 
esfuerzo en codigos mexicanos en hacer un diseno sismico por desempeno transparente ante 
los ojos de todo el mundo. En ella, se debe cumplir la ecuacion 6 por tratarse de marcos de 
acero con contravientos concentricos, donde S ANA2 es la distorsion elastica obtenida del 
analisis. 

S ANA2 QR < 0.015 (6) 

Buscando encontrar una relacion entre las dos estrategias de diseno permitidas en las 
NTCS-200, a partir de igualar los limites de distorsion permitidos para la revision del estado 
limite de colapso del metodo conventional (ec. 4) y del Apendice A (ec. 6) se tiene que: 

8 an A2 < 0-015(2 _ 1-25 ^ 

S MC “0.012 QR~ R 

La ecuacion 7 representa la proporcion entre las distorsiones elasticas entre ambos 
metodos para el intervalo de periodos considerado, despues de haber sido manipulada 
algebraica y aritmeticamente. De su estudio se deduce que, en el intervalo de periodos 
considerado, para que un modelo que satisface el limite de distorsion para el estado de 
prevencion de colapso establecido para el metodo conventional del cuerpo principal de las 
NTCS-2004 tambien satisfaga la restriccion de seguridad contra colapso del Apendice A, tal 
que S ANA2 /S M c- 1, es necesario que el factor de reduccion por sobrerresistencia R sea menor a 
1.25, lo cual no es posible conforme a las NTCS-2004, dado que el valor minimo establecido 
para el factor de sobrerresistencia es R= 2. Por lo tanto, en este caso la revision de seguridad 
contra colapso rige los disenos por deformacion lateral conforme al Apendice A. 
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Para ilustrar lo expuesto anteriormente, a la estructura de 15 niveles estudiada con el 
metodo convencional del cuerpo principal se le aplico el planteamiento propuesto en el ANA 
para un periodo dominante igual a T s - 2.0 s. La distorsion elastica maxima equivalente resulto 
ser S A na 2 - 0.002259, la que se debe multiplicar por el factor QR, que resulta ser QR= 8 para 
T a <T<T b . Sustituyendo en la ecuacion 6, se demuestra que la distorsion inelastica es igual a 
S’ana2- 0.0181 > Sperm- 0.015, por lo que no se satisface la desigualdad y, por tanto, esto nos 
ilustra que el diseno realizado anteriormente para el cuerpo principal no cumple con los limites 
de distorsion establecidos para el Apendice A, por lo que estrictamente hablando, ambos 
planteamientos de diseno no son totalmente equivalentes. 

Fue entonces necesario realizar un nuevo diseno con la misma configuracion 
geometrica. El diseno de la estructura se hizo de tal forma que la rigidez lateral satisficiera al 
limite los requisitos por deformacion establecidos en el Apendice A ( <0.015), por lo que no 
debe resultar extrano que las secciones transversales de los perfiles estructurales resultantes 
sean mas grandes. Los primeros tres periodos dominantes de la estructura resultaron ser de 
translacion: Ti= 1.193s, T 2 - 1.145 s y T 3 = 1.060 s. La distorsion elastica maxima se presento 
en el treceavo nivel y fue S ANA2 = 0.00186; tal que al multiplicarla por el factor QR= 8, dado que 
el primer periodo dominante de la estructura esta en el intervalo T a <T<T b , la distorsion 
inelastica es (figura 2b) S’ ANA2 =0.0'\49<Sp erm =0.0 / \5, que satisface la revision. 

Tabla 3 Secciones que satisface la revision de seguridad contra colapso, diseno conforme al 



Apendice A de las NTCS-2004 



Nivel 


Columna 


Trabe principal 


Trabe secundaria 


Contraviento 


1 


W24"x370 


W 18"x97 


W 18"x76 


OR 12"x12"x5/8" 


2 


W24"x370 


W 18"x97 


W 18"x76 


OR 12"x12"x5/8" 


3 


W24"x370 


W 18"x97 


W 18"x76 


OR 12"x12"x5/8" 


4 


W24"x250 


W 18"x97 


W 18"x76 


OR 12"x12"x5/8" 


5 


W24"x250 


W 18"x97 


W 18"x76 


OR 12"x12"x5/8" 


6 


W24"x250 


W 18"x97 


W 18"x76 


OR 12"x12"x1/2" 


7 


W24"x192 


W 18"x97 


W 18"x76 


OR 12"x12"x1/2" 


8 


W24"x192 


W 18"x97 


W 18"x76 


OR 12"x12"x1/2" 


9 


W24"x192 


W 18"x97 


W 18"x76 


OR 12"x12"x1/2" 


10 


W24"x146 


W 18"x97 


W 18"x76 


OR 12"x12"x1/2" 


11 


W24"x146 


W 18"x97 


W 18"x76 


OR 12"x12"x3/8" 


12 


W24"x146 


W 18"x97 


W 18"x76 


OR 12"x12"x3/8" 


13 


W24"x1 17 


W 18"x97 


W 18"x76 


OR 12"x12"x3/8" 


14 


W24"x1 17 


W 18"x97 


W 18"x76 


OR 12"x12"x3/8" 


15 


W24"x1 17 


W 18"x97 


W 18"x76 


OR 12"x12"x3/8" 



En la tabla 3 se resumen los perfiles de diseno obtenidos para los principales 
elementos estructurales. Adicionalmente, en la figura 11 se resumen los resultados de diseno 
obtenidos de la aplicacion de la formula de interaccion para la revision de las columnas (fig 
11a), la revision de las trabes principales y secundarias (fig 11b) y, por ultimo, la revision de los 
contravientos (fig 11c). Cabe senalarque los elementos estructurales obtenidos de este diseno 
conforme al Apendice A estan sometidos a menores relaciones de esfuerzos actuantes contra 
esfuerzos resistentes, en comparacion con los disenos conforme al metodo convencional del 
cuerpo principal reportados anteriormente (figura 9), lo que no significa que las secciones 
propuestas estan sobradas, sino que mediante ellas se aporto la rigidez necesaria para salvar 
el requisito de deformacion respetando la configuracion elegida. 
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a) Columnas b) Vigas c) Contravientos 

Figurall Resumen de las relaciones de esfuerzo de diseno de los principales elementos 
estructurales, modelo disenado conforme al Apendice A de las NTCS-2004 



3. ANALISIS ESTATICOS NO LINEALES 



Se presentan los resultados de analisis estaticos no lineales ante cargas laterales 
mondtonas y crecientes (tambien conocidos coloquialmente con el anglicismo “pushover”) con 
el programa DRAIN-2DX (Prakash et al. 1993), para estudiar el comportamiento de los modelos 
incluyendo los efectos P-A y sin incluirlos en contravientos y columnas. 



Se utilizaron patrones de carga lateral asociadas al mecanismo de columna fuerte - 
viga debil - contraviento mas debil, que se denomina como metodo convencional por ser en 
teoria el asociado a un diseno conforme a las NTCS-2004, asi como patrones de carga lateral 
identico al utilizado en analisis al limite para un mecanismo de onceavo piso debil (no 
mostrados), por ser el que rigid el criterio con los analisis al limite (Tapia 2005), como se 
ilustran esquematicamente en la figura 12. 



Los analisis estaticos no lineales se realizaron modelando a los marcos tanto de 
manera independiente (figura 12) como considerando que todos los marcos en una direccion 
de analisis trabajan en conjunto, para lo cual fue necesario modelar la accion del diafragma 
rigido por medio de elementos elasticos axialmente rigidos, como se ilustra en la figura 13. 
Ademas, para observar el impacto de la asimetria de los marcos centrales (B, D, 2 y 4, figura 
13) se realizaron analisis a los que se cambio el signo de la aplicacion de la carga que en el eje 
X para observar su diferencia. Los detalles de todos estos analisis se reportan en Tapia (2005). 
Por fines de espacio y claridad, se presentan los resultados de los modelos mas completos, 
que son aquellos donde se considera la accion e interaccion de todos los marcos ante carga 
lateral (figura 13) y donde se incluyeron los efectos de segundo orden. 








12 



Edgar Tapia Hernandez y Arturo Tena Colunga 





b) Mecanismo supuesto en las NTCS-2004 



Figura 12 Mecanismo de falla supuesto en los analisis estaticos no lineales ante carga 
monotona creciente 




Marco A 6 1 Marco B 6 2 Marco C 6 3 Marco D 6 4 Marco E 6 5 

Figura 13 Modelo con distribucion de las fuerzas laterales entre los marcos resistentes en la 
direccion de analisis 



Como se reporta con mucho mayor detalle en Tapia (2005), los resultados obtenidos 
considerando marcos independientes son similares, pero no iguales, a los que se obtienen 
considerando la interaccion de todos marcos, dado que el no modelar la interaccion que en 
ellos impone la rigidez del diafragma no fomenta la redistribucion de fuerzas, que es benefica 
ante carga lateral. Como se ilustra en la figura 14, obtenida para los modelos disenados 
conforme al metodo convencional del cuerpo principal de las NTCS-2004, la diferencia es 
mayor para marcos exteriores (marcos A y E) que para marcos interiores (marco C). El cortante 
basal que resisten en conjunto todos los marcos varla ligeramente entre ambas modelaciones, 
resultando mayor en un 8.7% cuando se analizan los marcos en conjunto (figura 13) a cuando 
se calculan de manera independiente y posteriormente se suman (Tapia 2005). 
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Figura 14 Comparacion de los cortantes basales obtenidos para los modelos disenados 
conforme al cuerpo principal de las NTCS-2004 



Por simplicidad, en las siguientes secciones se identificara como MC15 a los resultados 
de los analisis para el edificio disenado conforme al metodo convencional del cuerpo principal 
de las NTCS-2004 y ANA15 a los resultados del edificio disenado conforme al Apendice A del 
mismo cuerpo normativo. 



2.4.3 Mapeo de fluencias 

En las figuras 15 a 19 se presentan el mapeo final de fluencias de los analisis 
realizados para los modelos MCI 5 y ANA15. En estas figuras se ilustran con una esfera la 
magnitud y ubicacion de las rotaciones inelasticas presentadas en el paso en que se presenta 
el colapso. Las rotaciones graficadas se normalizaron respeto al mayor giro (o extension o 
contraccion en el caso de los contravientos) presentado en cada tipo de elemento, es decir, la 
esfera de mayor tamano representa la articulacion con un diametro igual a uno en las trabes, 
en las columnas y en los contravientos (extension o contraccion), tal que el resto de las 
rotaciones inelasticas tendran un diametro menor, aunque no implica que la magnitud de las 
rotaciones en esos tres grupos de los elementos sea comparable. Esta presentacion obedece a 
que la magnitud de algunas rotaciones en trabes y columnas resultaba tan grande en 
comparacion con sus analogas, que obligaba a que el diametro del resto de las esferas 
resultara poco apreciable. 

En la figura 15 se presenta el mapeo final de fluencias del analisis del modelo MCI 5 en 
la direccion global X con el efecto P-A cuando la carga lateral se aplica de derecha a izquierda 
(D-l). En general, en vez de presentarse el mecanismo de columna fuerte - viga debil - 
contraviento mas debil supuesto en un diseno conforme a las NTCS-2004, se identifica un 
mecanismo de colapso mixto y poco ductil, favoreciendo entrepisos debiles entre el quinto y el 
sexto nivel. Al parecer, despues de la cedencia de los contravientos, la redistribucion de carga 
no puede ser absorbida por las columnas, las que despues de articularse en ambos extremos 
ocasionan el colapso. 
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Figura 15 Mapa de fluencias del modelo MCI 5 con P-A, eje X con carga de D-l 



En la figura 16 se presenta el mapeo de fluencias de los mismos marcos del modelo 
MCI 5 con la aplicacion de la carga de izquierda a derecha (l-D). Se observa que la ubicacion 
de las articulaciones no resulta simetrica a los resultados mostrados anteriormente. Cabe 
senalar que existe un mayor numero de trabes articuladas y que el colapso de la estructura no 
se puede asociar al mecanismo formado entre el quinto y sexto nivel como en el caso en el que 
la carga se aplica con signo positivo. 
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Figura 16 Mapa de fluencias del modelo MC15 con P-A, eje X con cargas de l-D. 



En la figura 17 se presenta el mapeo de fluencias del modelo MC15 en la direccion 
global Y cuando la carga se aplica de derecha a izquierda (D-l), que es la que rige. En general, 
presenta una distribucion de rotaciones muy similar al analisis de la direccion global X, con un 
mecanismo de colapso porfalla local en la distribucion de cargas en las columnas. 
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Figura 17 Mapa defluencias del modelo MCI 5 con P-A, eje Y, con cargas de D-l. 



En la figura 18 se muestran los mapeos de fluencia del modelo ANA15 en la direccion 
global X considerando los efectos de segundo orden, aplicando la carga de derecha a 
izquierda. Se aprecia que sigue rigiendo un mecanismo de colapso por falla local en las 
columnas del lado derecho de la ultima crujia contraventeada del marco central, aunque su 
ubicacion ha variado, lo que parece estar ligado a la direccion de aplicacion de la carga y la 
asimetria del modelo. Cabe senalar que las articulaciones de los marcos internos B y D 
alcanzan rotaciones mucho menores al marco central y marcos perimetrales, que a su vez se 
articulan en los niveles superiores entre las crujfas no contraventeadas. 
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Figura 18 Mapas de fluencia del modelo ANA15 con P-A, eje X, con cargas de D-l. 



En la figura 19 se presenta el mapeo de fluencias del modelo ANA15 en la direccion Y 
con el efecto P-A, aplicando la carga de derecha a izquierda. Como se esperaba, la ubicacion y 
la proporcion de las rotaciones plasticas son similares a la direccion X, presentadas 
anteriormente. 
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Figura 19 Mapa de fluencias del modelo ANA15 con P-A, eje Y, con cargas de D-l. 



De la comparacion de los mapeos de fluencia de los modelos MCI 5 (figuras 15 y 17) 
con los respectivos modelos ANA15 (figuras 18 y 19), se aprecia que aunque ambos sugieren 
un mecanismo de colapso mixto poco ductil con fluencias y articulaciones plasticas en 
columnas mas que en vigas y contravientos, los mecanismos de los disenos conforme al 
Apendice A (ANA15) presentan un mayor numero de fluencias y articulaciones plasticas que 
los asociados al diseno conforme al cuerpo principal (MCI 5) y, en este sentido, disipan mayor 
energla y, por ello, son ligeramente superior. En cualquier caso, lo realmente alarmante es que 
con ninguno de los dos procedimientos de diseno se logra obtener el mecanismo de falla 
supuesto en un diseno conforme a las NTCS-2004: columna fuerte - viga debil - contraviento 
mas debil. 



En la figura 20 se comparan las distorsiones de entrepiso calculadas para ambas 
direcciones (X y Y) para los modelos disenados con el metodo convencional y el disenado con 
el Apendice A, cuando se incluye el efecto P-A. En general, el modelo disenado con el metodo 
convencional reporta distorsiones mayores que las del modelo analizado con los criterios del 
Apendice A; sin embargo, en los niveles cercanos al sexto piso en la direccion X, el modelo 
disenado con el metodo convencional alcanza distorsiones sensiblemente mayores. Las 
distorsiones de la direccion global Y son menores porque es la direccion mas rigida, aunque en 
todos los casos las distorsiones son menores a la distorsion maxima permisible perm = 0.015 
para marcos contraventeados de acero que se establece en el Apendice A de las NTCS-2004. 

En la figura 21 se comparan los cortantes por entrepiso normalizados con respecto al 
peso total de la estructura para ambos disenos en ambas direcciones. Se aprecia que los 
modelos disenados con el Apendice A soportan cortantes de mayor magnitud en todos los 
niveles, particularmente en los niveles inferiores. 

En las figuras 22 y 23 se comparan para los modelos MC15 y ANA15 las curvas 
cortante normalizado contra distorsion por nivel en las direcciones X y Y y la curva cortante 
basal normalizado contra la distorsion total de la estructura calculada con la altura total del 
edificio, denominada analisis general. Se observa de ambas figuras que los primeros siete 
niveles se comportan inelasticos, mientras que los ocho niveles superiores tiene un 
comportamiento elastico, lo que concuerda con los mapeos de fluencia presentados 
anteriormente. Sin embargo, cabe senalar que estas graficas incluyen los entrepisos de todos 
los marcos, por lo que al promediar no reflejan las articulaciones de los pisos superiores en los 
marcos perimetrales y se concentran en las articulaciones de mayor magnitud del marco central 
en los pisos cercanos al sexto nivel y que, la curva global cortante-distorsion (“analisis general “) 
no parece ser una buena representacion de las curvas de los niveles inferiores, pues en el 
promedio pesado domina el comportamiento de los entrepisos superiores (niveles 6 al 15). 
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Figura 21 Cortantes normalizados de entrepiso para los modelos MC15 y ANA15 
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Figura 22 Comparacion de las curvas cortante-distorsion, modelos MCI 5 y ANA15, Dir X 
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Figura 23 Comparacion de las curvas cortante-distorsion, modelos MCI 5 y ANA15, Dir Y 
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Como ya se habia ilustrado y comentado, tambien queda claro de la observacion de las 
curvas presentadas en las figuras 22 y 23 que el modelo disenado con el Apendice A desarrolla 
mayores cortantes de entrepiso, pero menores distorsiones y capacidad de deformacion ultima 
en los entrepisos con comportamiento inelastico (inferiores). 

En la tabla 4 se presenta el calculo de la sobrerresistencia maxima que, en teoria, 
desarrollarlan los disenos realizados conforme al metodo convencional del cuerpo principal 
(MCI 5) y el Apendice A (ANA15) de las NTCS-2004. Como se aprecia, en ambos casos esta 
resulta bastante superior al valor R = 2 que propone el Apendice A de las NTCS-2004 para 
estructuras en este intervalo de periodos fundamentales y que, de hecho, fue el factor utilizado 
en el diseno del modelo ANA15. Este resultado era esperado, dado que la propuesta del factor 
de sobrerresistencia R de las NTCS-2004 se calibro inicialmente para obtener resistencias 
similares a las demandadas por el reglamento anterior (NTCS-96) para el diseno de marcos 
ductiles disenados con Q=3 y Q=4 (Ordaz et al. 2003) y, de hecho, el valor R= 2 esta basado en 
algunos estudios hechos para marcos ductiles de concreto reforzado empleando analisis 
estaticos no lineales ante carga monotona creciente (“pushover”). Por lo tanto, antes de este 
estudio no existe ningun estudio que indicara la sobrerresistencia que este sistema estructural 
puede desarrollar. 



Tabla 4 Calculo de la sobrerresistencia desarrollada por los modelos MC15 y ANA15 



Modelo 


Coeficiente sismico 
de diseno 


Cortante resistente 
ultimo 


Sobrerresistencia 
desarrollada (R) 


MC15, Dir X, P- 


0.1125 


0.5310 


4.72 


MC15, DirY, P- 


0.1125 


0.5119 


4.55 


ANA15, Dir X, P- 


0.1976 


0.6078 


3.08 


ANA15, DirY, P- 


0.1976 


0.6058 


3.07 



Como tambien era de esperarse, dados todos los factores de seguridad (un poco 
obscuros y ocultos) inherentes a un diseno conforme al metodo convencional del cuerpo 
principal de las NTCS-2004, la sobrerresistencia que podrian desarrollar los modelos MCI 5 es 
muy grande, notablemente mayor a la de los modelos disenados conforme al Apendice A. En 
este sentido, estos resultados sugieren que el nuevo planteamiento de diseno establecido en el 
Apendice A es mas racional, dado que se tiene una mejor estimacion de la resistencia que 
realmente se suministra, y que es perfectible en la medida que se pueda determinar, con 
estudios adicionales, cual debe ser el factor de sobrerresistencia R que debiera emplearse en 
el diseno de este sistema estructural. 



4. ANALISIS DINAMICOS NO LINEALES 

Los analisis dinamicos no lineales se realizaron tambien con el programa DRAIN-2DX, 
considerando inicialmente el acelerograma registrado durante el sismo del 19 de septiembre de 
1985 con una magnitud igual a M s - 8.1 del registro de la Secretaria de Comunicaciones y 
Transportes (SCT) en su componente E-W, con una duracion de un poco mas de 160 s y una 
aceleracion maxima del terreno de 0.18 g. aproximadamente, con una respuesta espectral 
maxima de 1.0 g. para un periodo dominante de dos segundos (figura 24). Los analisis se 
realizaron con el modelo de los marcos continuos unidos con elementos rigidos (figura 13), con 
un nivel de amortiguamiento viscoso equivalente igual a C= 2%. 



Ademas, los modelos se sometieron a un registro artificial que se genera a partir de un 
sismo postulado de M s - 8.2 ocurrido en la brecha de Guerrero, que es el sismo de diseno 
considerado en las NTCS-2004, utilizando funciones de transferencia empiricas promedio 
calculadas a partir de 1988 para la estacion 56 de la red acelerometrica de la Ciudad de 
Mexico. La estacion 56 se ubica en la calle de Cordoba, en la colonia Roma, con un periodo 
dominante de 2.2 segundos. El espectro de pseudoaceleraciones de este registro artificial es 
mayor al espectro calculado para el acelerograma registrado en SCT-EW (figura 24). 
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Figura 24 Acelerogramas empleados en este estudio 



En la figura 25 se presenta la definicion de los parametros que se discuten durante la 
exposicion de los resultados. En ella, Afi uencia es la distorsion de entrepiso asociada a la primera 
fluencia de cualquier elemento estructural que forma parte del entrepiso, A maxima es la distorsion 
de entrepiso maxima asociada al medio ciclo de histeresis de carga maxima, k imin es el cociente 
entre la rigidez de entrepiso secante minima inelastica, /r eHsf , ca es la rigidez elastica de entrepiso 
y k mc es la rigidez asociada al medio ciclo histeretico inelastico; ademas, de la ductilidad de 
entrepiso ju, las envolventes de distorsiones de entrepiso maximas A, envolventes de cortantes 
de entrepiso maxima (V/W T ) y los cocientes del numero de ciclos donde se presento respuesta 
inelastica MCj ne entre el numero total de ciclos a que fue sujeta la estructura ante el 
acelerograma dado MC tot . 

A . . 

maxima 




Figura 25 Definicion de parametros 

Dado que los respuestas maximas obtenidas con el acelerograma S56-EW resultaron 
similares, pero ligeramente mayores a las calculadas con el registro SCT85-EW, a continuacion 
solo se presentaran los resultados ante el acelerograma S56-EW y solo en la direccion X, dado 
que los resultados obtenidos en la direccion Y proporcionan informacion similar. Los detalles de 
todo este estudio ante ambos registros acelerometricos y en ambas direcciones se pueden 
consultar en Tapia (2005) y Tapia y Tena (2005). 
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En la figura 26 se presentan las envolventes de respuesta maximas de la direccion X 
de ambos modelos ante la excitacion del acelerograma artificial S56-EW. Como se observa, se 
presenta una respuesta no lineal reducida con demandas de ductilidad menores a =1.64 en 
el treceavo nivel del modelo disenado con el metodo convencional del cuerpo principal de las 
NTCS-2004 (modelo MC15). Los pisos inferiores se comportaron sensiblemente elasticos, con 
excepcion del primer nivel del modelo MC15. Ademas, se presentaron repentinos cambios de 
ductilidad en entrepisos continuos en el modelo MC15. 





0.8 1.0 1.6 
DISTORSION, A 





2.00.00 0.00 0.16 0.24 0.32 0.40 



fluencia 



maxima 



v/w 




MC. /TOTAL MC DUCTILIDAD K . . /K . K .. /K . 

me minimi elastica promedio elastica 

Figura 26 Envolventes de respuesta maxima de los modelos MC15 y ANA15X en direccion X 
cuando son sujetos al acelerograma S56-EW 



El modelo ANA15 desarrolla proporcionalmente un mayor numero de ciclos inelasticos 
(MC ine /MC total) que el modelo MCI 5, incluso cercanas al 10%, entre el noveno al doceavo 
nivel. Las rigideces secantes minimas inelasticas de entrepiso k min resultaron practicamente del 
100% de su rigidez elastica k e ia S ti Ca , en ambos casos en todos los entrepisos. Cabe senalar que 
la distorsion asociada a la primera fluencia A f i uencia de los entrepisos que fluyen ( >1), 
practicamente en todos los niveles superiores al sexto nivel, se encuentran cercanas al 0.4%, 
que es el limite actualmente establecido en el Apendice A de las NTCS-2004 para revisar a la 
estructura ante el sismo de servicio. 

En las figuras 27 y 28 se presentan las envolventes de fluencia de ambos modelos en 
la direccion X. Para el modelo MCI 5, la rotacion maxima en vigas fue de 0.0045 rad y en 
columnas de 0.0003 rad, mientras que el acortamiento maximo en contravientos fue de 2.41 
cm. Para el modelo ANA15, la rotacion maxima en vigas fue de 0.0015 rad y en columnas de 
0.0001 rad, mientras que el acortamiento maximo en contravientos fue de 1.06 cm. Se aprecia 
que se presentan acortamientos pequenos que sugieren el pandeo incipiente en los 
contravientos de todos los niveles y que estos son mas importantes en los dos ultimos niveles. 
Ademas, tambien se presenta la fluencia incipiente de algunas trabes en los pisos superiores. 
Cabe senalar que en el diseno conforme al metodo convencional (MCI 5, figura 27) se 
presentan fluencias incipientes en ambos extremos de las columnas del marco central en el 
cuarto nivel, ademas de la fluencia de las columnas en la base del edificio de ese mismo 
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marco, mientras que en el modelo ANA15 solo se presenta una fluencia incipiente en la 
columna exterior del primer nivel del marco perimetral (figura 28). 
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Figura 27 Envolventes de fluencia del modelo MCI 5, dir. X, ante el acelerograma S56-EW 
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Figura 28 Envolventes de fluencia del modelo ANA15, dir. X, ante el acelerograma S56-EW 



Los lazos de histeresis de entrepiso de ambos modelos se presentan en las figuras 29 
y 30. Como se aprecia, resulta apenas perceptible la no linealidad en los niveles superiores, 
principalmente en el modelo disenado con el metodo convencional (MCI 5) un analisis modal 
espectral, asi como el incremento de ruido en la grafica que muestra el comportamiento global 
del edificio con la distorsion promedio y el cortante basal normalizado, que debiera ser 
sensiblemente lineal, conforme a los resultados mostrados, que es fruto de errores por 
truncacion en su calculo. 
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Figura 29 Curvas de histeresis de entrepiso del modelo MCI 5 en direccion X cuando es sujeto 
al acelerograma S56-EW 
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Figura 30 Curvas de histeresis de entrepiso del modelo ANA15 en direccion X cuando es 
sujeto al acelerograma S56-EW 



Como se aprecia, los modelos MCI 5 y ANA 15 presentaron un comportamiento no 
lineal muy reducido (en algunos entrepisos nulo) ante los acelerogramas SCT-EW y S56-EW 
que estan asociados al sismo maximo de diseno considerado en las NTCS-2004. Al respecto, 
ambos modelos se disenaron considerando que el periodo fundamental de vibracion del suelo 
blando de desplante es T s - 2 s, tal que la aplicacion de ambos acelerogramas representa las 
peores condiciones de excitacion para el sitio. La razon del por que se obtuvieron estas 
respuestas maximas tan reducidas se explica con apoyo de la figura 31, donde se presenta el 
espectro de respuesta elastico del acelerograma S56-EW, en el cual tambien se indica el 
periodo fundamental de vibracion de los modelos estudiados (Tj). Si se considera que el 
periodo fundamental de vibracion de los modelos estudiados son 27s para el modelo 
MCI 5 y Tj=1 .18s para el modelo ANA15, se observa de las lineas verticales en la figura 31 que 
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estos periodos se ubican en la rama ascendente de los espectros, pero en una zona donde no 
se presentan las maximas aceleraciones espectrales. Dado que el cortante basal resistente de 
estos edificios es relativamente alto, esto explica el por que se obtuvieron demandas 
inelasticas tan reducidas, liderando a las estructuras a un comportamiento dinamico 
practicamente elastico. 




Figura 31 Espectro de respuesta de pseudoaceleraciones del acelerograma S56-EW con la 
ubicacion de los periodos fundamentales de los modelos MCI 5 y ANA15 



Tabla 5 Sobrerresistencia demandada en los modelos ante la accion del acelerograma 
artificial S56-EW 



Modelo 


Direccion 


Cortante basal maximo normalizado 
V/W T 


Sobrerresistencia demandada 


MCI 5 


X 


0.3678 


3.27 


Y 


0.3019 


2.68 


ANA15 


X 


0.3402 


1.72 


Y 


0.2711 


1.37 



Finalmente, en la tabla 5 se presentan los cortantes basales dinamicos maximos 
(normalizados con respecto al peso de la estructura) obtenidos en los analisis realizados ante 
el acelerograma artificial S56-EW. Considerando que el comportamiento no lineal demandado 
es muy reducido, resulta claro que las sobrerresistencias demandadas en ambos modelos son 
relativamente bajas, dado que poseen reservas de resistencias aun mayores antes de ser 
llevados a un mecanismo de colapso, conforme ya se reporto anteriormente en la tabla 4 
cuando se determino su capacidad maxima a partir de analisis estaticos no lineales ante carga 
monotona creciente (pushover). Aun asi, cabe senalar que el modelo MC15 desarrollo una 
sobrerresistencia R> 2.50, que es mayor a la propuesta actualmente en el Apendice A. 



5 COMENTARIOS FINALES 

Se presento el estudio del comportamiento no lineal de dos edificios estructurados con 
base en marcos ductiles de acero estructural con contraventeo concentrico, uno de ellos 
disehado con el metodo convencional del cuerpo principal y el otro con el Apendice A de las 
NTCS-2004, ambos disenados de tal forma que cumplen al limite las restricciones por 
desplazamiento lateral de cada criterio. 

La capacidad ultima de los edificios se estudio utilizando analisis no lineales ante carga 
estatica monotona creciente (“pushover”), despreciando los efectos P - A (no mostrados) e 
incluyendo los efectos P - A (mostrados). La inclusion del efecto de segundo orden afecta de 
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manera definitiva los resultados de los analisis estaticos, obteniendose en la mayoria de los 
casos distorsiones ultimas de entrepiso menores que cuando no se incluyen, y con magnitudes 
maximas mas cercanas a los limites que proponen los reglamentos, los cuales estan en funcion 
de resultados obtenidos a partir del ensaye de modelos experimentales. Es por ello que, para el 
analisis de estructuras de acero con marcos contraventeados concentricos, se deben incluir los 
efectos P- A en estos analisis tipo pushover. 

De los resultados obtenidos con los analisis estaticos no lineales ante carga monotona 
creciente se pueden resaltar tres aspectos relevantes. 

El primero tiene que ver con el mecanismo ultimo de colapso de los modelos. En 
ambos casos el mecanismo de colapso desarrollado es de naturaleza mixta, poco ductil con 
fluencias y articulaciones plasticas importantes en columnas mas que en vigas y contravientos. 
El mecanismo de colapso de los disenos conforme al Apendice A presentan un mayor numero 
de fluencias y articulaciones plasticas que el asociado al diseno conforme al cuerpo principal y, 
en este sentido, disipan mayor energia y, por ello, es ligeramente superior. En cualquier caso, 
lo realmente alarmante es que con ninguno de los dos procedimientos de disefio se logra 
obtener el mecanismo de falla supuesto en su diseno conforme a las NTCS-2004: columna 
fuerte - viga debil - contraviento mas debil. 

El segundo tiene que ver con las curvas carga-distorsion obtenidas de los analisis. 
Dentro del grupo de gente que disena “por desempeno” es una costumbre muy generalizada 
que realicen disenos con analisis pushover viendo exclusivamente la informacion de la curva 
general cortante basal vs desplazamiento de azotea y sus variantes normalizadas. De hecho, 
este procedimiento lo establecen documentos como FEMA-273 (1997) o FEMA-368 (2001) 
para todos los sistemas estructurales, aunque se basa exclusivamente en el estudio de marcos 
a momenta. Los resultados de este estudio muestran claramente la grafica cortante basal 
normalizado contra distorsion global promedio (distorsion entre azotea y el nivel de desplante) 
parece no representar adecuadamente el comportamiento general de marcos de acero con 
contravientos concentricos, ya que a partir de ella diflcilmente puede uno intuir las grandes 
diferencias que existen en las demandas de ductilidad y la magnitud de las distorsiones de los 
distintos entrepisos, que dista mucho de ser uniforme. Ademas, es bien sabido que seria 
incorrecto comparar esta distorsion promedio que se obtiene con esta curva con las 
distorsiones ultimas que permiten los Reglamentos, que estan asociadas a las distorsiones de 
entrepiso. Para los modelos en estudio, esta grafica global no resulta una buena 
representacion del comportamiento global de este tipo de estructuracion, por lo que se propone 
estudiar para esta y otras estructuraciones las graficas distorsion contra cortante normalizado 
de cada entrepiso, para poder observar el comportamiento no lineal del modelo con mayor 
certidumbre, revisar adecuadamente los limites de distorsion que establecen los reglamentos 
de construccion y, ademas, detectar posibles deficiencias en la estructuracion con mayor 
facilidad. 



Este estudio muestra la imperiosa necesidad de confrontar siempre las curvas carga- 
distorsion obtenidas de un analisis no lineal estatico ante carga monotona creciente en funcion 
de los limites de deformacion que se han reportado experimentalmente en la literatura y de las 
capacidades ultimas de deformacion que tienen los elementos estructurales, que no son 
necesariamente las que se calculan con el software que se utiliza, dado que esto normalmente 
las modelan tomando en cuenta varias hipotesis simplificadoras. Es por ello que se presentaron 
en este trabajo algunos resultados que aun deben refinarse, dado que es comun que una 
buena cantidad de la gente (mas no todos obviamente) que esta trabajando en criterios diseno 
por desempeno utilizando analisis no lineales estaticos ante carga monotona creciente, no 
reflexionan de manera suficiente sobre estos aspectos, que son muy importantes y, por ello, 
despues reportan resultados donde se obtienen capacidades de deformacion que no 
corresponden a los sistemas estructurales que estan estudiando, por considerar que el colapso 
de la estructura esta dado con bastante certidumbre por el software de analisis no lineal que 
estan utilizando, y no que es posible obtener numericamente respuestas mayores a las que en 
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realidad la estructura puede desarrollar en funcion de la capacidad ultima de deformacion “real” 
de cada uno de sus elementos, que es funcion de sus caracteristicas propias y su detallado. 

El tercer aspecto relevante esta relacionado con el nivel de sobrerresistencia que 
pueden desarrollar estos sistemas estructurales a pesar de que su mecanismo de falla ultimo 
no sea ductil. Los factores de sobrerresistencia obtenidos para los modelos disenados 
sobrepasan notablemente el factor de sobrerresistencia R propuesto actualmente en las NTCS- 
2004. 



Se realizaron tambien analisis dinamicos no lineales de los modelos ante la accion del 
acelerograma registrado durante el sismo del 19 de septiembre de 1985 en la Secretaria de 
Comunicaciones y Transportes en su componente este-oeste (SCT85-EW) y con un 
acelerograma artificial de sitio generado para la estacion 56 ubicada en la colonia Roma (S56- 
EW), donde los efectos del sismo del 19 de septiembre de 1985 fueron devastadores. Estos 
dos acelerogramas estan asociados a los espectros de disefio que fueron empleados para el 
disefio de ambas estructuras. Los resultados de estos analisis sugieren que la respuesta de 
estas estructuras ante acelerogramas asociados al espectro de diseno seria esencialmente 
elastica, con una moderada no linealidad, debido a dos razones fundamentales: (1) su alta 
resistencia lateral y sus reservas de resistencia y (2) el desacoplamiento dinamico existente 
entre sus periodos fundamentales de vibracion y el periodo dominante del sitio, que lidera a 
que estas estructuras sean sujetas a demandas menores. 

Una informacion que resulto un poco diferente de los analisis dinamicos no lineales con 
respecto a los analisis estaticos no lineales es que en el dinamico el mapeo de las fluencias 
incipientes a bajos niveles de demandas si sugiere que los disenos liderarian a que los 
primeros elementos en responder inelasticamente serian los contravientos y despues vigas, 
aunque esto no evita el que fluyeran algunas columnas en los entrepisos intermedios o 
superiores en proporciones similares a las vigas en los disenos conforme al metodo 
convencional de las NTCS-2004. Para ello se requieren realizar estudios adicionales que 
permitan valorar si el mapeo de fluencias ante excitaciones dinamicas mayores, que sean 
realistas para el sitio en estudio, indica que el mecanismo de falla ultimo es muy distinto al 
obtenido en un analisis pushover. 

Conforme a los resultados de los analisis no lineales se puede concluir que las buenas 
noticias (entre comillas) son que los marcos de acero con contraventeo ductil de mediana a 
baja altura ubicados en la zona de terreno blando del Distrito Federal (en concreto, la zona lllb) 
y disenados conforme a las NTCS-2004 seguramente tendran un comportamiento satisfactorio 
ante el sismo de diseno, debido a las dos razones expuestas anteriormente: (a) 
desacoplamiento dinamico y, (b) su alta resistencia lateral. 

Sin embargo, desde otra optica, estos resultados denotan que la aplicacion de las 
NTCS-2004 conlleva a un disefio ineficiente, donde el mecanismo de falla ultimo supuesto en el 
diseno de columna fuerte - viga debil - contraviento debil aparentemente no se cumple. En 
este sentido, el procedimiento de diseno establecido en el Apendice A resulta ligeramente 
superior al metodo convencional del cuerpo principal de las NTCS-2004. Resulta claro 
entonces que los resultados de estos analisis dinamicos no lineales, y sobre todo los analisis 
pushover, sugieren la imperiosa necesidad de realizar una revision integral de los criterios de 
diseno de las NTCS-2004 para este tipo de estructuras, pues ante acelerogramas asociados al 
sismo de diseno se tienen reservas de resistencia muy grandes que lideran a la estructura a 
comportarse de manera casi elastica, y los analisis pushover sugieren mecanismos ultimos de 
falla diferentes e incompatibles a los que se suponen en su disefio. Actualmente los autores 
estan realizando estudios que estan abordando esta problematica y que intentan darle una 
solucion practica. 
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RESUMEN 

La simulacion numerica del comportamiento de los materiales compuestos es un 
campo en desarrollo cuya aplicacion a estructuras ha experimentado avances importantes que 
han conducido a mejoras en el refuerzo de las mismas. En relacion con este tema, el articulo 
parte de la teoria de de mezclas serie/paralelo que permite deducir el comportamiento de los 
compuestos a partir de las ecuaciones constitutivas de los materiales componentes. Con el fin 
de destacar las ventajas que tiene esta herramienta en el analisis y el diseno de estructuras 
reparadas o rehabilitadas con polimeros reforzados con fibras largas (FRP), en el articulo se 
realizan analisis numericos mediante modelos de elementos finitos de un conjunto de porticos 
con distintas configuraciones de refuerzo. Los resultados obtenidos muestran que el uso esta 
teoria en los analisis de las estructuras reforzadas con FRP contribuye a mejorar el diseno y el 
comportamiento de las estructuras. 

ABSTRACT 

The numerical simulation of the behaviour of composite materials is a developing field 
whose application to structures evidenced important advances allowing improvements in their 
strengthening. Related to this topic, this paper starts from the serial/parallel mixing theory which 
permits deducing the behaviour of the composites based on the constitutive equations of the 
component materials. In order to reveal the advantages this computational tool for the analysis 
and design of structures repaired o rehabilitated with Fiber Reinforced Polymers (FRP), 
numerical analyses using finite element models have been performed in this article to a set of 
frames with different strengthening configurations. The obtained results show that the use of 
this theory into analysis of the fiber reinforced structures contributes to improve the design and 
the behaviour of the structures. 



1. INTRODUCTION 

El uso de los materiales compuestos como materiales de refuerzo para estructuras ya 
construidas es una de las tecnologfas que estan teniendo gran aplicacion en la industria de la 
construccion debido a las ventajas que presentan. Entre estas ventajas estan las altas 
relaciones resistencia-peso y rigidez-peso, las cuales mejoran el comportamiento de la 
estructura existente sin alterar su configuracion geometrica. Asimismo, son materiales livianos 
que no demandan cambios en la distribucion del sistema estructural o en la cimentacion. 
Ademas para el caso de construcciones sometidas a ambientes especiales, presentan alta 
resistencia a la corrosion. 

No obstante, para optimizar el diseno estructural de los materiales compuestos en las 
estructuras, es necesario el analisis del comportamiento de los materiales compuestos donde 
se identifique la forma de participacion de los diferentes componentes del compuesto y su 
interaccion con otros materiales como el hormigon armado. De igual forma, se requiere el 
analisis del comportamiento global de la estructura, en el que se establezca la incidencia de los 
materiales compuestos como parte integral de los elementos estructurales reforzados. Un 
procedimiento eficiente para realizar estos analisis es la simulacion numerica con elementos 
finitos. 
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Teniendo en cuenta que la simulacion numerica en el analisis de los materiales compuestos 
es un campo que esta en desarrollo, y que su uso en obras civiles es relativamente nuevo, en 
este artfculo se presenta la teoria de mezclas serie/paralelo [i] como una potente herramienta 
para el analisis numerico. Con el fin de destacar las ventajas que tiene esta teoria en el analisis 
y diseno estructural de estructuras reforzadas con polimeros reforzados con fibras largas 
(FRP), como ejemplo se Simula un conjunto de porticos con diferentes configuraciones de 
refuerzo. 

2. ESTRUCTURAS TIPO PORTICO DE HORMIGON ARMADO REFORZADAS CON FRP 

El refuerzo en las estructuras como tecnica de reparacion mejora las deficiencias 
estructurales que han conducido al deterioro y a la reduccion su funcionalidad. Asimismo el 
refuerzo como tecnica de rehabilitacion es eficiente, en las estructuras que no cumplen con los 
requisitos necesarios para garantizar un buen comportamiento durante eventos extremos bien 
sea por inadecuados disenos o por baja calidad en la construccion, por cambio de las 
condiciones de uso o por la necesidad de adaptar la estructura a los requerimientos de diseno 
actuales [ii]. 

La rehabilitacion y/o la reparacion de las estructuras en hormigon armado se han realizado 
tradicionalmente con laminas de acero. No obstante, debido a la evidencia de la vulnerabilidad 
de las estructuras mostrada durante los sismos de Loma Prieta (1989), de Northridge (1994) y 
de Kobe (1995), se realizaron un gran numero de investigaciones para mejorar las tecnicas 
existentes y estudiar otros materiales como refuerzo con el fin de garantizar un adecuado 
funcionamiento de las estructuras existentes sin causar sobrecostes [iii; iv], 

Un area del refuerzo de estructuras se ha enfocado a investigar los materiales compuestos, 
se han realizado ensayos a diferentes escalas de vigas, columnas, uniones viga-columna, 
muros y losas de entrepiso. Las investigaciones enfatizan que el uso de los materiales 
compuestos en el refuerzo de estructuras son una buena tecnica de refuerzo; sin embargo, 
advierten que se requieren adecuados conocimientos de diseno y de construccion para 
estandarizar metodologfas que garanticen el apropiado uso de este tipo de materiales [v; vi], 

Los materiales compuestos mas utilizados como refuerzo son los polimeros reforzados con 
fibras largas (FRP), predominando el uso de la fibra de vidrio (GFRP) y la fibra de carbono 
(CFRP); en menor proporcion han sido utilizados los materiales compuestos con fibra de 
aramida (AFRP). En los FRP, las fibras largas soportan las acciones mecanicas en una 
direccion predeterminada y la resina o matriz actua como medio para transferir las tensiones 
entre las fibras cercanas garantizando de paso la uniformidad de las deformaciones de las 
mismas [vii]. 

Debido al creciente uso del FRP en la rehabilitacion y reparacion de estructuras de 
hormigon, en algunos pafses se han desarrollado gufas de diseno y construccion del FRP 
como refuerzo para edificios y puentes (ACI440.2R-08 [viii] en Estados Unidos; JSCE-1997 [ix] 
en Japon; Standard S806-02 [x] y CAN/CSA-S6-00 [xi] en Canada; FPI-CEB-2001 [xii] en la 
Union Europea; CNR DT 200/200 4[xiii] en Italia). No obstante, frente a la mayor incertidumbre 
del comportamiento de los sistemas reforzados mediante FRP con respecto a la de los 
sistemas reforzados con laminas acero, las gufas de diseno, ademas de los factores de 
reduccion de resistencia, aplican otros factores de reduccion a la contribucion del FRP 
dependiendo de la solicitacion [xiv], 

A pesar de que los compuestos tienen una resistencia mas alta que la resistencia el acero, y 
que estan conformados por fibras, que como las de carbono, tienen un modulo de elasticidad 
similar o superior al del acero, su uso esta limitado por las deformaciones admisibles de los 
materiales. Por ello, generalmente el FRP en las estructuras de hormigon presenta un 
comportamiento elastico [xv], Esta limitacion ha conducido a que en las diferentes gufas de 
diseno con FRP de refuerzo se deba partir de diferentes hipotesis, por ejemplo, no se 
considera la contribucion del FRP en el diseno a compresion para evitar problemas de pandeo 
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en las fibras mientras que a flexion se limitan las deformaciones del FRP para prevenir su 
delaminacion [xvi]. 

Como en el ejemplo de este artfculo se estudia el comportamiento de los porticos reforzados 
con FRP bajo cargas laterales, a continuacion se hace una breve resena del estado de arte de 
las investigaciones en lo que respecta a columnas y uniones de hormigon reforzados con FRP. 

2.1. Refuerzo en columnas de hormigon armado 

El FRP se aplica como refuerzo a columnas con deficiencias en la ductilidad a flexion y/o 
con insuficiente capacidad a cortante. En estos casos el FRP induce tensiones de 
confinamiento que restringen la dilatacion del hormigon para aumentar la capacidad de carga y 
de deformacion del elemento, ademas de generar un mecanismo que incremente la resistencia 
a cortante [xvii], El incremento del uso del encamisado de las columnas con FRP en lugar del 
encamisado con acero se debe a que su instalacion es mas rapida, es liviano, es adaptable a 
cualquier geometrfa que tenga la columna; ademas su alta resistencia a la corrosion reduce el 
mantenimiento e incrementa la durabilidad del hormigon. 

El comportamiento de las columnas reforzadas con FRP, depende de multiples factores 
entre ellos pueden enumerarse: la configuracion de la columna (geometrfa y armado), la 
configuracion del encamisado (espesor, angulo de la fibra, numero de capas y tipo de 
encamisado), las propiedades mecanicas de los materiales (acero, hormigon y FRP) y el 
estado de dano del hormigon y del acero. 

Es importante considerar los siguientes aspectos del confinamiento lateral dado por el FRP 
[xviii] : (i) El encamisado por su esbeltez, no aumenta la rigidez lateral de la columna; al confinar 
el hormigon solo mitiga los fallos locales prematuros. (ii) Los FRP son susceptibles a la rotura 
en los puntos localizados donde hay demanda de deformacion; por ello, el encamisado con 
FRP no evita el pandeo de la armadura longitudinal, pero disminuye la probabilidad de que 
suceda. (iii) El encamisado reduce la fisuracion del hormigon por code en la zona donde se 
forman las rotulas plasticas, conduciendo a que las fisuras ocurran por flexion cerca a la base; 
no obstante, al incrementarse la demanda, la columna puede colapsar si su armadura no 
cuenta con la suficiente longitud de anclaje. 

Los ensayos sobre columnas reforzadas con diferentes tipos de FRP han mostrado que la 
forma de la seccion transversal de las columnas incide en el grado de confinamiento producido 
por el refuerzo, y que la forma optima es la seccion circular, mientras que con una seccion 
rectangular, el confinamiento tiende a ser bajo [xix; xx], Recientemente se ha estudiado el 
comportamiento de columnas huecas reparadas [xxi] y rehabilitadas [xxii] con FRP, en las 
cuales tambien hay un incremento de ductilidad y resistencia, sin embargo, aun no se ha 
evaluado la eficiencia del confinamiento para este tipo de secciones. 

Xiao y Wu [xxiii] establecieron experimentalmente que la resistencia y el modulo de 
elasticidad del hormigon confinado junto con el espesor del encamisado de FRP, son 
parametros importantes en el comportamiento de los elementos sometidos a cargas axiales, 
donde el fallo del hormigon confinado depende de la rotura del FRP. Parvin y Wang [xxiv] 
investigan el comportamiento en columnas rectangulares reforzadas con FRP bajo cargas 
excentricas, sus resultados experimentales y numericos muestran que se puede incrementar la 
resistencia y la ductilidad de las columnas, no obstante advierten que la eficiencia del refuerzo 
depende de la rigidez del FRP. Ozcan et al. [xxv] de acuerdo con los ensayos realizados 
concluyeron que el numero de capas de FRP utilizado para confinar la base de las columnas, 
mejora apreciablemente la ductilidad del sistema, mientras que el nivel de confinamiento no 
aumenta en la misma proporcion, ademas que el incremento de la resistencia fue pequeno 
(1 0%-1 5%). 

Algunas de las investigaciones experimentales se han complementado con modelaciones 
numericas. Por ejemplo: Parvin y Wang [xxvi] simularon columnas circulares reforzadas con 
FRP con carga axial y bajo dos condiciones de carga lateral, una monotonica y otra cfclica, 
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ademas aplicaron un analisis considerando la no linealidad material y geometrica para 
involucrar degradacion de la rigidez de la estructura. Sus resultados indican que con el refuerzo 
en la base de las columnas se mejora la rigidez y la ductilidad, ademas que el encamisado 
retrasa la degradacion de la rigidez de las columnas. Parvin y Jamwal [xxvii] por medio de un 
analisis no lineal similar al de Parvin y Wang [xxvi], modelaron columnas rectangulares 
encamisadas con FRP bajo fuerzas ciclicas laterales y carga axial, donde consideraron tres 
configuraciones en la variaron la direccion de las capas del refuerzo. Elios senalan que con el 
aumento de espesor del FRP se incrementa la resistencia a la carga axial y la ductilidad de las 
columnas, y que el porcentaje de este incremento depende de la direccion de las capas del 
laminado. 

Uno de los limitantes en la modelacion numerica para el analisis de las columnas reforzadas 
esta en simular el confinamiento dado por el acero y por el FRP, por ello en algunas 
investigaciones se han propuesto modelos para introducir en los elementos finitos, como es el 
caso de Mirmiran et al. [xxviii], Elios proponen un modelo no lineal con plasticidad no asociada 
de Drucker-Prager para tener en cuenta el confinamiento hormigon, mientras que el 
encamisado con FRP lo consideran como un material elastico y lineal. Sin embargo, advierten 
que con este modelo no se puede aplicarse para establecer la perdida de rigidez y de 
resistencia. Pese a que con este tipo de estudios se han realizado importantes avances, en la 
modelacion numerica aun se deben desarrollar diferentes herramientas para simular de una 
forma mas precisa el comportamiento de las columnas reforzadas con FRP. 

2.2 Refuerzo en uniones viga-columna de hormigon armado 

El comportamiento de las uniones viga-columna en una estructura tipo portico es un factor 
importante que influye en la capacidad de resistir cargas sfsmicas y otras cargas laterales. En 
los diferentes sismos se ha identificado que el colapso de estructuras de hormigon armado y de 
acero se ha debido a los fallos en las uniones, principalmente en estructuras construidas antes 
de que se incluyeran los requisitos de diseno sfsmico en las normas de diseno. Estas 
estructuras tienen una baja capacidad de resistencia lateral ya que carecen de detalles que 
desarrollen ductilidad; por ejemplo, tienen uniones sin o con escasa armadura transversal, 
insuficiente longitud de anclaje de la armadura de las vigas y sistemas de viga fuerte-columna 
debil, donde se generan en las uniones o cerca de ellas, inapropiados mecanismos de 
articulacion plastica [Hi]. 

Varias tecnicas de refuerzo han sido aplicadas a las uniones viga-columna deterioradas o 
con deficiencias de diseno utilizando laminas de acero y FRP; sin embargo, la dificultad del 
refuerzo esta en proveer un confinamiento efectivo teniendo en cuenta la complejidad del 
comportamiento de las uniones [ii]. 

La seleccion de una alternativa de refuerzo se debe hacer acorde con el aumento de 
capacidad que se requiera en la union, donde es imprescindible revisar la resistencia de las 
vigas y columnas adyacentes frente a este incremento; asimismo, es necesario tomar las 
medidas de refuerzo necesarias para evitar, tanto el fallo fragil, como la formacion de 
articulaciones plasticas prematuras. Todo ello, permite garantizar la efectividad del refuerzo y 
mejorar comportamiento global de la estructura. 

Dependiendo del comportamiento esperado de las estructuras existentes y del objetivo del 
refuerzo en las uniones, existen diferentes configuraciones de refuerzo con FRP en las uniones 
viga columna. Entre ellas pueden citarse: el refuerzo a cortante en la union [xxix; xxx]; el 
refuerzo en las esquinas de la union colocando el laminado en la cara superior y/o inferior de 
las vigas y en el lado adyacente de la columna [xxxi; xxxii] ; el encamisado con FRP en la union 
y en la columna [xxix; xxx]; asimismo, cuando es necesario, en la zona adyacente a las 
uniones, las vigas se refuerzan a flexion [xxxii; xxxiii] y a cortante [xxxiv; xxxv; xxxvi], incluso 
en algunos casos, las encamisan [xxxvii; xxxviii; xxxix] para mejorar su comportamiento ante 
cargas ciclicas. Segun la configuracion del refuerzo es posible incrementar la capacidad a 
cortante y/o a flexion de la union, aumentar su resistencia al giro con el fin de reducir la 
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probabilidad de deslizamiento de la armadura de las vigas y confinar el hormigon tanto de la 
union como de la columna. 

Las uniones reforzadas con FRP han sido estudiadas en menor proporcion que las vigas y 
las columnas reforzadas con FRP. La mayorfa de los estudios se han enfocado en las uniones 
externas, tal vez por ser este el tipo de union en el que han ocurrido mas fallos. Ademas la 
rehabilitacion y la reparacion en las uniones se han investigado mas a nivel experimental, pero 
en el campo de la simulacion numerica solo ha sido considerado la rehabilitacion. 

El comportamiento a cortante de las uniones reforzadas con FRP se ha estudiado con 
mayor interes, debido a que los fallos por cortante de las uniones viga-columna durante los 
sismos han sido la principal causa de colapso de estructuras porticadas. Entre las 
investigaciones que han estudiado el efecto del cortante en las uniones exteriores se 
encuentran las siguientes: Gergely et al. [xxxvii], senalan que el FRP mejora su capacidad a 
cortante, pero que la delaminacion reduce la eficiencia del refuerzo, por ello enfatizan en la 
importancia de la adherencia entre el laminado y el hormigon. Ghobarah y Said [xxix] 
investigaron uniones tipo viga fuerte-columna debil con tres configuraciones de refuerzo; ellos 
concluyeron que para garantizar la efectividad del refuerzo, ademas de requerirse una optima 
adherencia entre el hormigon y el FRP, es necesario que se proporcione confinamiento en la 
union. Said y Nehdi [xxx] resaltan que al reparar o rehabilitar las uniones con FRP se aumenta 
la resistencia y la ductilidad global de la estructura, pero que en la reparacion de las uniones es 
necesario considerar las limitaciones debidas al dano preexistente en el hormigon y en el 
acero. Antonopoulos et al. [xxxiv] ensayaron uniones tipo viga fuerte-columna debil con 6 
configuraciones de refuerzo, indicando que el deslizamiento del FRP afecta el comportamiento 
de la union, al menos que se mejore la adherencia con anclajes mecanicos. Aunque en menor 
proporcion, tambien se ha estudiado el comportamiento a cortante de las uniones internas 
reforzadas con FRP [xxxv; xxxix], donde el refuerzo aumenta la capacidad y la ductilidad de la 
union, asimismo, la eficiencia depende de la adherencia entre el FRP y el hormigon. 

Existen un numero reducido de investigaciones que han estudiado otros parametros: 
Granata y Parvin [xxxi; xl] ensayaron y modelaron uniones exteriores reforzadas donde 
mostraron que el FRP mejora la capacidad a flexion y disminuye la rotacion en la union; 
tambien identificaron que el espesor del refuerzo incide en la concentracion de tensiones que 
ocurren en las fibras alrededor de las esquinas de la union y que conllevan a la delaminacion. 
Ghorabah y El-Amoury [xxxii] estudiaron cuatro configuraciones de refuerzo de uniones 
externas con insuficiente longitud de anclaje en la armadura de las vigas, y senalan que 
ademas de evitar el fallo por code en la union, es necesario que se conserve la integridad del 
hormigon en la zona de anclaje de la armadura de la viga para que se desarrolle su resistencia 
y no se produzca un fallo prematuro. Pantelides et al. [xxxvi] ensayaron las uniones internas 
reforzadas con insuficiente longitud de anclaje de la armadura de las vigas e indican que el 
FRP reduce la propagacion del dano y la perdida de rigidez, sin embargo advierten que con el 
refuerzo se retrasa el deslizamiento de la armadura pero no lo evita. Parvin y Wu [xxxvii] 
hicieron una modelacion de cuatro uniones externas reforzadas con cuatro capas de CFRP 
colocadas en diferentes direcciones y concluyeron que el incremento de la resistencia lateral en 
las uniones depende de la direccion en que se coloque cada una de las capas y de su 
adherencia, lo que concuerda con lo estudios experimentales de otros investigadores [xxxiv; 
xxxix]. 

3. SIMULACION NUMERICA DEL COMPORTAMIENTO DE LOS MATERIALES 
COMPUESTOS 

El comportamiento de las estructuras en hormigon armado reforzadas con FRP, depende de 
multiples parametros y el grado de influencia de muchos de ellos no ha sido aun determinado, 
lo que dificulta que experimentalmente se puedan estudiar todas las variables. Por otra parte, la 
heterogeneidad y la anisotropfa propias de los compuestos han hecho a que no haya metodos 
de analisis sencillos y efectivos que permitan determinar su influencia sobre el comportamiento 
en las estructuras [vi]. 
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Como consecuencia el desarrollo para optimizar el uso del FRP en las estructuras se debe 
realizar de forma conjunta entre el campo experimental y el campo de la simulacion numerica. 
No obstante, es necesario el desarrollo de nuevas herramientas y mejorar las existentes, en 
busca de que la simulacion numerica a nivel global y local represente de una forma mas 
aproximada el comportamiento real de las estructuras reforzadas con FRP, con la finalidad 
obtener informacion que no se puede medir en los ensayos. 

3.1 Teonas de simulacion del comportamiento de los materiales compuestos 

Los materiales compuestos estan conformados por diferentes tipos de sustancias 
inorganicas u organicas, cada uno de los materiales componentes teniendo su ley constitutiva 
que condiciona el comportamiento del conjunto en funcion de la proporcion volumetrica y de la 
distribucion morfologica que tenga dentro del compuesto [xli], 

En la actualidad, gran parte de las investigaciones acerca del comportamiento de los 
materiales compuestos corresponden al campo experimental. Su estudio por medio de 
modelaciones numericas esta en desarrollo pero presenta algunas restricciones; incluso en 
simulaciones recientes de elementos reforzados con FRP, se representa el comportamiento 
constitutivo del material compuesto como el de un material homogeneo, sin tener en cuenta el 
aporte de sus componentes. 

Con el objetivo de solucionar esta limitacion, se han propuesto diferentes teorfas de 
simulacion del comportamiento de los materiales compuestos como gestores del uso de los 
modelos constitutivos que representan el comportamiento de los materiales simples que 
componen los materiales compuestos. Estas teorfas de simulacion son herramientas que 
pueden ser utilizadas dentro de una tecnica de aproximacion de elementos finitos para simular 
apropiadamente el comportamiento de los materiales compuestos, entre ellas las mas 
relevantes son: 

3.1.1 Teoria de la homogenizacion. Dentro de la mecanica de medios continuos 
requiere del analisis desde dos escalas diferentes: una macroscopica en la que 
se determina el comportamiento de la estructura; y una microscopica en la que 
se obtiene el comportamiento del compuesto partiendo de la respuesta de sus 
componentes [xlii], 

3.1.2 Teoria de las mezclas. Considera que el comportamiento de cada 
componente define el comportamiento global del compuesto. Parte de la 
mecanica del continuo local bajo el principio de interaccion de las sustancias 
que componen el material, suponiendo que en el volumen infinitesimal del 
compuesto participan en conjunto todos sus componentes. Asimismo, 
considera que cada uno de ellos contribuye al comportamiento del compuesto 
en la misma proporcion que su participacion volumetrica [xliii]. 

Este artfculo esta orientado bajo el lineamiento de la Teoria de Mezclas, tomando como 
base la teoria de mezclas serie/paralelo propuesta por Rastellini [i], y que ha sido validada a 
traves de la comparacion de los resultados con diferentes pruebas experimentales [i; xliv], Esta 
teoria se fundamenta en la teoria de mezclas clasica inicialmente estudiada por Trusdell y 
Topin (1960). 

3.2. Anisotropia en los materiales compuestos 

Debido a las diferencias en el comportamiento entre los componentes y a la forma en la que 
estan distribuidas las fibras y la matriz, los materiales compuestos tienen una elevada 
anisotropia, ademas de un comportamiento no lineal [xli]. La elevada anisotropia de los 
materiales compuestos tambien conduce a que su simulacion a traves del metodo de los 
elementos finitos sea compleja, siendo esta una de las limitaciones del uso de la simulacion 
numerica. No obstante, Car et al. [vii], con base en la teoria del mapeo de espacios y la teorfas 
de mezclas de los materiales compuestos logran involucrar cada uno de los componentes del 
compuesto en el analisis de una estructura, donde la anisotropia del compuesto es el resultado 
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de la participacion de cada uno de los materiales que lo componen con su respectivo 
comportamiento (isotropo, ortotropo o anisotropo). 

Oiler et al. [xlv] proponen la teorfa del mapeo de espacios para modelizar el complejo 
mecanismo de fallo de los materiales anisotropos, el cual con base en el concepto de mapeo 
de los tensores de tensiones y de deformaciones desde un espacio anisotropo real a un 
espacio ficticio isotropo, estableciendo una relacion biunfvoca entre el comportamiento 
anisotropo del material real y el comportamiento de un material isotropo ficticio cualquiera. Con 
la generalizacion de esta teorfa para el campo de grandes deformaciones propuesta por Car et 
al. [vii], se puede formular en la configuracion referencial o espacial utilizando la cinematica 
lagrangeana total o actualizada respectivamente. Esta formulacion parte del concepto de 
mapeo de tensiones propuesto por Betten, el cual permite el transporte en forma lineal del 
tensor de tensiones de un espacio real a un espacio ficticio. Cabe anotar que al hacer la 
trasformacion lineal se garantiza la convexidad de la funcion de fluencia y del potencial plastico, 
cumpliendose con ello la segunda ley de la termodinamica. Sin embargo, la propuesta de hacer 
solo el mapeo de tensiones parte de la hipotesis que las deformaciones elasticas son identicas 
en ambos espacios, lo cual es una limitacion por excluir el analisis del comportamiento de los 
materiales no proporcionales como los materiales compuestos [xli], Para evitar esta limitacion, 
la teorfa del mapeo de espacios, aparte del transporte en el espacio de tensiones, propone el 
transporte en el espacio de deformaciones. 

En la figura 1 se observan los espacios de tensiones y de deformaciones anisotropo real e 
isotropo ficticio; ademas se indica la relacion entre ellos a traves de dos tensores simetricos de 
cuarto orden ( A* 1 , A E ), los cuales contienen la informacion sobre las propiedades del material 
anisotropo obtenidas mediante ensayos experimentales [xlvi]. 

La teorfa de mapeo de espacios para materiales anisotropos permite generalizar la 
formulacion clasica isotropa en el analisis del comportamiento de materiales anisotropos. 
Basicamente, el analisis sobre cada punto de Gauss consiste en transportar el tensor predictor 
de tension y el tensor de deformacion de un material anisotropo de un espacio real anisotropo a 
un espacio ficticio isotropo; luego, en el espacio ficticio isotropo se realiza la integracion de la 
ecuacion constitutiva isotropa para obtener el nuevo estado tensional y el tensor constitutivo 
tangente correspondiente; una vez calculados los tensores resultantes, los tensores de tension, 
de deformacion y constitutivo tangente, se transportan al espacio real anisotropo [xlvii]. 




Figura 1 Relacion entre el espacio isotropo ficticio y el espacio anisotropo real para pequenas 
deformaciones [xlvi]. 
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La implementacion del mapeo de espacios de tensiones y de deformaciones es una 
herramienta numerica util para trabajar con materiales anisotropos, ya que permite el uso de 
las funciones de fluencia, potenciales plasticos y algoritmos de integracion de las ecuaciones 
constitutivas para materiales isotropos en el analisis de comportamiento de materiales 
anisotropos [xli], Ademas, al integrarlo con la teoria de mezclas, permite que el estudio de la 
anisotropia de los compuestos no sea una restriccion. 

3.3 Teona de mezclas serie/paralelo 

En la teoria serie/paralelo, se considera que en una direccion particular los componentes se 
comportan en paralelo y en las otras direcciones en serie [i]. En el caso de los polimeros 
reforzados externamente con fibras de carbono (CFRP), la direccion en paralelo corresponde a 
la direccion de la fibra y, en el caso del hormigon armado, a la direccion de la armadura. La 
teoria parte de las siguientes hipotesis: 

• En cada volumen infinitesimal del compuesto participan en conjunto todos sus 
componentes, es decir que la distribucion de los componentes es homogenea. 

• Los componentes tienen una misma deformacion en la direccion en paralelo (condicion 
de isodeformacion). Los componentes tienen una misma tension en la direccion en 
serie (condicion de isotension). 

• La adherencia entre los componentes es perfecta. 

• La contribucion de los componentes en la respuesta del compuesto es proporcional a 
su participacion volumetrica de cada material componente. 

• El volumen ocupado de cada material componente es inferior al volumen total del 
compuesto. 

Esta formulacion plantea combinar el comportamiento de los materiales componentes con 
el fin de obtener la respuesta del material compuesto. Por ello, la anisotropia global de los 
materiales compuestos se considera como el resultado de la interaccion de los componentes. 
Adicionalmente, la teoria serie/paralelo permite que se analicen los materiales componentes 
que presenten cualquier tipo de no linealidad como el dafio y la plasticidad. 

3.3.1 Definicion de los componentes en serie y en paralelo de los tensores 
deformacion y de tension [xlviii] 

El tensor deformacion s se descompone una parte en serie e s y otra en paralelo e p por 
medio de los tensores de proyeccion (P P ,P S ) de cuarto orden en paralelo y en serie 
respectivamente: 

£ = £p £ s ( 1 ) 

s P = Pp : s = P s : e (2) 

El tensor de tensiones a tambien se separa en sus componentes en serie a s y en paralelo cr P : 

CT = <3 p + <J S (3) 

= P p : ct a s = P s : a (4) 

donde los tensores de proyeccion de cuarto orden se hallan a traves del tensor de proyeccion 
paralelo de segundo orden N P y del tensor identidad de cuarto orden I , 

P F = /Vp®/Vp P S =I-P p (5) 

N P es funcion del vector unidad e ? que indica la direccion del comportamiento en 
paralelo, es decir, la direccion de la fibra 

N p = e l ®e 1 (6) 

3.3.2 Ecuaciones de equilibrio y de compatibilidad en las capas del compuesto 

[xlix] 

Con el fin de minimizar el costo computacional, la implementacion numerica de la teoria fue 
desarrollada para descomponer el compuesto c en un determinado numero de capas neap, tal 
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que cada capa / este conformada por una matriz m y una fibra f. Consecuentemente, de 
acuerdo con las hipotesis enunciadas, en cada capa del compuesto / se plantean las siguientes 
ecuaciones de equilibrio y de compatibilidad de deformaciones: 



• Comportamiento en paralelo 



jE p - j£ p 



P 

j b P 



• Comportamiento en serie 



c mi m , f 7 f 

j<*P=J k ’ j°P+J k -J°r 



( 7 ) 

( 8 ) 



c m f 

j°s = J°s - jVs 



0) 



) e s = "j k • "j e s + f j k ' f j e s ( 1 °) 

Para garantizar el principio de conservacion de masa en cada capa, los porcentajes de 
participacion volumetrica de la matriz y de la fibra de la capa / ( "‘k , jk ) deben cumplir con la 
siguiente ecuacion 



d "V 

"Jc - 1 

1 dV; 






lk+jk = l 



( 11 ] 



donde"y representa el volumen de la matriz m o de la fibra f de la capa / , y V ; es el volumen 
de la capa /del compuesto analizado c. 



3.3.3 Ecuacion constitutiva de los materiales de las capas del compuesto 

Como la implementacion de la teorfa serie/paralelo se realiza a nivel constitutivo a partir del 
estado de deformacion en un punto de Gauss, una vez obtenida la deformacion del compuesto 
c e, se calcula el estado tenso-deformacional de cada componente cumpliendo con las 
ecuaciones de equilibrio y de compatibilidad, para luego hallar el estado de tensiones y de 
deformaciones del compuesto. El estado de tensiones y de deformaciones de los componentes 
se determina a traves de la ecuacion constitutiva que rige el comportamiento de cada uno: 



"<J= "C: "£ e = "C: "(E-E P ) n = {mj) (12) 

donde "C es el tensor constitutivo de cada componente de la capa /del compuesto. 

La descomposicion del operador tangente inducida por la separacion del comportamiento en 
serie y paralelo de cada componente [i], queda establecida como 

~".C 



n C = 



j -'pp 

n 

j^SP 



j^PS 

n 

j^SS 



{»!,/} 



n 

j^pp - 



n 

j^SP ~ 



UjUp 


— P 


. n , 


■ P 


n 


UjUp 


d'fr 


— * / • 


' j 




j^PS 






— P 


:”C: 


• P 


n 

^SS 


Fi°s 


d’]B p 


— *v • 







= P„ 



= P 



"C • P 

j ^ " r s 



"C • P 

jK, . r s 



(13) 



3.3.4 Algoritmo de la teorla serie/paralelo 

En la figura 2 se observa el algoritmo de desarrollo de la teorfa serie/paralelo. Se observa 
que, una vez que se determina el tensor deformacion del compuesto, en cada capa / se 
descompone el tensor de deformacion en sus partes en serie y en paralelo, con el fin de 
calcular las correspondientes deformaciones en la matriz y en la fibra. 
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Deformaciones de la capa de compuesto j en un PG 
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Prediccion deformaciones en serie de la matriz : ,■ S s 



Deformaciones en 
serie de la fibra: 
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Determination elastica 
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Figura 2 Esquema de solucion de la teoria de mezclas serie/paralelo en pequenas 
deformaciones para una capa de un compuesto en un punto de gauss para un paso i+1 [li] 



Segun la teoria serie/paralelo, la deformacion en paralelo de los componentes es la 
misma, mientras que la deformacion en serie es diferente para cada componente; por ello, en 
el analisis de la parte en serie requiere de una primera aproximacion de la deformacion en serie 
de alguno de los componentes. En este caso, en la implementacion se tomo como predictor, el 
tensor de deformacion en serie de la matriz del compuesto para determinar el tensor de 
deformacion en serie de la fibra con base en la ecuacion (1 0) [xliv], 

El predictor del tensor de deformacion para el paso de carga i+1 en la iteracion k, parte 
de los valores convergidos del paso anterior / y la prediccion del incremento del tensor de 
deformacion en la iteracion k 

< 14 > 



Para la primera prediccion se propone la ecuacion (15), donde se considera que las 
componentes en serie y paralelo son distribuidas de acuerdo con el tensor constitutive de los 
componentes determinado en el paso anterior convergido / [i]. 
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,A=(^'C SS + 7 *"C») ’ 



m ,c SP ): 



/+1 Vc 1 ‘ 


r i 




L j b *1 



(15) 



La ecuacion (9) de equilibrio de tensiones, se verifica calculando la tension residual de 
la capa jA<j s , comprobando que sea menor que una tolerancia Toler 



A°s = "As - f j°s £ Toler 



(16) 



Cuando la tension residual es mayor, la prediccion de la deformacion de la matriz se corrige 
en forma iterativa. Una vez obtenida una tension residual menor que la tolerancia, se hace la 
recomposicion de los tensores de tension y de deformacion de la capa j. 



3.3.5 Ecuacion de equilibrio en el compuesto 



El analisis de cada capa proporciona el tensor de tension del compuesto C cr como a la suma 
de los tensores de tension de las capas del compuesto neap ponderados por el porcentaje de 
participacion volumetrica de cada capa f k 



(17) 



C(7 = 1L J k > = Ya j k ^ C Pr+ C Ps) 

7=1 7=1 

Para garantizar el principio de conservacion de masa en el compuesto, los porcentajes 
de participacion volumetrica de las capas del compuesto deben cumplir la condicion 

r ]\z neap 

i k = => k = l (18) 



dV„ 



7=1 



donde V, es el volumen de la capa / del material y V r corresponde el volumen total del material 

J ^ 

compuesto c. 



4. APLICACION DE LA TEORIA SERIE/PARALELO A LA SIMULACION NUMERICA EN 
2 DIMENSIONES 



El estudio de los efectos que han producido los grandes sismos en estructuras tipo portico 
pone en evidencia que las zonas mas susceptibles a dano son las uniones viga columna y las 
bases de las columnas [iii]. Por ello, con el proposito de garantizar la estabilidad de las 
estructuras durante un evento extremo, en muchos estudios se hace entasis en la necesidad 
de la rehabilitacion y reparacion de las estructuras antiguas o de las estructuras construidas 
antes de las actuales normas de diseno, siendo una de las alternativas de refuerzo el uso los 
polimeros reforzados con fibras largas (FRP) [v; vi; xiv; xv; xvi; xvii]. 



Utilizando la teorfa de de mezclas serie/paralelo en el programa de elementos finitos 
PLCDYN [I], se muestra en este apartado un estudio de estructuras porticadas reforzadas con 
FRP. Se realizo un analisis no lineal estatico incremental (pushover analysis) de diez 
estructuras planas con una misma geometria. Cinco de ellas son de hormigon simple y las 
otras son de hormigon armado, con distintas configuraciones de refuerzo utilizando polimeros 
reforzados con fibras de carbono CFRP [li]. 



Aunque los FRP como refuerzo se emplea en estructuras de hormigon armado o de acero, 
se estudio tambien, el comportamiento de las estructuras de hormigon simple reforzadas, con 
el proposito de analizar unicamente la influencia del refuerzo de FRP sobre las estructuras. 
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4.1. Descripcion de los modelos estudiados 
4.1.1 Geometria y configuraciones de refuerzo 



5.00 



Seccion transversal 
de la viga 



Seccion transversal 
del pilar 



REFUERZO DELOS 
PORTICOS EN HORMIGON 
ARMADO 



104)16 



l$10c/0.20m 



Refuerzo de la viga 



16 ( 1)16 

\ -'/\y 



l())8c/0.20m 



Refuerzo del pilar 



Figura 3 Geometria tipo de las estructuras porticadas [xliv] 



En la figura 3 se indica la geometria de los diez modelos junto con las armaduras de la 
viga y de las columnas para los porticos en hormigon armado. En la figura 4 se presentan las 
configuraciones del refuerzo con CFRP, y en la tabla 1 se indica la nomenclatura utilizada para 
identificar los modelos. 





Figura 4 Modelos de los porticos con las diferentes configuraciones de refuerzo [li] 
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Tabla 1 Nomenclatura de los modelos de los portico. 



Tipo de hormigbn 

Tipo de portico _ 


Sin 

armadura 
de acero 


Con 

armadura 
de acero 


Portico sin CFRP (Figura 4 (a)) 


SASF 


ASF 


Portico con CFRP en la union (Figura 4(b)) 


SAF 


AF 


Portico con CFRP en la union y en la base de la columna (Figura 4 (c)) 


SAFC 


AFC 


Portico con CFRP en la union y extension en la viga (Figura 4 (d)) 


SAFV 


AFV 


Portico con CFRP en la union, extension en la viga y en la base de la 
columna (Figura 4(e)) 


SAFCV 


AFCV 



El refuerzo en las columnas de estos modelos corresponde a dos capas orientadas a 0 
y a 90 grados para tener en cuenta que la eficiencia del encamisado de la columna depende de 
las direcciones en que se coloque la fibra [xxvii], Asimismo, cuatro de los diez porticos 
analizados tienen CFRP en la base de las columnas dado que segun los resultados 
experimentales [xxv; xxvi] su nivel de influencia es notorio en la capacidad y en la ductilidad en 
las columnas. Con el fin de tener un mejor comportamiento de la viga ante un desplazamiento 
lateral, en los modelos con refuerzo se encamiso la viga en las zonas cercanas a las uniones 
viga-columna de acuerdo con los estudios experimentales [xxxvii; xxxviii]. Se seleccionaron dos 
longitudes de refuerzo en la viga para hacer una comparacion del comportamiento cuando el 
refuerzo se extiende desde el borde de la union dos y cuatro veces la altura efectiva de la viga. 

DETALLE DE LA VIGA DETALLE DEL PILAR DETALLE DE LA UNION 

EN HORMIGON EN HORMIGON EN HORMIGON 












Seccion Seccion 
sin CRFP con CRFP 



DETALLE DE LA VIGA 
EN HORMIGON ARMADO 





P 







Seccion Seccion 
sin CRFP con CRFP 



MATERIALES COMPUESTOS 
EN LA VIGA 




Material 1 
Material 5 
Material 6 
Material 7 
Material 12 
Material 13 
Material 15 



T/ 



id. 




Seccion 
sin CRFP 

Seccion 
con CRFP 



DETALLE DEL PILAR 
EN HORMIGON ARMADO 




Seccion 
sin CRFP 



Seccion 
con CRFP 



MATERIALES COMPUESTOS 
EN EL PILAR 

Material 1 
Material 2 
Material 3 
Material 7 
Material 8 
Material 9 
Material 14 





DETALLE DE LA UNION 
EN HORMIGON ARMADO 




MATERIALES COMPUESTOS 
EN LA UNION VIGA -PILAR 





Material 4 
Material 7 
Material 9 
Material 10 
Material 11 
Material 14 
Material 15 



Figura 5 Configuracion de los materiales compuestos en los porticos [li] 
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Tabla 2 Porcentajes de los componentes en los materiales compuestos de los porticos [li] 



Material 

compuesto 


Matriz de 
Hormigon 


Acero 

longitudinal 


Acero 

Vertical 


Cercos 

Horizontales 


Cercos 

Verticales 


Matriz 

Polimerica 


CFRP 0 s 
Horizontal 


CFRP 

90 s 

vertical 


1 


100 
















2 


98.70 




1.10 


0.20 










3 


97.53 




1.10 


0.17 




0.40 


0.40 


0.40 


4 


97.71 




1.09 






0.40 


0.40 


0.40 


5 


98.70 








0.10 


0.40 


0.40 


0.40 


6 


99.90 








0.10 








7 


98.80 










0.40 


0.40 


0.40 


8 


87.13 




12.67 


0.20 










9 


86.04 




12.59 


0.17 




0.40 


0.40 


0.40 


10 


81.51 


4.84 


12.45 






0.40 


0.40 


0.40 


11 


92.87 


4.84 


1.09 






0.40 


0.40 


0.40 


12 


93.86 


4.87 






0.07 


0.40 


0.40 


0.40 


13 


95.00 


4.90 






0.10 








14 












34.00 


66.00 




15 












34.00 




66.00 



En estos ejemplos, con la teoria serie/paralelo las matrices de los compuestos son el 
hormigon y la resina polimerica, mientras que el acero y la fibra de carbono constituyen el 
refuerzo. En la figura 5 se indican los materiales compuestos asignados a los elementos 
estructurales de los modelos, dependiendo de si tienen o no armadura. Los porcentajes 
volumetricos en los que participan los componentes en cada material compuesto se senalan en 
la tabla 2, donde se especifica la direccion de las fibras del compuesto considerada como parte 
de la anisotropia del compuesto. Las propiedades de los materiales simples se muestran en la 
tabla 3. 



Tabla 3 Propiedades mecanicas de los materiales componentes de los compuestos en los 
porticos 



Material 


Criterio de 

fluencia 


Ex=Ey=Ez 

[MPa] 


V 


(JC 

[MPa] 


at 

[MPa] 


Gc 

[kPa-m] 


Gt 

[kPa-m] 


Hormigon 


Mohr-Coulomb 


2.5-10 4 


0.20 


30.0 


3.0 


50.0 


5.0 


Acero 


Von Mises 


2.1-1 0^ 1 


0.00 


270.0 


270.0 


2000.0 


2000.0 


Matriz 

polimerica 


Mohr-Coulomb 


1.2-10 4 


0.20 


87.5 


29.2 


36.0 


3.0 


Fibra de 

carbono 


Von Mises 


1.5-10 5 


0.00 


2300.0 


2300.0 


2000.0 


2000.0 



4.1.2 Descripcion de los modelos constitutivos aplicados a los materiales 

Los modelos analizados por control de desplazamientos, tienen una malla de elementos 
finitos rectangulares de 4 nodos. En lo que respecta al analisis del comportamiento de los 
materiales simples, para determinar el dano en el hormigon y la matriz polimerica del 
compuesto se aplico el modelo de dano de Kachanov, el acero se considero como un material 
elasto-plastico utilizando el algoritmo Euler, mientras que la fibra de carbono se analizo como 
un material elastico y lineal. En todos los casos se ha partido de la hipotesis que no hay dano 
inicial en el portico, por lo que los ejemplos corresponden al caso rehabilitacion con FRP. 
Ademas se supuso la adherencia inicial entre el refuerzo y el hormigon es perfecta. 
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4.2 Analisis del comportamiento de los modelos de hormigon simple y de hormigon 
armado 

4.2.1 Resistencia y rigidez 

En la figura 6 se observa que los diferentes modelos tienen una rigidez inicial similar, 
su comportamiento se conserva lineal aproximadamente hasta los 0.0015m; luego, segun la 
configuracion del refuerzo del sistema, la perdida rigidez varia de acuerdo con el desarrollo del 
dano y de plasticidad. La resistencia maxima de cada caso se indica en la tabla 4. 

En lo que respecta al comportamiento de las estructuras de hormigon simple, la 
resistencia del sistema aumenta casi un 40% cuando se coloca CFRP en la base de las 
columnas Por el contrario, en los casos con CFRP en la union y en la viga (SAF y SAFV), el 
aumento de la resistencia del sistema es pequeno, un 5% y un 6%, respectivamente. 
Asimismo, se observa que los modelos con refuerzo adicional en la viga (SAFV y SAFVC) 
muestran un comportamiento ductil despues de alcanzar su capacidad maxima y, en los dos 
casos conservan una perdida de resistencia menor al 24%. Se aprecia como el modelo SAFC 
tiene una tendencia de perdida de resistencia similar al modelo SAF, donde el refuerzo en las 
columnas incrementa la resistencia del sistema pero no incide en su ductilidad. 

De otra parte, en los porticos de hormigon armado, al reforzarse solo la union (AF) hay 
un pequeno incremento en la resistencia del orden de 7%. Cuando el refuerzo se coloca en la 
base de las columnas (AFC) o mas alia de la zona de confinamiento en las vigas (AFV), el 
aumento de la resistencia Mega a un 27%. Al reforzar tanto la base de las columnas como la 
viga (AFCV), la resistencia se incrementa un 50%. A diferencia de los modelos con hormigon 
simple, en todos los casos con armadura hay ductilidad independientemente del refuerzo. En 
los porticos de hormigon armado, el FRP como refuerzo, aumenta la capacidad de la 
estructura, dado que al disminuir el nivel de dano en el hormigon, posterga la plastificacion del 
acero. 




Figura 6 P-A porticos de hormigon simple y de hormigon con armadura de acero [li] 
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Tabla 4 Fuerzas maximas y desplazamientos en la mitad del claro de la viga. 



Hormigon sin 

armadura 


SASF 


SAF 


SAFC 


SAFV* 


SAFCV* 


Fuerza maxima (kN) 


62.36 


65.48 


87.65 


66.06 


89.73 


Apmax (m) 


0.0030 


0.0030 


0.0048 


0.0031 


0.0048 




Hormigon con 

armadura 


ASF 


AF 


AFC 


AFV 


AFCV 


Fuerza maxima (kN) 


183.62 


196.54 


233.98 


234.49 


274.96 


Apmax (m) 


0.0354 


0.0318 


0.0312 


0.0233 


0.0222 



* Modelo con perdida convergencia 



Es relevante notar que el refuerzo no aumenta la rigidez lateral de los porticos, lo que 
coincide con la observacion de Tastani y Pantazopoulou [xviii]. Con el refuerzo hay un 
incremento de la resistencia a cortante y a flexion [xxiv; xxv; xxix; xxx; xxxi], no obstante, el 
nivel de incremento de resistencia depende de la configuracion del refuerzo que se seleccione 
y de la armadura de los porticos. En lo que respecta a la ductilidad, cabe anotar que su 
aumento depende de la distribucion de la armadura; si la armadura es insuficiente, el refuerzo 
con FRP incrementa considerablemente la ductilidad del sistema; pero si la armadura aporta 
por si misma ductilidad al sistema, al reforzarlo el incremento en la ductilidad es imperceptible. 

4.2.2 Tensiones y deformaciones en los componentes de los materiales 
compuestos 

Una de las ventajas de la teorfa serie/paralelo es que permite identificar durante el 
proceso de carga, el estado de tensiones y de deformaciones tanto del compuesto como de 
sus componentes. Para observar la evolucion de las tensiones en los diferentes materiales, se 
tomaron dos puntos de referencia que se indican en la figura 7: uno de ellos, en la esquina 
izquierda, interna donde se analiza las deformaciones y las tensiones en x; y el otro, en la base 
de la columna izquierda, para el analisis de las deformaciones y las tensiones en y. 




En la figura 8(a) se muestran las tensiones en xdel hormigon. Se aprecia que en los 
modelos SASF y ASF se tiene una perdida de resistencia inicial similar, con la diferencia que 
en el caso SASF, debido a la mayor degradacion del hormigon, hay una descarga total, 
mientras que en el modelo ASF, el elemento pierde resistencia sin presentar descarga. Los 
modelos de hormigon sin armadura (SAF, SAFV, SAFC y SAFCV) tienen una perdida de 
resistencia inicial similar; los modelos SAF y SAFC presentan descarga total, posiblemente, por 
la perdida de resistencia del sistema; entretanto los casos SAFV y SAFCV no evidencian 
descarga, aunque cabe anotar que estos modelos perdieron convergencia antes de aplicar el 
40% del desplazamiento total. En lo que respecta a los porticos de hormigon armado con 
refuerzo, los modelos tienen una curva de perdida de resistencia similar. La diferencia consiste 
en que los casos AF y AFC muestran una descarga parcial equivalente a un 28% de la 
deformacion maxima cuando se refuerza la viga (AFV y AFCV). 
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Referente a las tensiones en yen el hormigon indicadas en la figura 8 (b), los modelos 
de hormigon simple que presentan perdida de resistencia son los reforzados en la base de la 
columna (SAFC y SAFCV); los otros tienen un comportamiento lineal y en todos los casos hay 
una descarga total. Los modelos de hormigon con armadura pero sin refuerzo en la base de la 
columna (AF y AFV) muestran una perdida de resistencia similar alcanzando un porcentaje de 
perdida del 75%; asimismo, el comportamiento de los casos AFC y AFCV es similar, con un 
15% perdida de resistencia. 




-t- Sff — 5«FC — 5»FV -^S*FCV \=F -a- V AFC -p-^FV -^-^FCV| 

(a) Tension -deformacion en x 

IS — 




o in: dim soce oioe am a. mi 



| -t- 5«F ~^S«FC SlfV -*-S*FCV -a- V -g-^FV -^-^PCV| 

(b) Tension deformacion en y 

Figura 8 Tension deformacion a traccion en el hormigon. [li] 
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(a) Tension deformacion en x 
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(b) Tension deformacion en y 

Figura 9 Tension deformacion a traccion en la matriz del CFRP. [li] 



En la figura 9 se observa el comportamiento de la matriz del CFRP. En la direccion xen 
la esquina de la union, pese a que la fibra de carbono se mantiene en el rango lineal elastico y 
que solo se ha alcanzado un 30% de la resistencia a traccion de la matriz (Figura 9 (a)), su 
comportamiento se ve afectado por el dafio y perdida de resistencia del hormigon. El 
comportamiento no lineal de la matriz se debe a que con el dano del hormigon, en esta zona 
hay grandes desplazamientos a nivel local. Dependiendo de lo que sucede en el hormigon, la 
matriz en los porticos con hormigon simple, los casos SAF y SAFC se descargan totalmente, 
mientras que los modelos SAFV y SAFV solo presentan perdida de resistencia. La perdida de 
resistencia inicial de los casos SAFC y SAFCV es similar a la de los modelos de hormigon 
armado hasta alcanzar un 50% de la deformacion maxima de los modelos AFV y AFCV. En los 
porticos con hormigon armado, los modelos AF y AFC se descargaron parcialmente con un 
80% de la deformacion maxima de los casos AFV y AFCV. 
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Por contrario en la direccion y de la base de la columna, como se ve en la figura 9 (b), 
aunque hay dano y perdida de resistencia en el hormigon, el comportamiento de la matriz 
polimerica es lineal y elastico para todos los casos. Esto es debido a que en esta zona no hay 
desplazamientos locales que afecten el comportamiento lineal de la matriz. Los porticos de 
hormigon simple, por su baja capacidad carga lateral, tienen una tension maxima alrededor del 
10% de la resistencia a traccion, mientras que los porticos de hormigon armado alcanzan un 
80%. 



4.2.3 Dario 

Otra de las ventajas que tiene el uso de la teorfa serie/paralelo, es que al asignar un 
modelo constitutive a cada componente de los compuestos, se puede observar la evolucion de 
las variables internas tales como el dano o el endurecimiento plastico en el compuesto y en los 
componentes. En la figura 10 se muestra la distribucion de dano en los porticos para el estado 
de carga ultima. 

En las estructuras de hormigon simple, el dano es localizado y la perdida de capacidad 
es considerable, mientras que en las estructuras de hormigon armado el dano se distribuye a lo 
largo de los elementos y la perdida de capacidad del sistema es pequena. 

Dependiente del dano que se produce en los diferentes casos, se aprecia que en las 
columnas y en las vigas es necesario que la longitud del refuerzo externo sea suficiente, para 
evitar el dano localizado en el hormigon y retrasar, de esta manera, la plastificacion del acero, 
previniendose el fallo prematuro del sistema. Como un ejemplo, se observa en las estructuras 
de hormigon reforzado, al comparar los modelos AFV con AF y AFCV con AFC, que los 
porticos con mayor longitud de refuerzo en las vigas tienen menor perdida de rigidez, aunque 
presentan mayor dano en la union y tienen una resistencia lateral un 12% superior. Asimismo, 
al contrastar los casos AF con AFC o SAF con SAFC, se observa que al colocar el refuerzo en 
la base de las columnas, la zona con dano en toda la seccion transversal se traslada del apoyo 
de la columna al borde en el que termina el FRP, con lo cual, el refuerzo aumenta la resistencia 
del sistema. 

Se corrobora que al reforzar las estructuras con CFRP se aumenta la ductilidad y se 
reduce la perdida de resistencia del sistema estructural. Sin embargo, el nivel de eficiencia del 
refuerzo depende tanto de la configuracion de refuerzo que se seleccione, como del estado y 
de la configuracion de la estructura original. 
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Figura 10 Dano en los modelos. Aplicacion de A en el extremo izquierdo 
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4.3 Comportamiento de los porticos de hormigon armado en funcion del 
comportamiento del acero 

Aprovechando que la teorfa serie/paralelo permite asignar un modelo constitutivo a cada 
componente del compuesto, se hace una comparacion del comportamiento de de los porticos 
en hormigon armado, cuando el acero se analiza como un material elastico (EL) y cuando se 
considera elasto-plastico (PL). En la figura 11 se muestra la relacion carga-desplazamiento de 
los porticos de hormigon armado, con los dos tipos de comportamiento del acero. Se observa el 
aumento de resistencia del sistema de acuerdo con la configuracion del refuerzo. En los 
modelos con elasticidad y con plasticidad del acero, se observa la similitud de la rigidez inicial y 
su perdida inicial cuando solo influye el dano en el hormigon. Asimismo, se evidencia la 
influencia de la plastificacion del acero. 




A On) 

-^asc-ci. -w-ws-SL _»_Agcy-s. -o-«gv-ci — Ascv-ei~| 

EL: Modelo elastico PL: Modelo elasto-plastico 

Figura 11 Comparacion de algoritmos. P-A . Aplicacion de A en el extremo izquierdo. 



Tabla 5 Relacion entre reacciones horizontales 


para un desp 


azamiento de 0.0274m 


Portico 


ASF 


AF 


AFC 


AF V 


AFCV 














LLI 

CL 

"'Ll 

CL 

CL 


64.7 


62.7 


68.0 


67.8 


73.9 



En la tabla 5 se muestra la relacion de las reacciones horizontales entre los modelos 
con plasticidad P PL y los casos con elasticidad Pel , fuerzas que son producidas por un 
desplazamiento horizontal de 0.0274m en la viga (maximo desplazamiento con el que se 
pueden comparar los 10 casos). En esta comparacion se observa que, al considerar la 
plasticidad, las cargas laterales representan entre un 62.7% y un 73.9% de las fuerzas de los 
modelos elasticos. Esto implica que, en los modelos con elasticidad en el acero, aunque se 
considere el dano en el hormigon, la capacidad de la estructura se sobreestima entre un 30% y 
un 40% debido a que, en los casos elasticos no se tiene en cuenta el efecto de plasticidad del 
refuerzo. 

De acuerdo con lo observado, en los modelos con elasticidad en el acero, la no linealidad 
del comportamiento se debe unicamente al dano en el hormigon. Por el contrario, en los 
modelos con plasticidad, la no linealidad esta gobernada por el dano en el hormigon y su 
interaccion con la plastificacion del acero. 
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5. CONCLUSIONES 



• La teoria serie/paralelo es una herramienta versatil para analizar el comportamiento de 
materiales compuestos y de sus componentes dentro de una estructura. Por una parte 
permite que cada componente del compuesto se analice utilizando la ecuacion 
constitutiva que sea mas conveniente para predecir su comportamiento (elasticidad, 
plasticidad, daiio etc.) y, por otra, al obtener el comportamiento del material compuesto 
en cada punto de Gauss, lo acopla al resto de la estructura definiendo el 
comportamiento global de la misma. Por tanto, esta herramienta numerica conlleva a 
un mejor analisis de estructuras con materiales compuestos ya que permite considerar 
en la simulacion diferentes tipos de no linealidades. Al mismo tiempo, los costos 
computacionales son reducidos al utilizar un solo mallado continuo para el analisis de 
la estructura, de los materiales compuestos y de sus componentes. 

• El refuerzo con FRP en los porticos de hormigon armado y simple, es una tecnica 
viable para incrementar la capacidad de la estructura. Ademas, siempre que se 
configure adecuadamente el refuerzo en los elementos, se mejora el control de dano y 
se logra aumentar la ductilidad global del sistema. Segun los resultados obtenidos en 
este analisis, tanto en los modelos de hormigon simple, como en los de hormigon 
armado, el aumento de la resistencia global del sistema estructural depende de la 
configuracion de refuerzo que se coloque. Por ello, es necesario que, a traves de la 
simulacion numerica utilizando la teoria de mezclas serie/paralelo, se optimicen los 
analisis de las estructuras con materiales compuestos para mejorar el diseno del 
refuerzo de los edificios y de los puentes que necesitan rehabilitacion o reparacion. 

• Al no tenerse en cuenta la plasticidad del acero en una estructura de hormigon armado, 
no se esta considerando su influencia en la resistencia y en la perdida de rigidez del 
sistema. En estos casos se sobreestima la capacidad estructural, con lo que es posible 
que se disene un refuerzo exterior para unas solicitaciones que no se desarrollaran en 
la estructura por superar su capacidad real. 

Por otra parte, considerar unicamente el fenomeno de dano, tiene su aplicacion en el 
analisis de estructuras construidas unicamente con materiales fragiles, como es el caso 
del hormigon simple, de la mamposteba o de la piedra, entre otros. Esto permite 
estudiar el comportamiento de las estructuras fragiles y su mejora al ser reforzadas con 
FRP. Esto puede tener una interesante aplicacion en la reparacion y rehabilitacion de 
algunas estructuras antiguas. 

• La simulacion numerica, es una herramienta util y de gran aplicacion en el diseno de 
soluciones de refuerzo de estructuras de hormigon armado con materiales compuestos. 
La simulacion numerica con elementos finitos junto con la teoria serie paralelo, permite 
analizar una estructura con diferentes alternativas de refuerzo, involucrando los 
modelos constitutivos de los componentes de los materiales compuestos; asimismo, 
permite que se seleccionen las configuraciones mas eficientes y se determine la 
capacidad de resistencia del sistema estructural reforzado, identificando el nivel de 
daiio que se puede alcanzar. Esta aplicacion es posible a dos escalas: una local, en 
donde se estudian los elementos estructurales, en la que los materiales compuestos se 
utilizan como refuerzo a cortante, a flexion o a compresion; y una escala global, en la 
que se analiza el comportamiento del sistema estructural al reforzar los elementos que 
lo requieran de acuerdo con las normativas existentes. 
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DESEMPENO, RIESGO, RESILENCIA Y SOSTENIBILIDAD 
PERFORMANCE, RISK, RESILIENCE AND SUSTAINABILITY 



William Lobo Quintero 

Universidad de Los Andes, Merida, Venezuela 

RESUMEN 

En este trabajo se integran los conceptos de desempeno, riesgo, resilencia y 
sostenibilidad para ampliar las posibilidades de la ingenieria Estructural y la Ingenieria 
Sismorresistente de intervenir en la Gestion Integral de Los Riesgos. Estos elementos se 
consideran como parte de los sistemas estructurales, pero tambien pueden interpretarse en 
funcion de los sistemas sociales. El desempeno se trata como una concepcion de ingenieria 
estructural o social que predice y evalua con aceptable nivel de confianza demandas y 
capacidades. El riesgo es la esperanza matematica de las perdidas como una convolucion de 
la amenaza y la vulnerabilidad. La resilencia se interpreta como la capacidad de los sistemas 
para mantener su funcion o ser rehabilitados en corto tiempo. La sostenibilidad es el 
compromise de perdurar las infraestructuras manteniendo el equilibrio de los sistemas 
naturales. En el texto se consideran las amenazas multiples, la desagregacion de las 
amenazas, los eventos de baja probabilidad de excedencia, el marco de Hyogo 2007 para la 
reduccion de desastres, los niveles de diseno por desempeno, los metodos de la primera y la 
segunda generacion, las curvas de fragilidad, la cuantificacion de la resilencia, la campana de 
Gauss multidimensional, el calculo de las perdidas estructurales y no estructuales y un 
decalogo para las ciudades sostenibles. 



ABSTRACT 

In this paper the concepts of performance, risk, resilience and sustainability are integrated 
to increase the structural and earthquake engineering possibilities to participate on the Integral 
Risk Management. These elements are considered as a part of the structural systems but also, 
they can be interpreted as a function of the social systems. Performance is considered as an 
structural or social engineering concept that predicts and evaluates demands and capacities 
with acceptable confidence acceptable levels. Risk is the mathematical expectation of the 
losses as a convolution of the hazard and vulnerability. Resilience is interpreted as the system 
capacity to maintain its function or to be rehabilitated in short time. Sustainability is the 
compromise for infrastructures to last while preserving the equilibrium of natural systems in 
order to attain quality of life. In this paper the concepts of multihazards, hazard dessagregation, 
low-probability of exceeding events are considered. The framework of Hyogo 2007 for disaster 
reduction, the levels of performance based design, the methods of first and second generation, 
fragility curves, quantification of resilience, the multidimensional Gauss bell, the computation of 
structural and non structural losses are also reviewed. Finally, a decalogue for sustainability 
cities is proposed. 

1. INTRODUCCION. 

La ocurrencia persistente de los terremotos ha venido generando perdidas crecientes en 
las ciudades y en sus entornos que obligan a considerar modificaciones en los conceptos, en 
las filosofias y en las normativas del diseno sismorresistente, ampliando los campos de 
investigacion y generando mayores expectativas en la sociedad. Desde el terremoto de Mexico 
de 1975 hasta el gran sismo de 1985, aumentaron las perdidas entre 10 y 67 millones de U.S. 
dolares. Luego, entre los terremotos de Loma Prieta (1989), Northridge (1994) y Kobe (1995), 
se contabilizaron perdidas de 182 millones de dolares. Posteriormente, estas cifras se 
incrementaron durante los sismos de Kocaeli, Turquia (11-08-99) a 25 millardos, Chi Chi, 
Taiwan, (21-09-99) a 11.5 millardos y Bhuj, India ( 26-01-01) en 4.6 millardos de dolares. Para 
el esperado terremoto de San Francisco “The Big One”, se esperan perdidas entre 100 y 150 
millardos de dolares y Tokyo: 1200 millardos. El sismo de Italia, 2009, dejo 289 muertes, 1500 
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heridos y 1.000 millones de perdidas economicas y de patrimonio historico. En general, se 
atribuyen algunas razones para justificar estos incrementos: 1. Crecimiento de la poblacion. 2. 
Diseno que solo protege vidas. 3. Utilizacion de zonas desechadas para la construccion. 4. Que 
cada evento es unico y trae sorpresas. En los palses bolivarianos se presentaron terremotos 
importantes entre 1996 y 1999, con un saldo de vlctimas y perdidas materiales resumidas en la 
Tabla 1 y en donde el Terremoto “cafetero” de Armenia, Colombia, fue calificado como un gran 
destructor (Aguiar, 2001). 



Tabla 1. Sismos Bolivarianos (Aguiar, 2001). 



INT. 


PAIS 


FECHA 


SITIO 


MAGNITUD 


PROF. 

(KM) 


MUERTOS 


PERDIDAS 

($) 


VII 


Peru 


12/11/1996 


Nazca 


Md=6.4 


14 


17 


43M 


VIII 


Venezuela 


09/07/1996 


Cariaco 


Ms=6.8 


10 


74 


100M 


VIII 


Bolivia 


22/07/1997 


Cochabamba 


Mw=6.8 


24 


105 


13M 


VIII 


Ecuador 


04/08/1998 


Caraquez 


Ms=7.1 


39 


3 


30M 


VIII 


Colombia 


25/01/1999 


Armenia 


Ms=6.2 


15 


1230(528) 


180M 



Ampliando a otros riesgos probables, una nota del diario El Universal de Caracas de 
fecha 04-09-05, ademas de mostrar el drama humano, estimaba las perdidas del Huracan 
Katrina, 2005 en 100 millardos de dolares, pudiendo las aseguradoras soportar solo entre 20 y 
25 millardos, con un esfuerzo de reconstruccion que podrla durar mas de 10 anos. El caso del 
Huracan Mitch 1998, dejo una estela de desastres en toda la region de America Central: se 
tuvieron seis millardos de dolares en perdidas directas ademas de las perdidas humanas, lo 
que genera una polemica, para quedar los orfgenes en las pobres condiciones economicas, 
sociales y ambientales existentes. Se observo una falta de resilencia economica asociada a la 
pobreza debido a los procesos de exclusion social, falta de organizacion y de participacion 
ciudadana (Lavell, 1989). 

Quizas por el crecimiento demografico, la vulnerabilidad manifiesta, las elevadas cifras 
de pobreza y la desmedida degradacion ambiental, estas catastrofes ocurren ahora con 
mayor frecuencia e intensidad, y surgen como producto de imprevisiones, corto-placismo, 
desconocimiento, imposicion, lucro, olvido, descuido o marginalidad. Es mas, la presencia 
irracional del hombre contribuye a la construccion y reconstruccion de las vulnerabilidades 
mediante la deforestacion, la alteracion de los cursos de agua, el uso de espacios vulnerables, 
la agricultura en laderas, la insistencia en ocupar zonas de comprobada peligrosidad, la 
urbanizacion caotica y la construccion de asentamientos humanos y llneas vitales inseguras. El 
control de los diferentes niveles de danos, ha hecho surgir el diseno basado en desempeho, 
que se desarrolla sobre las variables del riesgo, asegura la resilencia de las infraestructuras y 
de los sistemas ecologicos y sociales, para que sean sostenibles en el tiempo. 

2. LAS DEFINICIONES 

2.1 Desempeno. Es una concepcion de ingenieria estructural que predice y evalua 
con un aceptable nivel de confianza las demandas y capacidades estructurales, para 
asegurar los diferentes niveles de desempeho que se requieran, de acuerdo a las 
caracteristicas de una determinada estructura, de sus componentes y de sus contenidos, 
asegurando ademas la calidad de la construccion y el mantenimiento a largo plazo, para que 
sean realmente construcciones sostenibles. Un diseno basado en desempeho DS-BD aporta 
una metodologla sistematica para determinar la capacidad de desempeho de sistemas 
estructurales o componentes, verificar el desempeho de alternativas, reducir costos o confirmar 
desempenos mayores en edificaciones esenciales. 

2.2 Riesgo. El riesgo se define como la esperanza matematica de la perdida o el grado de 
perdidas esperadas en el espacio de una region y en un tiempo dado, como una convolucion 
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de dos variables: la amenaza H (hazard) y la vulnerabilidad V, ambas con una connotacion 
probabilistica (Esteves, 2007). El riesgo solo es posible si ante una amenaza existen 
condiciones de vulnerabilidad o viceversa, haciendo que ambas variables sean mutuamente 
condicionadas. 

La Gestion Integral de los Riesgos es una nueva vision-accion que orienta la 
planificacion del desarrollo sostenible, porque “toda amenaza se construye socialmente”, y el 
riesgo se obtiene para la concurrencia de determinadas condiciones de vulnerabilidad. Las 
fases de la gestion son la prevencion, la mitigacion, la preparacion, la respuesta y la 
rehabilitacion. Hoy para la construccion de un proyecto importante se exige un “Informe de 
Riesgos” y sabemos que muchos barrios y urbanizaciones necesitan de una evaluacion de los 
riesgos, trabajo que en principio maneja las variables locales. La gestion del riesgo, tiene su 
temporalidad como “Gestion Compensatoria” para superar la vulnerabilidad existente. En el 
nuevo riesgo que la sociedad construira, esta la “Gestion Prospectiva” para estudiar las 
repercusiones futuras, un reto insoslayable e impostergable (Allan Lavell, 1996). 

2.3 Resilencia. Es la habilidad de las unidades fisicas y sociales de mitigar las amenazas, 
contener los efectos de los desastres cuando ocurren y realizar actividades de recuperacion 
que minimicen el impacto social y mitiguen los efectos de los proximos sismos. 

La resilencia sismica se ha definido como la habilidad de un sistema estructural para 
reducir los efectos de un evento sismico, absorberlo si ocurre y recuperarse rapidamente 
despues. Especlficamente, un sistema resilente reduce las probabilidades de falla, aminora las 
consecuencias de falla, en terminos de perdida de vidas, danos y consecuencias negativas 
economicas y sociales, y reduce el tiempo de recuperacion hasta un nivel normal de 
desempeno (Bruneau y Reinhorn, 2003). 

2.4 Sostenibilidad. El desarrollo sostenible es un proceso que va mas alia de la proteccion 
del medio ambiente, y se ocupa de las generaciones actuales y futuras, de la salud, de la 
igualdad entre las personas, de los aspectos sociales y eticos que conlleva el bienestar 
humano. No se puede hacer desarrollo sostenible sin lograr una gestion del riesgo aceptable 
que integre seguridad, educacion, democracia, salud, clima y pobreza. 

Para reafirmar los conceptos de sostenibilidad, el camino es remontarse hasta el Informe 
Brundtland para aceptar una definicion de desarrollo sostenible que “satisfaga nuestras 
necesidades actuales sin comprometer la capacidad de las generaciones futuras para 
satisfacer sus propias necesidades” (Nuestro Futuro Comun, ONU 1987). Aceptando que existe 
un real compromiso para aunar voluntades que puedan crear un futuro sostenible, haciendo 
que las aldeas, los pueblos y principalmente las ciudades se conviertan en verdaderos 
asentamientos humanos, se ampllan las definiciones conocidas, y se aproxima una ciudad 
sostenible la que deber ser equitativa, segura, saludable, ordenada, educativa, ambiental, 
cultural, informada, diversa, gobernable, eficiente y competitiva, que se integran en el sistema 
Ciudad como los subsistemas cultural, ecologica, social y funcional y los elementos son las 
cualidades senaladas arriba (Fig. 1). (Kuroiwa, 2002; Roger, 1998, Lobo-Quintero, 2008). 
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SISTEMA GLOBAL 
Figura 1 La Ciudad como un gran sistema. 



2.5 Relaciones entre las variables. 



Se deduce de las definiciones que: 

2.5.1. No es casual que se logren integrar cuatro condiciones que definen el 
comportamiento de los sistemas estructurales y las infraestructuras: Desempeno, 
Riesgo, Resilencia y Sostenibilidad. En la medida en que podamos tomarlos en cuenta 
estaremos mejorando nuestras capacidades de analisis, diseno, construccion y 
operacion, y asi poder cumplir los requerimientos de los usuarios y el compromiso que 
tienen nuestras profesiones con la sociedad. 

2.5.2. Un diseno basado en desempeno controla los diferentes niveles de 
desempeno y estados limites de los sistemas estructurales y de sus componentes, 
cubrir las amenazas latentes mediante un control adecuado de los danos y de las 
perdidas. 

2.5.3 La gestion integral de los riesgos detecta las amenazas y reduce las 
vulnerabilidades, introduciendo los conocimientos y las tecnicas apropiadas, buscando 
desarrollar la comprension y la concientizacion ciudadana, hacia una mejor calidad de 
vida. 

2.5.4 Evaluar la resilencia de los sistemas estructurales significa saber que estos tienen 
capacidad para recuperarse en un tiempo corto despues de un evento y ademas, que 
puedan guardar parte de la capacidad instalada para utilizarla en eventos distintos. 

2.5.5 Ante la presencia de las diferentes amenazas socio-naturales y antropicas, es 
necesario tener la conciencia profesional de que los sistemas e infraestructuras que se 
construyan, deben ser sostenibles en el tiempo y no cubrir solo la utilidad de las 
presentes generaciones. 

2.5.6 En el Resumen Ejecutivo del Estado del futuro - 2008 ( ONU ), este 
ofrece “condiciones para mejorar en la mayor parte del mundo, pero una serie de 
cuestiones puntuales podrian alterar drasticamente las perspectivas globales. Sin 
embargo, los avances de la ciencia, la tecnologia, la educacion, la economia, y la 
gestion parecieran ser capaces de mejorar mucho el trabajo del mundo en relacion de 
lo que se hace hoy en dia”. 
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3. AMENAZAS MULTIPLES. 

Generalmente cuando las construcciones quedan expuestas a la presencia de diversas 
amenazas, existe en ellas una resilencia inherente, que implica una interrelacion entre las 
diferentes maneras como el sistema responde a las diferentes amenazas. Las ventajas de 
tratar con multi-amenazas estan en el gran potencial que da para un diseno y construccion 
economica, se hace una estimacion mas segura de la resilencia inherente del sistema, un 
tratamiento mas seguro del costo de la vida util, un analisis mas seguro del sistema estructural 
y una optimizacion de la “salud estructural” para incrementar si eficiencia experimental 
(Ettouney, 2007). La salud estructural se mide por instrumentacion, "Tal como un doctor 
monitorea la salud de sus pacientes, quien mide y registra la presion, los ingenieros pueden 
monitorear la salud de los edificios. Los parametros dinamicos de las estructuras permiten 
conocer su estado y su comportamiento, para de esta forma validar modelos teoricos. Este 
analisis da como resultado construcciones y disenos de estructuras mas seguras... y se puede 
hacer por internet", (Boroshek y Gonzalez, 2009). 

Entre las amenazas posibles podran ser terremotos, problemas socio ambientales 
crecientes, deslizamientos, erupciones volcanicas, tsunamis, huracanes, incendios urbanos y 
forestales, uso y manipulacion de gas licuado, terrorismo y/o tecnologia, las cuales se 
recomienda considerar en la evaluacion del riesgo en forma conjunta, en cuenta que 
reduce los costos y aumentan los niveles de seguridad. Tambien, es posible que varias 
amenazas actuen simultaneamente. Cuando se analizan la accion sismica y la carga explosiva, 
se aprecia que las acciones son diferentes en las secciones de la estructura y la consideracion 
de una sola amenaza no suple la accion de la otra tal como se aprecia en la Fig.2. En el 
trabajo del Ettouney et al, 2007, se aprecian los niveles de amplitud y las escalas de frecuencia 
de varias amenazas una grafica bastante util (Figura 3). 



• Rotulas plasticas 





Fig. 2. Explosion vs. Sismo (Ettouney et al, 2007). 
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Figura 3. Distribucion cualitativa de la escala de frecuencias respecto a la Amplitud para 
diferentes amenazas (Ettouney et al, 2007). 

Algunas veces, conviene desagregar las amenazas, para determinar la importancia de 
las fuentes que su curva de amenaza sfsmica se desagrega en diversas curvas referidas a 
cada fuente, tal como se muestra en la Fig. 4. Esto nos permite decir que la mayor amenaza 
para la ciudad, corresponde a la fuente mas cercana (A4) por donde pasa la falla de Bocono, y 
las fuentes A2 y Al contribuyen mas a la amenaza que la A3, especlficamente, para niveles 
bajos de aceleracion. Todas las ciudades deberlan tener su curva de amenaza para precisar 
mejor el riesgo sobre todo cuando se desea manejar un diseno basado en desempeno. 




Aceleracion (g) 

Figura 4. Desagregacion de la amenaza slsmica en Merida (Rivero et al, 2006). 

Realmente las motivaciones para cuantificar el riesgo esta en la posibilidad de las 
nuevas tecnologlas que llevan a nuevos sistemas o a infraestructuras criticas, las demandas 
por desempenos para que cubran mas alia de los mmirnos normativos, la percepcion que el 
riesgo aumenta para las construcciones esenciales, las mayores demandas del colectivo por 
aumentar la seguridad y los graves problemas que generan los eventos con baja probabilidad 
de excedencia, que dejan a quienes manejan el riesgo involucrados en tareas de ensayo y 
error. Para estas amenazas extremas se muestra en la Tabla 2 el orden de magnitud de las 
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frecuencias anuales y los periodos de retorno, tomando en cuenta el maximo evento creible 
MCE (Ellingwood, 2007). 

Tabla 2. Orden de Magnitud de frecuencias anuales y periodos de Retorno en 
Amenazas con baja probabilidad de excedencia ((Ellingwood, 2007). 



N°. 


AMENAZA 


Frecuencia anual 


Periodo de 

Retorno, anos 


1 


Terremoto MCE 


4 x 10-4 


2.500 


Terremoto (2/3 de MCE) 


7 x 10-4 


1.429 


Normativo Canada 


2 x 10-3 


500 


2 


Huracan 


1 x 10 _b 


1.000.000 


3 


Explosion de Gas (por vivienda) 


2 x 10 s 


50.000 


4 


Transporte de gas licuado 


1 x 10 b 


1.000.000 


5 


Explosion de bomba 


2 x 1 0" b 


500.000 


6 


Impacto de avion 


1 x 10 _B 


100.000.000 


7 


Fuego total (por edificio) 


1 x 1 0" 8 /m2 


100.000.000 



En Venezuela, la Norma Covenin 1756-2001 disena para un sismo con una 
probabilidad de excedencia de un 10% con una vida util de 50 anos, un periodo de retorno de 
475 anos que corresponde a una frecuencia anual de 2.1 x 10" 4 , evento que busca la proteccion 
de vidas. Para edificaciones esenciales se tiene un factor de importancia a = 1.3, donde “el 
periodo medio de retorno puede alcanzar mil a mil trescientos anos segun la region”. Tambien 
se informa que los eventos extremos se encuentran asociados a probabilidades de excedencia 
que pueden ser menores que 1% en 50 anos, lo que corresponde a un periodo de retorno de 
4975 anos. En el caso de las instalaciones petroleras, las probabilidades de excedencia 
toleradas son muy pequenas de 1 .5% a 3% en treinta anos de vida util con frecuencias anuales 
entre 5 x 10" 4 y 1 x 10" 3 . Se dan criterios para control de inundaciones en Vargas, Accion del 
viento, aliviaderos de grandes presas y obras excepcionales, los cuales aparecen senalados en 
la grafica de la Fig. 5 como “Riesgos aceptados en los criterios de Diseno”, graficados para 
periodos medios de retorno de los eventos de diseno y periodos medios de retorno de eventos 
extremos, quedando por mencionar otras amenazas de la ingenieria (J. Grases, 2002). 




Figura 5. Riesgos aceptados en los criterios de Diseno (Grases, 2000). 



En el marco de Hyogo, 2007, se asegura que el enfoque de amenazas multiples puede 
mejorar la eficacia: “Las comunidades suelen ser expuestas a una variedad de riesgos debido a 
multiples amenazas, las que pueden ser de origen natural o inducidas por el hombre tales 
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como las amenazas hidrometeorologicas, geologicas, biologicas, tecnologicas o 
medioambientales. Los riesgos acumulados no pueden ser abordados eficazmente si 
simplemente se consideran algunas amenazas. Un enfoque de riesgos multiples considera 
vincular los conocimientos de las distintas amenazas en la gestion de desastres y de riesgos de 
desastres, en las estrategias pollticas, evaluaciones profesionales y analisis tecnicos, en las 
capacidades operativas y en las de comprension del publico, que conduzca a una mejor 
efectividad y mayor costo-eficiencia” (A Guide for implementing the Hyogo Framework 
International Strategy for Disaster Reduction, 2007). 

4. DISENO BASADO EN DESEMPENO. 

La version de FEMA 445- NEHRP nos permite delinear un proceso actualizado de diseno 
basado en desempeno, dada la amenaza potencial y considerando las incertidumbres al riesgo 
y a la respuesta estructural, aceptando que es un proceso iterativo en el diseno preliminar 
hasta que se logre el desempeno adecuado. La asociacion de un nivel de desempeno (estado 
de danos) a un novel de amenaza se llama objetivo de desempeno. Los niveles de desempeno 
se identifican mejor sobre una curva de capacidad, desde la publicacion del documento 
SEAOC, Vision 2000 (Fig. 6), que quedan especificados en la Tabla 3 (Gutierrez, 2007). 




Figura 6. Niveles de desempeno segun SEAOC, Vision 2000 (Gutierrez, 2007). 



.Tabla 3. Niveles de desempeno. SEAOC, Vision 2000, (Gutierrez, 2007). 



Nivel de 
Desempeno 
Estructural 


Descripcion 

cualitativa 


Desplazamiento 
inelastico (%) 


Desempeno de 
componentes 
no 

estructurales 


Dano (%) 


SP-1 


Operacional 


0 


NP-1 


0-10 


SP-2 


Ocupable 


30 


NP-2 


5-30 


SP-3 


Proteccion de 
vidas 


69 


NP-3 


20-50 


SP-4 


Cercano al 
derrumbe 


80 


NP-4 


40-80 


SP-5 


Derrumbe 


100 


NP-5 


> 70 
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Los seis aspectos recomendados para su estudio en Vision 2000, han venido 
marcando el rumbo de la investigacion en esta area, en orden decreciente de sofisticacion: 

A). El enfoque de un diseno completo. B). El diseno basado en desplazamientos. C). Un 
enfoque de diseno basado en Energia. D). Enfoque de diseno general fuerza/resistencia. E). 
Enfoque simplificado fuerza/resistencia. F). El enfoque de diseno prescriptivo. El diseno 
definitivo hace cumplir que cada objetivo de desempeno corresponda a un riesgo aceptable y a 
un nivel de dafios estructurales y/o no estructurales, asi como de las perdidas 
correspondientes. All! esta el corazon de. Metodo DS-BD, que podra utilizarse para: 

• Disenar edificaciones con mayor nivel de confianza que el establecido por las normas. 

• Disenar edificaciones capaces de cumplir las normas presentes pero a costos bajos. 

• Disenar edificaciones para alcanzar el mayor desempeno (perdidas mlnimas) 
aceptando las normativas presentes 

• Disenar edificaciones que esten fuera de las prescripciones de las normas, respecto al 
sistema estructural, la configuracion y los materiales, pero ajustandose a la normativa. 

• Evaluar el potencial sismorresistente de las estructuras existentes y estimar el potencial 
de perdidas probables ante un evento sismico. 

• Evaluar el potencial prescriptivo de una norma y servir de instrumento para mejorar los 
criterios de diseno y edificaciones mas consistentes y confiables. 

Es conveniente acompanar la Tabla 3 con las propiedades de excedencia, los periodos de 
retorno y la frecuencia anual que corresponda a los niveles de desempeno, para definir los 
riesgos correspondientes. En la Tabla 4, para una vida util de 50 anos, se amplfa la presentada 
por Gutierrez, 2007, proveniente de Rose School, materia que debe discutirse hacia una nueva 
norma. 



Tabla 4. Niveles de Diseno, Probabilidades de Excedencia, Periodos de retorno 



(Gutierrez, 2007) 



Nivel 

de 

Des. 


Estructuras normales 


Estructuras esenciales 


Estructuras con contenidos 
peligrosos 


Prob. 

exc. 


Tr 

anos 


Free. 

anual 


Prob 

exc. 


Tr 

anos 


Free. 

anual 


Prob. 

exc. 


Tr 

anos 


Free. 

anual 


ND-1 


50% 


73 


0.013767 


20% 


224.6 


0.004453 


20% 


224.6 


0.004453 


ND-2 


10% 


475 


0.002105 


4% 


1225 


0.000816 


2% 


2475 


0.000404 


ND-3 


2% 


2475 


0.000404 


1% 


4975 


0.000201 


1% 


4975 


0.000201 



La determinacion del desempeno se puede hacer interceptando las curvas de 
capacidad con las demandas espectrales, dentro de los que se han llamado metodos basados 
en desempeno de la primera generacion. En Venezuela, no se ha generado experiencia 
suficiente con estos metodos, pero deberlan ser incorporados en las normas sismorresistentes. 

Una muestra no exhaustiva de los metodos revelan la forma como se ha adelantado el 
conocimiento; mayor informacion la da el Prof. Sigmund Freeman en Internet cuya direccion se 
da en las referencias. En la fig. 7(a), se muestra su Metodo propuesto en 1978, y basado en 
espectros amortiguados que con porcentajes de dafios y que sirvio de base a los desarrollos 
posteriores. En 1998, Freeman propuso su version acerca del espectro de capacidad. En la 
Fig. 7(b), esta ilustrado la determinacion del punto de desempeno, mediante el Metodo del 
Espectro de Capacidad propuesto por en el ATC 40, un esfuerzo importante que dejo huellas 
en la ingenierla sismorresistente. Este metodo determina un perlodo equivalente T eq . y un 
factor de amortiguamiento (3 eff . como funcion de la ductilidad p, el factor de endurecimiento a y 
un factor k cercano a uno, segun las expresiones (1). 
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Eip«ctro 





Fig. 7. (a). Metodo de Freeman et al, 1978. (b).Punto de desempeno (ATC 40,1996). 



( 1 ) 



Otras de los metodos conocidos, son el de Aschheim y Black, 2000 mostrado en la Fig. 
8 para sistemas elasto-plasticos, en el piano (Cy, Ay), referido a la ductilidad p cumpliendo la 
regia de los desplazamientos iguales R = p. La capacidad portante de la estructura se corta 
con el cumplimiento de la ductilidad y el metodo ha tenido relativo exito. No hay duda que estas 
metodologlas han introducido una simplicidad, una reduccion de calculo y una forma grafica 
de visualizar el problema de demandas y capacidades incluso mas alia del comportamiento 
elastico, a pesar de las limitaciones de las curvas pushover. Estas curvas son de la primera 
generacion y se ajustan al modo fundamental, haciendo mejor estimacion en 
edificaciones bajas o medianas. En tal caso, se han ensayado soluciones para ajustar 
mejor la demanda y mejorar las condiciones del espectro de capacidad, hacia la consideracion 
de los modos superiores y el uso del desplazamiento cedente (Kalkan y Kunnath, 2004; 
Aschheim y Black, 2000). 




Figura 8. Metodo de Aschheim y Black, 2000 
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4.1. Analisis pushover modal (APM). 

Hoy se aplica la propuesta del APM (Analisis Pushover Modal) para tomar en cuenta los 
modos superiores, particularmente en el rango inelastico, a pesar de que la superposicion 
modal no es valida, ya que los modos no estan siempre desacoplados, pero, se ha demostrado 
que el error cometido es muy pequeno si se compara con la historia de aceleraciones (Chopra, 
2004). Este metodo consiste en resolver el analisis elastico modal de la estructura y para 
realizar una curva pushover estatica, calcular las fuerzas horizontales efectivas asociadas a un 
registro acelerografico que se incrementa y permite registrar los deplazamientos modales 
correspondientes. Posteriormente, se transforma cada curva modal en una bilineal que se 
convierte en un espectro de capacidad modal. Los valores definitivos se determinan aplicando 
la raiz cuadrada de la suma de cuadrados incrementada por el efecto de las cargas 
gravitatorias. El uso posterior es la intercepcion con los espectros inelasticos (Chopra y Goel, 
2000 y 2001). El metodo APM es mas seguro que los metodos FEMA para estimar las 
demandas sismicas. La metodologia de transicion hacia los nuevos metodos del Prof. Chopra, 
2001, es la siguiente: 



1 . Hacer el analisis dinamico modal y determinar los parametros para el modo n: 

T n (periodo), cj n (frecuencia circular), cjj n (modo), T n (factor de participacion modal), L n , 
(factor de masa). Otros simbolos son la matriz de masas m y la masa efectiva modal 

M* n = C}D T n .m. C}j n 

2. Calcular: S n = T n .m. <£„ ; S = Z S n (n = 1,...n) ; S* n =m. 4> n 

Peff,n(t) = - s n iig (t) = -r n .m. c}) n = - m 1 u g (t) ; 1 = vector unitario (2) 

3. Efectuar la curva pushover para cada modo n en terminos del corte basal V bn y el 
desplazamiento del piso r en el modo n: u rn 



4. Idealizar la curva obtenida mediante una bilineal, encontrando el punto V bny , u rny . 



5. Convertir la curva bilineal (V bn - u rn ) en un espectro de capacidad (F sn / L n - D n ), 
aplicando. ( F sny /L n — V bny / M n ) y ( D ny — u rny / r n . cf) rn ), 

6. Calcular: el periodo de vibracion elastico, el desplazamiento en el piso r y la respuesta 
total incluida la respuesta bajo carga gravitatoria r g . 




Fig. 9. (a). Tres periodos y formas modales. (b). Distribucion de fuerzas S* n para los 

tres modos de un edificio para SAC en 
Los Angeles (Chopra, 2007). 



( 3 ) 

Una ilustracion de este metodo esta muy bien detallada en la referenda Chopra, 2007, 
de la cual se toman solo algunos detalles de un edificio de 9 pisos para SAC en Los Angeles, 
con efecto P-delta (Fig. 9 y 10). Segun Chopra, 2007 
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a) . Este metodo mantiene la simplicidad conceptual de los procedimientos basados en 
pushover. 

b) . Cuando son validos solo los tres primeros modos, el metodo requiere un calculo comparable 
al de FEMA-356. 

c) . Este es un metodo aproximado que no da una estimacion tan exacta con el analisis historia- 
tiempo de aceleraciones. 

d) . Con este metodo se pueden mejorar los procedimientos senalados como de la primera 
generacion. 

e) . Se han dado soluciones de APM para estructuras asimetricas (Goel R. y Chopra A. K, 2005) 
y para puentes (Kappos et al, 2005). 



Los Angeles 




Roof Disp acement' Hei ght (%■} 

"Made" 1 

9-Story "Mode" 2 

"Made" 3 

Figura 10. Curvas pushover modales con efecto P-delta. 

4.2. Metodos de diseno por desplazamiento. 

Entre los pioneros en considerar el diseno incorporando los espectros de desplazamiento 
fue Moehle, 1992, dejando un metodo sencillo para controlar las derivas y determinar los 
detalles estructurales y no estructurales que aseguren un desempeno adecuado. 

Posteriormente se presento el “Diseno Directo Basado en Desplazamientos” controlando 
dos estados Ifmites: operativo y de danos (Priestley, 2000). Ultimamente, se ha publicado un 
texto (Displacement-based design of structures DB-SD) cuya filosofla de diseno se basa en la 
determinacion de una resistencia estructural optima para alcanzar el estado Ifmite de 
desempeno, referido a un nivel de danos bajo una intensidad sismica. Tambien se hace enfasis 
en la distribucion de la resistencia sobre la estructura, mediante metodos de analisis estructural 
y metodos de diseno por capacidad. Se dejan las herramientas del DB-SD mediante tres 
programas: Analisis de momento-curvatura (Cumbia), Analisis historia-tiempo (Ruaumoko) y 
Analisis dinamico con elementos de fibra (Seismo Struct) (Priestley et al, 2007). Este grupo ha 
revisado los metodos de diseno por desplazamiento, que son tambien los metodos ya 
mencionados del diseno basado en desempeno. Para que los metodos sean valiosos para 
cubrir este procedimiento, se desea que tengan: 

• Simplicidad. Que puedan aplicarse relativamente facil y aceptados por la comunidad 
de disenadores. 

• Versatilidad. Que sean capaces para disehar la variedad de formas estructurales para 
el rango de niveles de desempeno. 

• Desempeno. Que sean capaces de disehar eficientemente una estructura para 
deformarse dentro de los parametros establecidos. 

• Completitud. Que proporcionen suficiente informacion al disenador, para que pueda 
proceder sin hacer demasiadas hipotesis, que luego son dificiles de interpretar 
(Sullivan, 2002) 



Desempeno, riesgo, resilencia y sostenibilidad 



67 



En Rose School, Universidad de Pavia, ocho metodos propuestos se aplicaron a 5 formas 
estructurales, con cambios en las resistencias de disefio, poca influencia en el desempeno 
y limitaciones que pueden ser resueltas. Quedo claro que disehar para EQ-I = SP-1 (Fig. 7) 
para modelos no agrietados o secciones gruesas, por su rigidez controla el diseno y mantiene 
la deriva para EQ-IV = SP- 4; que cuando la deriva controla el desplazamiento, reduce el 
desplazamiento cedente y aumenta la ductilidad y que trabajando con el espectro inelastico 
con buena estimacion del punto cedente, se puede disehar para varios estados limites Se 
concluye que el futuro para los metodos de diseno basados en desempeno es brillante. 
Los cuatro valores mencionados arriba y aplicados a los metodos, se muestran en la Fig. 11 
(Sullivan, 2002). 



Evaluation tie los metotlos tie Diseno basatlos en Desplazamiento 




o Simplkitlad ■ Versatilidad □ Desempeno t Claiidad f oinpletitud 



Figura 11. Propiedades de los metodos basados en desempeno (Sullivan, 2002). 

4.3. Metodos basados en Energia. 

Un metodo desarrollado recientemente basado en los conceptos de balance de 
energia, llamado “Diseno plastico basado en desempeno-DPBD”, convierte la curva de 
capacidad estructural en otra de capacidad de energia, para resolver favorablemente sistemas 
elasto-plasticos de acero simples y de multiples grados de libertad. Se parte de que la energia 
calculada de la respuesta monotonica de un sistema inelastico es igual a la energia del sistema 
elastico E correspondiente, donde M es la masa del sistema, S v es la pseudo-velocidad, V y es 
el cortante cedente, D y es el desplazamiento cedente y D m es el desplazamiento maximo. 
(Leelataviwat et al, 2008). 

Para sistemas simples la energia vale (4) 

Se toma un factor y para ampliar el balance de energia a todos los periodos posibles, 
la ecuacion de energia queda modificada y y es funcion de la ductilidad p y del factor de 
respuesta Ry: 

yE = r^MS; = Ir,Z>, + V,(D a - D, ) 

( 5 ) 



y.ltLzl 

x; 



( 6 ) 
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Sistema inelastico 

Sistema conespondiente 

elastico 




Fig. 12. Concepto de balance de Energia modificado (Leelataviwat et al, 2008). 




Fig. 13. (a). Balance de energia en modo n. (b). Oscilador equivalente. (Leelataviwat et 
al, 2008 ). 




Fig. 14. Diseno plastico basado en desempeno - DPBD: a). Curva pushover, b). Energia- 
desplazamiento. c). Demanda de energia. d). Demanda de desplazamiento (Leelataviwat et al, 
2008). 



Estas energias se representan en la Fig. 12. Cuando el balance de Energia se extiende a 
multiples grados, se usa el concepto del oscilador equivalente modal, conocida la forma del 
modo se obtiene la masa modal M*. Mediante la superposicion modal se podra obtener la curva 
pushover, segun el balance de energia modal mostrado en las Fig. 13 y 14. 

En la Fig. 14, para el modo n, se muestra la curva pushover elasto plastica, el diagrama 
energia desplazamiento, el diagrama de demanda de energia y el punto de intercepcion o de 
desempeno que determina la capacidad del sistema. Esta propuesta puede aplicarse al analisis 
de sistemas simples y multiples, y fue calibrada con el metodo APM, aunque no esta 
totalmente acabada. La ventaja esta en el uso de la energia, es una representacion global que 
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evita la relacion entre corte basal y desplazamiento en el tope usada en el espectro de 
capacidad. 

4.3. Analisis Dinamico Incremental (I.D.A.) 

Es un metodo que se ha creado como herramienta para el Diseno Sismico Basado en 
Desempeno (DS-BD), que ofrece la prediccion de la demanda y de la capacidad 
sismorresistente, utilizando uno o mas registros acelerograficos escalados a diferentes niveles 
de intensidad, produciendo curvas de respuesta en cada piso y para cada nivel de intensidad. 
Sabemos que el analisis estatico estructural ha evolucionado hacia el analisis dinamico, el 
estatico no lineal, el de historia-tiempo de aceleraciones y finalmente, al analisis dinamico no 
lineal incremental IDA (Vamvatsikos y Cornell, 2002). 

Este metodo cuando se presento tenia multiples propositos, que por sus potencialidades, 
se han superado en el tiempo, como es el desarrollo de la investigacion en las siguientes 
areas: 

* Entender el rango de demandas respecto al rango potencial de niveles de un registro 
sismico. 

* Comprender las implicaciones estructurales de los movimientos sismicos raros o severos. 

* Comprender los cambios en la respuesta en la medida que el movimiento aumenta. 

* Producir estimaciones a la capacidad dinamica de un sistema global. 

* Ver que tan estables o variables son los parametros de un multi registro IDA. 

CTJRVA IDA SIMPLE Y CURVA PUSHOVER 




Fig. 15. Curva IDA y Curva Pushover de un conjunto de tres niveles de aceleracion 
(Vamvatsikos y Cornell, 2002). 

Una curva IDA se muestra en la Fig. 15, para el Terremoto de El Centro de 1940 en 
una estructura de acero de 20 pisos, con periodo fundamental de 4 seg. incluyendo el efecto 
P-delta, para tres niveles de aceleracion: 0.01 g, O.lg y 0.2g. Como comparacion se muestra la 
curva pushover estatica respecto a primer modo (Vamvatsikos y Cornell, 2002). Los valores de 
las derivas por nivel de aceleracion se observan en la Fig. 16 y 17 son suficientemente 
ilustrativos. Igualmente en la Fig. 18 se muestra la determinacion de la mediana, el 16 y el 84 
percentil de todas las curvas (Krawinkler y Miranda, 2004). 



70 



William Lobo Quintero 



Relation <le la deriva inter-pisos v los niveles de cada piso 




Relation de deriva inter pisos O- 

Fig. 16. Curvas de respuesta por piso para cada incremento de la aceleracion (Vamvatsikos y 

Cornell, 2002). 
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Relation de deiiva inaxiina inter pisos 
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Fig. 17. Muestra de curvas IDA para cada entrepiso (Vamvatsikos y Cornell, 2002). 
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Fig. 18. Aplicacion de curvas IDA para determinar la mediana y las respuestas de 16 y 
84 percentil. T1 =0.9 seg. (Disertacion de Medina, 2002; Krawinkler y Miranda, 2004). 

4.4.Curvas de Fragilidad. 

Una curva de fragilidad describe la probabilidad (P ik ) de alcanzar o exceder un estado de 
danos (d,) en una construccion para varios niveles determinlsticos del movimiento de tierra Y k : 

P ik = P [ D > di] / Y = y k ] ( 7) 




Siendo D (variable aleatoria o vector de danos) = ( d 0 , di, d 2 ,d 3 ,... d n ) 

Es necesario modelar las incertidumbres asociadas con las capacidades estructurales y 
sus demandas. En principio las curvas de fragilidad se pueden obtener mediante la ecuacion 
fundamental de las probabilidades falla expresada en funcion de las acciones S y las 
resistencias R, tomando en consideracion la funcion de probabilidades de S (F s (r)) y la funcion 
de densidad de probabilidades de f R (r) : 



( 8 ) 



Para una funcion normal estandarizada, donde son conocidas las medias p R y p s y 
las desviaciones estandar cr R y a s esta probabilidad de falla p f se expresa en la ecuacion (9), 
donde p es el Indice de confiabilidad. La ecuacion similar (10) trabaja con el logaritmo 
neperiano del cociente de las medias p R y p s con las desviaciones estandar Pr y Ps, y la 
ecuacion (11) se refiere a la probabilidad de danos d s correspondientes a un pseudo 
desplazamiento S d como la funcion de distribucion acumulativa lognormal dependiente de la 
desviacion estandar p c y del valor medio de los pseudo desplazamientos S d ds . 



P [ D > <li / Si ] 
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Fig. 19. Curvas de Fragilidad para edificios bajos (Singhai y Kiremidjian, 1996) 

En la Fig. 19 se muestran las curvas de fragilidad obtenidas para cuatro estados de 
danos, segun el indice de Park y Hang: Menor (0.1 - 0.2), Moderado (0.2 - 0.5), Severo (0.5 - 
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1.0) y Colapso (> 1). Los datos se manejan con el metodo de Montecarlo o con las reglas del 
hipercubo latino. Con los valores obtenidos se crean las matrices de probabilidad de danos. 
Algunas veces se interpretan las curvas de fragilidad como areas referidas a los danos, tal 




Fig. 20. Curvas y Areas de Fragilidad (Hazus, Fema, 1999). 




Fig. 21. Determinacion de Danos estructurales y no estructurales (Hazus, Fema 1999). 



La tendencia actual es obtener funciones de fragilidad cuando la funcionalidad y la 
seguridad se describen con parametros multiples. Generalmente para las edificaciones se 
utilizan dos variables: la aceleracion pico efectiva APE y la deriva inter-pisos, aunque tambien 
se toma la probabilidad de excedencia. La clave esta en la correcta evaluacion de los estados 
limites, teniendo presente parametros demasiado influyentes como el amortiguamiento, la 
rigidez, el modelo de comportamiento histeretico y considerar las variables aleatorias de multi- 
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respuesta, especialmente cuando tienen mucha importancia las componentes no estructurales. 
Para este caso hay la necesidad de determinar cada estado limite a traves de analisis y 
experimento. Sin esta cuantificacion los estados limites de fragilidad no reflejan con seguridad 
el desempeno de una estructura y las decisiones que se tomen pueden ser incorrectas 
(Cimerallo et al, 2006). Como parte final en la Fig. 21 se muestra la aplicacion de las curvas de 
fragilidad en Sd y Sa, para la determinacion de los danos para cinco estados limites de danos 
(Hazus, Fema, 1999). 

5. EVALUACION DEL RIESGO. 

Con la determinacion de las curvas de fragilidad se tiene un mecanismo para evaluar el 
riesgo antes y despues de un evento, llamese geologico, sismico, geomorfico, atmosferico, 
ecologico, contaminante, toxico, por incendio u otra actividad humana. La Ley de Gestion 
Integral de los Riesgos decretada por la Asamblea Nacional el 08-01-2009 debera ejecutarse 
en todas sus facetas y de la misma se deducen acciones que debe participar la ingenieria, las 
cuales se citan algunas importantes: 

• Esta Ley se circunscribe a los riesgos de caracter socionatural y tecnologico, originados 
por la probabilidad de ocurrencia de fenomenos naturales o accidentes tecnologicos 
potenciados por la accion humana que puedan generar danos sobre la poblacion y la 
calidad del ambiente. 

• El ente rector del Sistema Nacional de Ciencia y Tecnologia, promovera la generacion 
de conocimientos relativos a los niveles de amenaza, vulnerabilidad y riesgos en los 
distintos espacios geograficos y el libre acceso a dicha informacion. 

• Todos los proyectos para obras de infraestructura deberan contemplar criterios de 
reduccion de riesgos a fin de garantizar la preservacion de la poblacion y la 
sustentabilidad de dichas inversiones. 

• Los organos de administracion de desastres coordinaran la ejecucion de las acciones 
de rehabilitacion de servicios basicos en el menortiempo posible. 

• Promover el desarrollo de la normativa requerida para la instrumentacion de la politica 
nacional de gestion integral de riesgos socionaturales y tecnologicos. 

En el 2005 se celebro en Hyogo, Japon, la Conferencia Mundial sobre Reduccion de 
desastres en donde 168 gobiernos aceptaron un Plan de 10 anos para lograr un mundo mas 
seguro frente a las amenazas naturales (La gestion del Riesgo de Hoy, ONU, 2005). El 
llamado “Marco de Hyogo” establecio los siguientes objetivos estrategicos: 1). La integracion de 
la reduccion del riesgo de desastres en las politicas y la planificacion del desarrollo sostenible. 

2). El desarrollo y fortalecimiento de las instituciones, mecanismos y capacidades para 
aumentar la resilencia ante las amenazas. 3). La incorporacion sistematica de los enfoques de 
la reduccion del riesgo en la implementacion de programas de preparacion, atencion y 
recuperacion de emergencias. 

En la Fig. 22 se muestra la interrelacion de los conceptos de Amenaza, Vulnerabilidad, 
Riesgo y desastres, dentro de un ciclo continuo, en el cual estos se generan mutuamente y 
son fenomenos de interaccion en un territorio concreto de la dinamica de la naturaleza con la 
dinamica de la comunidad (La Gestion de Riesgo de Desastres 2008, EIRD-ONU, El Marco de 
Accion de Hyogo). 
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Fig. 22. Manejo del riesgo como interaccion entre comunidad y naturaleza (El Marco de 

Accion de Hyogo, 2008). 

5. CUANTIFICACION DE LA RESILENCIA ESTIMACION DE PERDIDAS. 

Para la gestion de los riesgos, la Resilencia utiliza las siguientes variables: a. Tecnica: 
habilidad de los sistemas fisicos para desempenar niveles aceptables. b. Organizativa: 
capacidad y responsabilidad para tomar acciones, funciones y decisiones. c. Social: evitar 
las consecuencias negativas de la perdida de los servicios criticos. d. Economica: 
capacidad para reducir las perdidas directas e indirectas (Todas constituyen el T.O.S.E). La 
resilencia R se representa en el area de un diagrama (Q, t) (Fig. 23) que muestra el 
instante ti cuando ocurre un evento, la reduccion de la funcionalidad Q a causa de la 
vulnerabilidad, la ordenada restante es la Robustez y la Resilencia es la recuperacion 
ocurrida en el tiempo (t 2 - L) o area sombreada; Se toma a Qo como el 100%. La 
estructura resilente se recupera en el menor tiempo posible, pues revela el mayor grado de 
preparacion previa (Buckle y Lee, 2006). 




Fig. 23. Representacion de la Resilencia (Buckle and Lee, 2006). 

Especificamente, un sistema es resilente cuando: 

a) . Reduce las probabilidades de falla. 

b) . Reduce las consecuencias de la falla, en terminos de perdidas de vidas, danos y 
consecuencias sociales y economicas. 

c) . Reduce el tiempo de recuperacion a su nivel normal de desempeno. 

Las propiedades que refuerzan la resilencia, tanto en las areas fisicas y sociales, 
mostradas en la Fig. 24, son las siguientes (Bruneau y Reonhorn, 2004): 
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*. Robustez. Resistencia o habilidad de los elementos, sistemas y otras unidades de analisis 
para recomendar un nivel de esfuerzos o demanda sin experimentar una degradacion o perdida 
de funciones. 

*. Redundancia. Extension de elementos, sistemas y otros elementos y otras unidades de 
analisis que son sustituibles y pueden satisfacer requerimientos funcionales en caso de 
interrupcion, degradacion o perdida de funcionalidad. 

*. Recursos. Capacidad de identificar problemas, establecer prioridades y movilizar recursos 
cuando existen condiciones que perturban algun elemento, sistema o unidad de analisis. 

*. Rapidez. Capacidad de cumplir prioridades y alcanzar metas en tiempos minimos. 




Rapidez 

Fig. 24. Esquema de propiedades de la Resilencia (Bruneau y Reinhorn, 2004). 




Fig. 25. Resilencia desacoplada (Cimellaro et al, 2006). 

La definicion grafica de la resilencia para diversos eventos, donde la funcionalidad Q(t) 
se expresa en forma adimensional y permite proponer una expresion para evaluar 
matematicamente la resilencia (ecuacion 12) como un proceso estocastico que en cada 
proceso ensambla una funcion continua por pasos, tal como se muestra en la Fig. 25. 
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Donde NE es el numero de eventos extremos esperados durante la vida util del sistema 
T lc ; N, es el numero de eventos extremos distintos; T RE es el tiempo de recuperacion del 
evento E; t 0 E es el tiempo de ocurrencia del evento E; L(l, T RE ) es la funcion normalizada de 
perdidas; f REC (t, t 0 E, T RE ) es la funcion de recuperacion; P (I) es la probabilidad de que un 
evento (I) de intensidad dada ocurra en el intervalo de tiempo Tlc; Pe (0,Tlc) es la probabilidad 
de que un evento ocurra E veces en un intervalo de tiempo T LC ; a R es el factor de recuperacion 
y H (to) es la funcion involucrada. En la ecuacion hay una funcion de perdidas L (I,Tre), la 
funcion de recuperacion f RE (t,t 0 E, T RE ) y la funcion de fragilidad que no aparece explicita ( 
Cimellaro et al, 2006 ). 



Fig. 26. Representacion 3D de una funcion de densidad de probabilidad f R Q 

(Nowak y Collins, 2000). 



Si se tienen las variables Q que representan a las cargas y R a las resistencias, en un 
espacio bi-dimensional (Fig. 26a), se puede separar el dominio de la seguridad del dominio de 
la falla, y la frontera es una linea recta a 45° que parte del origen R-Q = Oy cada curva 
representa una probabilidad de excedencia. Como R y Q son variables aleatorias se puede 
definir la funcion de densidad de probabilidades f RQ (r,q) representada en el espacio 3D (Fig. 
26b) que en forma de una campana multidimensional de Gauss muestra mejor la densidad 
general de probabilidades (Nowak y Collins, 2000). 

Las curvas de recuperacion son funciones lineales, exponenciales y trigonometricas. La 
lineal se aplica cuando no hay suficiente informacion, la exponencial es para una respuesta 
rapida y organizada de la sociedad y la trigonometrica se refiere a un evento de baja 
probabilidad de excedencia sin ninguna preparacion. En la Fig 27 se muestran las funciones de 
recuperacion (Cimellaro et al, 2006). En la cuantificacion de la resilencia estructual y de los 
elementos no estructurales, se puede aplicar tres estados llmites tal como se muestran en la 
Fig. 28, donde a la Integridad o Funcionalidad, se aplican un estado llmite de servicio con 
perdidas pequenas (1), un umbral de colapso con perdidas maximas (2) y un estado intermedio 
de danos (3), todos con niveles aleatorios asi como las curvas de recuperacion. 




R 
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Fig. 27. Funciones de recuperacion de la Resilencia (Cimellaro et al, 2006). 





Fig. 28. Estados limites aleatorios (Bruneau y Reinhorn, 2004). 

La Campana multidimensional (Fig. 29) permite delimitar los espacios limites, encontrar 
las areas seguras en el espacio de aceleraciones pico Psa y desplazamientos o derivas pico 
Psd, as! como las areas crlticas cuando se cortan con los contornos isoprobabillsticos. Aqul se 
toman limites determinlsticos lo cual no es necesariamente asl. La probabilidad que excede la 
respuesta para un estado llmite se puede obtener directamente del volumen ubicado por 
debajo de la distribucion que excede el llmite especlfico. Las condiciones son diferentes para la 
estructura como para los elementos no-estructurales. 

BRUNEAU y REINHORN, 2004, MCEER 



ESTADO LEUITE I)E VELOC'IDAD 




Fig. 29. Uso de la campana multidimensional para controlar los limites de las respuestas 

(Bruneau y Reinhorn, 2004). 
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En la Fig. 30 se utiliza la resilencia para calcular las perdidas probables en elementos 
no estructurales determinando la probabilidad de exceder un estado limite PEL de una curva de 
fragilidad determinada para amenazas variables expresadas en probabilidad de excedencia 
(Bruneau y Reinhorn, 2004). 
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Fig. 30. Calculo de Perdidas en elementos no estructurales (Bruneau y Reinhorn, 2004). 
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fig. 31. Respuesta Sismorresistente. a). Curva de resilencia. b). Nueva curva de resilencia en 
estructura no reparada c). Probabilidad de perdidas estructurales antes de un sismo. d). 
Probabilidad de perdidas al repetir un sismo (Bruneau y Reinhorn, 2004). 



En la Fig. 31, se parte de la curva de resilencia y de la probabilidad estimada de 
perdidas antes de un sismo que en la curva de fragilidad tiene una probabilidad de excedencia 
de 10% en 50 anos de vida util. Si repitiera el sismo en el tiempo t E+ , la probabilidad de exceder 
el estado limite sera mayor y una perdida ulterior podria ocurrir, tal como lo muestra la segunda 
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curva de resilencia. Finalmente, cuando se hace una rehabilitation estructural, la resilencia 
puede cuantificarse aceptando el deslizamiento de las curvas de fragilidad para reducir las 
perdidas tal como lo expresan los diagramas de la Fig. 32. Se observa la rehabilitation pre- 
sismica y el reforzamiento de la curva de resilencia para reducir la probabilidad de perdidas 
(Bruneau y Reinhorn, 2004). 





Fig. 32. Rehabilitacion pre-sismica y deslizamiento de curvas de fragilidad 
(Bruneau y Reinhorn, 2204). 



6. CIUDADES SOSTENIBLES. 

Segun las recomendaciones de la O.N.U. establecidas en “El Marco de Accion de Hyogo, 
EIRD, 2008” la reduccion de los riesgos debe incorporarse a las actividades del desarrollo 
sostenible y por tanto los ingenieros deberiamos estar atentos a participar en los retos que 
implica la sostenibilidad. En Merida (Fig. 33) hemos propuesto “El Decalogo de la Merida 
Sostenible” que deja una doctrina de accion (Lobo Quintero, 2007). 

1°. Acordar un nombre y mantenerlo. Promover el proyecto “Merida Sostenible, una Ciudad 
para la Gente" con precisiones conceptuales, que permitan satisfacer nuestras necesidades sin 
comprometer la capacidad de las generaciones futuras. Sera inconveniente personalizarlo o 
presentarlo bajo otra denominacion como “Merida Sustentabte" , que implica otras ideas y 
principios. Cuando los proyectos generan dudas en su definicion o tienen metas ambiguas, no 
logran aplicaciones eficaces. 

2°. Adoptar la ideologia ambiental. Adoptar el “desarrollo sostenible”, como la politica de las 
polfticas, como una ideologia ambiental, plataforma comun para hacerla cumplir e irla viendo 
en sus avances, con horizontes en el mediano plazo y resultados innovadores parciales, que 
vayan creando el consenso y la participacion. El conocimiento ambiental debe ser difundido y 
practicado. 
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Fig. 33. Vista de Merida con un desarrollo urbano no ordenado. 



3°. Definir el paradigma. La sostenibilidad de la ciudad de Merida depende del uso mesurado 
de su capital natural, contar con un proyecto de ciudad cultural, cientifica, tecnologica, turistica 
y ambiental, que dependa de la posibilidad de abastecerse de suministros para subsistir y de 
aprender a deshacerse de los residuos, controlando la contaminacion del aire, aguas y tierras, 
aplicando la reduccion, reutilizacion, recuperacion y el reciclaje de los mismos. 

4°. Buscar la calidad de vida. La calidad de vida de los meridenos se lograra buscando un 
equilibrio entre el hombre y la naturaleza, para alcanzar un estado de bienestar fisico, mental, 
cultural y social, valorado en una dimension posible y equitativa, salud, trabajo, educacion, 
cultura, vivienda, servicios, ambiente, tiempo libre, economia, interaccion social, seguridad 
personal y juridica. Se va creando un nuevo contexto cultural alimentado por realizaciones 
sociales. 

5°. Tener una ciudad compacta. Merida debe responder a un sistema urbano compacto, que 
promueva la diversidad cultural, biologica, social, estetica y paisajistica, con asentamientos 
controlados que convoquen a la proximidad, a la conectividad ciudadana, al uso del 
transporte publico, al ahorro de energia, al amor y respeto por la ciudad, al triunfo de la 
serenidad sobre la velocidad, para formar ciudadanos responsables. El proceso de gestion 
urbana sostenible debe ser eficiente, es fundamentalmente politico, pero no es un fin, en si 
mismo. 

6°. Propiciar la participacion. Estimular la participacion de todas las capas sociales, para 
generar y legitimar las decisiones de gobernabilidad local, favorecer su implementacion, tener 
mejores mecanismos de comunicacion e informacion, aprender lecciones de ciudadania, 
hacerse mas sociables, definir una identidad individual y colectiva, precisar las prioridades, 
defender intereses y enriquecer la vida humana. Las comunidades deben prepararse para 
tomar con dignidad decisiones democraticas, aceptando voluntariamente limitaciones 
colectivas. 

7°. Fortalecer la localizacion. La ciudad debe revalorizar su entorno natural y estetico, 
vincular su desarrollo a una concepcion planetaria, dar respuestas al mundo globalizado, 
exaltar y ancestrar los valores locales del paisaje, la cultura, el arte, la historia, el agua limpia, 
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la fertilidad, la silvicultura, la productividad, el comercio, la fortaleza de las organizaciones 
locales, la competitividad y la participacion. Los barrios deberian convertirse en los 
“vecindarios urban os sostenibles". 

8°. Tener una ciudad segura. El desarrollo sostenible en los ecosistemas de montana es 
incierto, por los riesgos que causa la susceptibilidad a la erosion, los deslizamientos de tierra, 
la sensibilidad climatica, los movimientos sismicos y las amenazas socio-naturales, el deterioro 
ambiental, las precipitaciones, la sequia y la desertizacion. Hay que crear una cultura de 
prevencion y conservacion, encontrar soluciones seguras a los moradores en riesgo 
permanente, limpiar los cauces de rlos y quebradas, crear conciencia de la peligrosidad ciclica 
y hacer mas firmes los controles y la vigilancia ambiental. 

9°. Atacar la pobreza urbana. Hay que adelantar programas para aliviar la pobreza urbana, 
buscar con tenacidad la solucion al hambre, vivienda, malnutricion, calidad ambiental, 
salubridad, servicios, y educacion de pobres y marginales. Encontrar soluciones prioritarias a la 
feminizacion de la pobreza, los derechos de la mujer y del nino, la inequidad, la exclusion, la 
discriminacion, el sub-empleo, la inseguridad alimentaria y la violencia. Los proyectos de 
desarrollo sostenible no prosperaran sin la atencion de la pobreza presente. 

10°. Sensibilizar a los ejecutores. El desarrollo sostenible es un objetivo idealizado y 
complejo, y para alcanzarlo hay que basarse en la educacion civica y ambiental. Que sus 
ejecutores, los alcaldes del area metropolitana puedan abrir espacios de participacion, 
convocar a la ciudadania, vitalizar la sociabilidad, oir y convocar a personas e instituciones, 
aceptar la asesoria especializada, compartir los propositos y las decisiones, integrar 
voluntades, garantizar los servicios eficientes, la seguridad ciudadana y conciliar dentro de sus 
territorios de gestion, asumiendo cambios de mentalidad, nuevos conocimientos e 
intercambios, pensar mas en la gente, y buscar que Merida sea una ciudad ejemplar. 

Para seguir proyectos de ciudades sostenibles debemos preguntarnos ^Como se han 
comenzado a lograr ciudades sostenibles?. Las experiencias dicen que los metodos son 
distintos, pero en todos ellos ha privado la participacion, la educacion y la mistica ciudadana. 

*. Resolviendo un problema de transporte y transito. 

*. Con un programa efectivo de disposicion de residuos. 

*. Logrando mejorar la calidad de los servicios publicos. 

*. Desarrollo de un programa de seguridad y vulnerabilidad. 

*. Haciendo un plan de ordenamiento territorial. 

*. Adoptando a los rios como proyectos socio-ambientales. 

*. Adelantando programas de fortalecimiento social de los barrios. 

7. CONCLUSIONES. 

*. Respecto al Desempeno Estructural 

Considerar las amenazas multiples, la desagregacion de las amenazas presentes, dar mejor 
atencion a las amenazas de baja probabilidad de excedencia por su caracter caotico, 
establecer una normalizacion basada en desempeno, promover la evaluacion estructural de las 
edificaciones esenciales, considerar la relacion entre los sistemas estructurales y no 
estructurales, aplicar en las rutinas de calculo metodos basados en desempeno estructural y 
con mas frecuencia curvas de fragilidad, aplicar la salud estructural y estar pendiente de la 
evolucion de los nuevos metodos de calculo, de construccion y de mantenimiento de los 
sistemas estructurales. 

*. Respecto al Riesgo. 

Involucrarse en la gestion integral de los riesgos, intervenir en la planificacion urbana segura, 
entender las relaciones entre sociedad y naturaleza para detectar la presencia de los riesgos 
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permanentes fisicos y sociales. El ingeniero debe manejar los programas (GIS) del sistema de 
informacion geografica. 

*. Respecto a la Resilencia. 

Reconocer las potencialidades de la Resilencia fisica y social como alternativa integral y util 
para preparar a la poblacion, mitigar la pobreza y reducir las probabilidades de falla estructural, 
las grandes perdidas economicas y sociales, teniendo en cuenta que el analisis es 
probabilistico y multi-dimensional. 

*. Respecto a la Sostenibilidad. 

Manejar los conceptos de la sostenibilidad urbana y rural para responder a las 
responsabilidades profesionales y a la satisfaccion de la calidad de vida de los seres humanos. 
Un proyecto de ciudad debe contar con la cooperacion positiva de las instituciones 
universitarias, de su estructura academica y de sus personeros, los profesores con sus 
conocimientos, capacidad y experiencia, y los estudiantes, vigor de juventud y una oportunidad 
de formacion. La integracion de todo el subsistema de educacion superior presente, fortalecera 
la pertinencia en la definicion de nuevas obligaciones y compromisos con la sociedad, en las 
perspectivas de lograr el desarrollo sostenible, ambiental, social, economico y cultural, como 
una planificacion concertada entre las alcaldias, los especialistas y las comunidades, dando 
inmensas oportunidades para la cultura, la ciencia y la tecnologia. 

*. Respecto al total: 

Tener una concepcion holfstica para integrar los conceptos de desempeno, riesgo, resilencia y 
sostenibilidad para aplicarlos integralmente como ingenieria social y reconocer que este 
ejercicio obliga a ampliar los espacios de la especializacion sismorresistente. 
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RESUMEN 



A raiz de los trabajos de Lagrange en pandeo, se presenta primero una simplificacion 
de la situacion pos critica de una columna. Se simplifica la formula final, matematicamente 
impecable, para adaptarla a las situaciones pos criticas practicas. A continuacion, basandose 
de nuevo en los trabajos de Lagrange, geniales en sus planteos, no asi en sus resultados, 
seguidos por uno de Keller, impecable en su planteo matematico y en sus resultados, se 
presenta una pequena generalizacion de la busqueda de la mejor columna para pandeo, o sea 
de como debe variar la inercia a lo largo de la columna. Despues se analiza la mejor forma que 
debe tener la seccion de la columna. 

ON LAGRANGE’S WORK ON THE BUCKLING OF COLUMNS 

Works of Lagrange on the buckling of columns are considered in this paper. First, a 
simplification of the postcritical conditions in a column is discussed. The final equation due to 
Lagrange (which is a fine piece of mathematics) is simplified for its use in practical situations. 
Based on the high quality works of Lagrange, which on the other hand had problems in the 
results, together with the mathematically impeccable work by Keller, a new generalization is 
presented in the search for the best column to resist buckling; i.e. identification of inertia 
variations along the column to improve buckling strength. Finally, the problem of best column 
shape is addressed. 



1 INTRODUCCION 

Entusiasmado por unos trabajos de Lagrange y de Keller que pasaron por mis manos, 
tanto sobre el aspecto pos critico en pandeo como sobre la busqueda de la mejor forma de una 
columna, me parecio que en el primer caso la formulacion final necesitaba una simplificacion 
para adaptarla a las situaciones practicas y en el segundo caso una posible generalizacion. 

Enamorado por la genialidad de los trabajos de estos ilustres matematicos, este trabajo 
quiere ser una modesta contribucion sobre el tema de pandeo. 



2. LA TRAYECTORIA POS CRITICA SEGUN LAGRANGE 

Es bien conocida la formula de Euler para pandeo de una columna prismatica 
biarticulada: 



tt 2 EI 

i 1 



( 1 ) 
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Donde P cr es la carga critica, E el modulo de elasticidad, I el menor momento de inercia 
de la seccion de la columna, £ el largo de la columna. 

Lagrange (1770 -1773) reconoce la validez de dicha formula y parece ser el primero en 
preocuparse por el hecho que la indeterminacion de la flecha cuando llegamos a la carga 
critica no satisface una mente despierta, o dicho de otra manera, no convence del todo esta 
situacion de bifurcacion del estado de equilibrio. Lagrange parece asi ser el primero en analizar 
el fenomeno de pandeo para grandes deformaciones, con la ecuacion diferencial mas 
elaborada de la Ifnea elastica. Esta ecuacion no es sencilla y el tiene el merito de haber logrado 
su integracion y dar asi un aspecto de la situacion pos-critica por primera vez. 

En el diagrama P-f, la asintota horizontal de Euler en P cr se cambia asi en una leve 
curva, y la ecuacion de la misma fue proporcionada por Lagrange, con un habil desarrollo en 
serie. 



Si llamamos P cr la carga critica de Euler, P la carga pos-critica, f la flecha transversal 
de la columna, £ la luz entre apoyos articulados, la expresion de Lagrange, adaptada a 
nuestras notaciones, sera: 




= 1 + 
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— = 1,0561 ; Si — = 0,3 
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Notamos asi que dicha expresion no solamente vale para una carga P superior a la 

/ 

carga critica de Euler P cr , sino que para valor ya importantes de la relacion — , P es muy poco 
superior a P cr . 



En realidad desde un punto de vista practico, no se pueden aceptar relaciones de 
superior a 0,1 . 



/ 

£ 



Para graficar la funcion de Lagrange, vamos a considerar primero un solo termino del 
desarrollo en serie, ademas del termino del principio, con lo que f-i resulta de la siguiente 
manera: 



/i 



= M. . p 

n P P cr 



A continuacion consideramos un termino mas, obteniendo el valor f 2 siguiente. 



h 



4i P rr 

3.x' P 




- 8-1 



(3) 



( 4 ) 



De un libro de C. Massonnet (1960), obtenemos la formula para f vm atribuida a Von 
Mises, que es la siguiente. 
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En realidad, si realmente tenemos en cuenta que P es solamente muy ligeramente 
superior a P cr , todas estas formulas coinciden en una sola mas simple que proponemos en su 
reemplazo: 



/ 



4i 

K 




( 6 ) 



La obtencion de dicha formula se obtiene de la siguiente manera: Se considera que 
como P es muy poco superior a P cr , P/P cr puede considerarse como 1 + e, con e muy pequeno. 
Conociendo los habituales desarrollos: 

J g 

-JlVs &l-\- s 1 2 ; « !-<£■ ; y reemplazando en todas 

l + £ 1 + £■ 

las expresiones anteriores (3 - 4 - 5) se obtiene que cada una se simplifica en (6) 

Si graficamos estas cuatro funciones, en la Figura 1, vemos que coinciden 
perfectamente en el dominio de uso. 

Queda bien claro no obstante que dicho estudio realizado por Lagrange, no modifica el 
valor de la carga crftica de Euler y por lo tanto en todos sus calculos subsiguientes, Lagrange 
usa la ecuacion diferencial de Euler por considerarla suficiente. 



3. FORMA OPTIMA DE UNA COLUMNA, SEGUN LAGRANGE 

En un trabajo de su juventud, Lagrange (1770 - 1773) aborda la busqueda de la mejor 
forma de columna, entendiendo por ello la forma que conduce a la mayor carga crftica. Si bien 
sus resultados no son correctos, su metodo de trabajo es realmente genial y da pie, aun 200 
anos despues, a varios estudios sobre el tema. Su idea fundamental es elegir una variacion de 
inercia en forma de seccion conica y dar la forma general para integrar la ecuacion diferencial 
de Euler en este caso. Truesdell (1960) retoma el mismo trabajo, con ligeras correcciones. 



Evidentemente, si comenzamos con una columna prismatica tomada de referenda 
(Inercia = l 0 ) y modificamos la inercia para mejorar la carga crftica, debemos primero considerar 
que variamos el tamano de la seccion sin cambiar su forma, y en segundo lugar debemos 
mantener una referenda constante (l 0 ) para comparar. Para cumplir este segundo aspecto, una 
posibilidad es conservar el volumen de la columna. 



rl 1 r( 

Q n .t = Q ax o sea Q n = — \ Q dx 

u Jo 0 £ Jo 



(7) 



Como no cambiamos el largo de la columna, elegimos como referenda la columna con 
una seccion Q 0 que sea promedio de las secciones. 



En realidad, como el fenomeno de carga crftica esta intrfnsecamente ligado a la Inercia 
(minima) de la seccion de la columna, se podrfa proponer como referenda un promedio de la 
inercia: 





90 



Gerald Pirard 



Por otro lado, si se restringe el analisis a una columna circular, la variacion de inercia 
esta ligada directamente a la variacion del radio, y en este caso surge naturalmente de 
considerar como referencia la columna circular con un radio promedio. 



r( 

rdx 

Jo 



( 9 ) 



Vemos con estas consideraciones que existen varias posibilidades de referencia. 
Vamos en este trabajo a considerar cuatro casos de referencia. 

Teniendo en cuenta la primera consideracion que hicimos, cuando una seccion cambia 
de tamano sin cambiar de forma, utilizando el indice 0 para la seccion de referencia, tenemos 
por ejemplo que el promedio de seccion puede escribirse: 

a =\i iidx ( io ) 

y el promedio del radio en la columna circular seria: 

(/„)"" =fj on 

Generalizando asi, consideramos cuatro casos con k = 4, 3, 2, 1 y la referencia sera: 

z d 1 ft / d 

( h ) 4 =- J 0 ( 1 ) 4 dx ( 12 ) 

Notamos asi que el caso k = 4 es el promedio de inercia, k = 2 es promedio de la seccion, k 
= 1 seria el promedio del radio si la seccion fuera circular. El caso k = 3 no tiene sentido 
geometrico definido pero entra en la generalizacion. 



Estudiamos asi diferentes variaciones de inercia. En primer termino consideramos una 
variacion de inercia en seccion conica como propuesta por Lagrange, obteniendo asi los 
resultados de la primera columna de la Tabla 1 adjunta. 

Despues se estudio una columna con el tramo central de inercia constante y los extremos 
con variacion conica, obteniendose los resultados de la segunda columna de la Tabla 1. 

En una tercera etapa estudiamos una columna con dos tramos de inercia constante, 
resultados en la tercera columna de la Tabla 1. 

Finalmente damos en la ultima columna, para los diferentes casos, una adaptacion de la 
solucion de Keller al mismo problema. 

Se estudiaron tambien variaciones en forma tronco-conica, resultando que la mejor forma de 
la columna se obtiene cuando se transforma en prismatica. 

El estudio de los resultados muestra evidentemente que la maxima carga critica (S m a X ) 
corresponde a la solucion de Keller (1960), lo que no es de extranar ya que el planteo de 
dicha solucion preconiza que sera la maxima. No obstante la solucion conica se acerca 
mucho a la solucion de Keller y en cierta forma la aventaja ya que en los extremos superior 
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e inferior, Keller obtiene seccion e inercia nula, lo que de manera practica complica la 
colocacion de la carga, problema que no presenta la solucion en seccion conica. 

En la Tabla 1 podemos ver una representacion esquematica de la variacion de inercia a 
lo largo de la columna. 

Tambien es interesante introducir la nocion de eficiencia de la columna, de manera 
algo diferente de autores anteriores. Consideramos la eficiencia (Ef) de una columna como 
el cociente entre la carga critica y su volumen. 

En nuestro caso, como la carga critica de referenda se considera constante, 
podemos escribir, modificando el cociente senalado por factores constantes: 



Ef = o bien yy 



1 cr 0 

V 



V 0 _ s n Q i 

'V 

a Jc 



S 

i Q 

°^ n 



dx 



p Jc 



(13) 



e Jo /„ 



dx 



Evidentemente en el caso k = 2, la eficiencia se reduce a S, ya que en este caso 
tenemos conservacion del volumen. 

En la Tabla 2 se dan los valores de la eficiencia asi definida y tambien la relacion de 
volumen. Como los graficos de la primera tabla dan las variaciones de inercia, creemos 
conveniente representar en esta segunda tabla las mismas variaciones pero referidas a 
radios de columnas circulares equivalentes: 



r = 




y =4 



n 



(14) 



Notamos que las variaciones de volumen si bien existen no son muy grandes. La 
mayor eficiencia la Neva siempre el Caso 1 y la solucion de Keller (1960) evidentemente. La 
columna con menos eficiencia es la con tramos constantes. Las otras tres son muy 
parecidas en sus resultados. La columna con tramo central constante y extremos conicos es 
probablemente la mas simple a realizar practicamente. 



4. LA MEJOR FORMA DE UNA COLUMNA PARA RESISTIR CARGA AXIAL 



De los resultados anteriores notamos que para tratar de mejorar la carga critica de 
una columna, aumentamos la inercia en el tramo central y la disminuimos en los extremos 
por razones geometricas segun los cuatro casos analizados. 

Si consideramos en la Fig. 2, uno de los casos, notamos que desde un punto de 
vista practico, lo que hemos hecho aumentando la inercia al medio y disminuyendola en los 
extremos, es casi como indicado por la linea de puntos, es decir aumentar la inercia en el 
medio y disminuir la longitud. 

De hecho consideramos la expresion siguiente de la carga critica: 



n 2 El 

t 2 



(15) 



Y con variaciones diferenciales tenemos: 



dP cr _dl 2 

-y 



(16) 





92 



Gerald Pirard 



Notamos asf que como previsto un aumento de I, aumenta la carga crftica y una 
disminucion de i tambien. 



Segun los diferentes casos analizados anteriormente, tendremos: 
sea que reemplazando: 



M 

T 



k dl 
- — . — o 
4 I 



dP cr dl k dl k + 2 dl 



(17) 



Es decir, segun el caso tenemos la segunda columna de la Tabla 3. 

Por ejemplo si queremos un aumento de carga crftica del orden del 20%, tendremos 
la variacion relativa de inercia de la tercera columna de la Tabla 3, y la correlativa 
disminucion de largo de la cuarta columna. 

Notamos que la disminucion de longitud es muy pequena. En realidad vemos que 
buscar aumentar la carga crftica significa de hecho aumentar la inercia en el tramo central 
con una disminucion casi imperceptible de la longitud. 



Tabla 1 de resultados 
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Tabla 2 de Volumen y eficiencia: . 



P cr o V 



Conica 



Constante y conica 



Dos tramos 
constantes 



Keller 



Radio 
de la 
columna 
circular 
equivalente. 

4 : 
r = 4 .1 

K 








V/V 0 


Eficiencia 


V/V 0 


Eficiencia 


V/Vo 


Eficiencia 


V/Vo 


Eficiencia 


Caso 4 


0,965 


1,248 


0,947 


1,255 


0,957 


1,183 


0,962 


1,264 


Caso 3 


0,981 


1,279 


0,970 


1,286 


0,977 


1,193 


0,975 


1,294 


Caso 2 


1,000 


1,317 


1,000 


1,326 


1,000 


1,207 


1,000 


1,333 


Caso 1 


1,022 


1,365 


1,040 


1,379 


1,024 


1,229 


1,050 


1,388 



Tabla 3 




cr 

Per 


— % 

/ 


— % 

i 


Observacion 


3 d - 1 - 

I 


6,66 


-6,66 


Corresponde a disminucion de 3,33% 
en cada extremo 


5 dl 
2' I 


8,00 


-6,00 


Corresponde a disminucion de 3,00% 
en cada extremo 


2.— 

I 


10,00 


-5,00 


Corresponde a disminucion de 2,50% 
en cada extremo 


3 dl 
2 ' 1 


13,33 


-3,33 


Corresponde a disminucion de 1,66% 
en cada extremo 
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5. SOBRE LA FORMA DE LA SECCION 

En la seccion anterior se supuso que la inercia de la columna varia con una cierta 
ley, pero que la forma de la seccion se mantiene parecida a si misma y a la seccion de 
referencia. 

No hemos hecho ninguna restriccion en cuanto a la forma de la seccion, si bien 
queda claro que para el pandeo de la columna biarticulada, es la inercia minima la que 
debemos tener en cuenta. 



Ademas de los resultados anteriores que muestran como debe variar la inercia para 
aumentar la carga critica, podemos abocarnos ahora a ver cual es la mejor seccion o sea la 
que tenga mayor inercia minima. Por supuesto de nuevo debemos mantener algo fijo, sino 
podemos aumentar indefinidamente la inercia tanto maxima como minima. Elegimos de 
manera natural mantener fijo el valor del area de la seccion. 

Independiente de la forma de la seccion, el estudio de las inercias muestra que la unica 
manera de tener una inercia minima que sea maxima es cuando alcanza la inercia maxima, 
o sea cuando el circulo de Mohr de las inercias se reduce a un punto. 

Por ejemplo, todos los poligonos regulares entran en dicha categoria. La inercia de 
un poligono regular de n lados es, si llamamos “a” el radio circunscrito: 



_ n.a 2.71 

I„ = ——.sen( ). 

8 n 



7 \ 2 

7J7 

cos(— ) 

V n J 



r 



Y 



+ - 






Y la superficie del mismo es: 



Q 



n 



n.a 2 .2.71 

sein 

2 



(18) 



(19) 



Si calculamos la relacion entre la inercia de un poligono de n lados a la inercia de un 
circulo de misma area, tenemos: 

2.71 

q(n) = 

sen(—) 
n 

Graficando esta funcion obtenemos la Fig. 3, donde notamos que evidentemente 
hemos hecho variar n a partir de 3 y que para este valor que corresponde al triangulo 
equilatero tenemos la maxima relacion de 1,209. 

Elio significa que si usamos una columna en forma de triangulo equilatero en vez de 
una columna circular, ya tendremos casi un 21% de aumento de la carga critica, a mas del 
aumento que podemos tener por variar la inercia a lo largo de la columna. Dicha ganancia 
se reduce muy rapidamente a medida que aumentan los lados del poligono. Para un 
hexagono seria solamente de un 0,7%. 

En otras palabras podemos concluir tambien que desde el punto de vista de la carga 
critica, la columna circular es la menos indicada. 

Seguimos realizando consideraciones sobre la mejor forma sin variacion del area de 
la seccion. 

Sea una seccion cualquiera referida a sus ejes principales de inercia y, z, como se 
muestra en la Fig. 4. Una seccion asi puede siempre encerrarse en un rectangulo de 
dimension (a, b). 



( n V 

cos(— ) 

V n ) 



sen(—) 

V n J 



( 20 ) 
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El area Q 0 siempre puede expresarse como: 
Q 0 = a.b.a 



( 21 ) 



o sea por el area del rectangulo multiplicada por un coeficiente a de forma. Si modificamos 
la seccion proporcionalmente a si misma, a y b varian proporcionalmente al cambio, pero a 
se mantiene. 



En cuanto a las inercias podemos escribir: 



I 



zO 




_ b.a 
~ ~Vl‘‘ 



( 22 ) 



donde p y y son coeficientes de forma de la misma manera que a. 

Una forma de modificar la seccion sin cambiar su superficie Q 0 , es modificando 
todas las dimensiones en la direccion “y” en un factor k y en la direccion “z” en un factor 1/k. 
Asi “b” sera “b.k” y “a” sera “a/k”, de manera que la seccion no cambia: 



Q = — .b.k. a = a.b.a = Q n 
k 0 

En cuanto a las inercias: 



o sea: 



h = 



,(b.k) 3 

12 



■P 






b.k. 






I = 






12 



■r 



i = 



yo 



(23) 



(24) 

(25) 



Una inercia disminuye y la otra aumenta en la misma proporcion. Este procedimiento 
permite siempre igualar las inercias maxima y minima sin cambiar el area. Notaremos que el 
producto de las inercias no cambia. 



I t J y =I a JyO 



(26) 



Este procedimiento indicado permite siempre obtener la maxima inercia minima, que 

sera: 



^ Max I min ZO ' vO (27) 

De nuevo insistimos que este resultado no depende de la forma de la seccion. Este 
procedimiento aumenta la inercia minima con una correspondiente disminucion de la inercia 
maxima, hasta que se igualen. 

Una columna rectangular nunca podra tener mayor carga critica que la columna 
cuadrada de misma area. 

Podemos perfectamente tambien aumentar la inercia minima sin modificar el area 
de la seccion ni el valor de la inercia maxima y llegar asi a igualar las inercias principales. Si 
consideramos Fig. 5a, una seccion cualquiera referida a sus ejes principales “z” de maxima 
inercia e “y” de minima. Si consideramos un rectangulo elemental paralelo al eje “z”, Fig. 5b, 
podemos dividirlo en dos de igual superficie y deslizar cada una de estas de una misma 
cantidad, una a la izquierda y la otra a la derecha de una misma cantidad arbitraria. Elio no 
modifica ni el area ni la inercia relativa a “z” y aumenta la inercia relativa a “y”. Si cuidamos 
de mover estas superficies iguales de una misma distancia, los ejes “y, z” seguiran siendo 
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ejes principales de inercia. Si seguimos asi con los sucesivos rectangulos elementales 
superiores obtenemos la Fig. 5c. 

De esta forma podemos aumentar la inercia minima hasta igualarla a la maxima. 

Si analizamos una columna con seccion rectangular, primero consideramos que 
base y altura del rectangulo varian proporcionalmente. Ese caso entra directamente en los 
resultados de las Tablas 1 y 2. 

Si se mantiene la base constante, variando solamente la altura, no debemos olvidar 
que el pandeo se rige por la inercia minima. Por otro lado la hipotesis que hicimos en los 
calculos anteriores de que las secciones varian manteniendo su forma ya no tiene validez. 
En este caso puede estudiarse el problema por el metodo de Lagrange y obtener el mismo 
resultado que el caso 4 de la Tabla 1 con una ganancia sobre la carga critica de 20,7%. Se 
debera tener cuidado que la altura en los extremos sea mayor o a lo sumo igual a la base 
del rectangulo. El metodo de Keller (1960) no puede utilizarse en este caso, ya que en este 
metodo la inercia varia de un maximo al centra hasta cero en los extremos. Elio obliga a que 
en algun punto intermedio, la altura llegue a igualarse a la base (cuadrado) y mas alia de 
este punto la direccion de inercia minima gira de 90°, Fig. 6. 

6 CONCLUSIONES 

• Hemos consultados varios autores ilustres, y si bien muchas veces hemos 
introducido modificaciones y/o correcciones a sus desarrollos, nuestra 
meta no fue criticar sus trabajos sino agradecerles el trabajo que han 
realizado y el surco que emprendieron, tratando de prolongarlo y 
profundizarlo. 

• En cierta forma es un trabajo de recopilacion, sintesis y algo de 
generalizacion. 

• Hemos eliminado todos los pasos matematicos, algunas veces pesados 
para presentar solamente los resultados. Dichos calculos, si bien cambian 
un poco con respecto a los autores originales, no tienen nada de 
fundamentalmente nuevo. Las novedades que puedan existir estan en 
estas paginas en cuanto a sus resultados. 

• En primer lugar introducimos una simplificacion en la expresion de la carga 
de una columna en estado pos critico. Es un simple paso matematico, 
considerando que P solamente puede sobrepasar muy poco P cr . 

• Generalizamos la busqueda de la mejor variacion de inercia de una 
columna, con diferentes referencias. Evidentemente la solucion propuesta 
por Keller (1960) es siempre la mejor forma ya que su planteo analitico lo 
preconiza. No obstante la forma propuesta por Lagrange (1770 - 1773) se 
acerca mucho a la mejor forma, con la ventaja de no tener seccion nula en 
los extremos. 

• Del punto de vista de la eficiencia, como la hemos introducido, que cambia 
ligeramente con la definicion de autores anteriores, sigue siendo la forma 
propuesta por Keller (1960) que tiene la maxima eficiencia. El caso con 
tramo central constante y extremos con variacion de inercia conica es el 
que sigue, de muy cerca, con la ventaja de nuevo de no tener seccion nula 
en los extremos. 

• En cuanto a la mejor forma de la seccion, no hay duda que la igualdad de 
las inercias principales es el requisito mas relevante. Elio implica que 
ninguna columna con seccion rectangular pueda tener mayor carga critica 
que la viga cuadrada de misma area. De las secciones en forma de 
poligono regular, el triangulo equilatero proporciona la mejor forma... y el 
circulo la peor. 
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RESUMEN 

El uso de los materiales compuestos para la rehabilitacion y la reparacion de 
estructuras de hormigon y de acero se ha ido incrementando en las ultimas decadas. Frente a 
la necesidad de entender de una forma mas apropiada el comportamiento de las estructuras 
reforzadas, la simulacion numerica de los materiales compuestos ha surgido como una 
herramienta; sin embargo, es un campo que aun esta en desarrollo. En este trabajo se 
desarrolla la teoria de mezclas serie/paralelo como una alternativa eficiente para el analisis de 
estructuras reforzadas con polfmeros reforzados con fibras largas (FRP) y se estudia con 
detalle el comportamiento de la interfaz hormigon-adhesivo-FRP. Utilizando el programa de 
elementos finitos no lineal PLCDYN [iii], se hace una simulacion tridimensional de un ensayo 
de doble cortante y se compara con los resultados obtenidos en el laboratorio. Se han podido 
evidenciar las ventajas de esta teoria como herramienta para el analisis numerico de los 
materiales compuestos. 



ABSTRACT 

The use of composite materials in the rehabilitation and reparation of concrete and steel 
structures has increased in the last decades. However, the behaviour of the structures 
strengthened with composite materials is complex and still under study. Therefore, the 
numerical simulations of such structures appeared to be a necessary as tool, although it is still a 
developing field. This paper develops the serial/parallel mixing theory as an efficient alternative 
for the analysis and design of structures reinforced or rehabilitated with Fiber Reinforced 
Polymers (FRP) and studies in detail the behaviour of the FRP- concrete interface. The 
simulation of a double shear laboratory test was performed by using the non linear finite 
element code PLCDYN [iii], and the obtained results are compared whit experimental ones. The 
used theory and computer code proved to be excellent computational tools for the numerical 
analysis of composite materials. 

1. INTRODUCTION 

Frente al incremento del uso de los materiales compuestos conformados por polfmeros 
reforzados con fibras largas (FRP) como una alternativa de refuerzo de las estructuras, es 
importante buscar nuevas herramientas que permitan estudiar su comportamiento para 
optimizar su utilizacion diseno tanto en la reparacion de los elementos estructurales con cierto 
nivel de dano como en la rehabilitacion de los elementos que requieren tener un mejor 
comportamiento dentro de la estructura. 

Los polfmeros reforzados con fibras largas (FRP) mas utilizados como refuerzo son los 
reforzados con fibra de vidrio (GFRP) y con fibra de carbono (CFRP); en menor medida han 
sido utilizados los polfmeros reforzados con fibra de aramida (AFRP). 
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La interaccion entre la experimentacion y la simulacion numerica con elementos finitos 
representan una estrategia eficiente para estudiar el comportamiento de las estructuras 
reforzadas con materiales compuestos. Dado que la simulacion numerica de los materiales 
compuestos es un campo que esta en desarrollo, y que su uso en obras civiles es 
relativamente nuevo, aun hay aspectos de los FRP que, por su complejidad, no han sido 
suficientemente estudiados para su utilizacion como refuerzo en estructuras de hormigon 
armado. Entre estos aspectos hay que destacar, la perdida de adherencia entre el FRP y el 
hormigon. 

Entre los resultados experimentales existentes se ha encontrado que los mecanismos de 
fallo de la adherencia del FRP son factores criticos para el comportamiento de los elementos 
reforzados solicitados a cortante o flexion. Ademas, dejan claro que existe la necesidad de 
nuevas herramientas que permitan predecir, de una forma mas precisa la perdida de 
adherencia por el deslizamiento entre el FRP y el hormigon o por la delaminacion del FRP [i], 

En este artfculo, se comienza por hacer una revision del estado del conocimiento de las 
investigaciones realizadas en el campo del comportamiento de la interfaz hormigon-adhesivo- 
FRP, analizando la influencia del comportamiento de los materiales y su interaccion con la 
perdida de adherencia de la interfaz. Luego se describe la teorfa de mezclas serie/paralelo [ii], 
la cual es una nueva herramienta que permite la simulacion numerica de los materiales 
compuesto. Por ultimo, utilizando el programa de elementos finitos PLCDYN [Hi], se Simula un 
ensayo de doble cortante utilizando la teorfa de mezclas serie/paralelo, en el que se analiza el 
comportamiento de la interfaz hormigon-epoxi-FRP; por ultimo, se contrastan los resultados 
numericos obtenidos con los resultados experimentales. 

2. ESTADO DEL CONOCIMIENTO SOBRE LA SIMULACION NUMERICA DE LA 
INTERFAZ HORMIGON-EPOXI -REFUERZO CON FRP 

La eficiencia del refuerzo para las estructuras de hormigon armado depende tanto de la 
apropiada seleccion de la configuracion del refuerzo, como de la adherencia entre el FRP y el 
compuesto [iv], 

Los estudios experimentales realizados demuestran que uno de los parametros criticos para 
definir la resistencia del hormigon armado reforzado con FRP es la adherencia entre el 
hormigon y el FRP [v], de la que depende el mecanismo de transmision de tensiones y por 
tanto del comportamiento del compuesto. Consecuentemente, el comportamiento de la interfaz 
hormigon-epoxi-laminado en FRP es uno de los aspectos que ha sido mas estudiado; la 
mayorfa de los estudios han sido enfocados en la parte experimental, para determinar la 
resistencia de la adherencia; sin embargo en muy pocas estudios se considera la simulacion 
numerica como una herramienta para mejorar el analisis del comportamiento de la interfaz. 
Frente a la complejidad del comportamiento de la interfaz hormigon-epoxi-FRP es evidente 
que, ademas de la experimentacion, se debe utilizar la simulacion numerica para evaluar el 
comportamiento estructural de los elementos reforzados, para optimizar los disenos y para 
garantizar la eficiencia del FRP como refuerzo [vi], 

En la decada pasada, gran parte de las simulaciones se realizaron para la simulacion de 
vigas reforzadas y alguna simularon los modos de fallo de la interfaz por deslizamiento, lo que 
por el grado de complejidad del problema, ha implicado un gran coste computacional [vii; viiij. 
En los primeros estudios, se realizaron excesivas simplificaciones; por ejemplo no se considero 
el dano en el hormigon y se supuso la adherencia perfecta entre el hormigon y el FRP durante 
todo el proceso de carga; ademas, en algunos casos no se modelo la capa del adhesivo o se 
considero como un material lineal y elastico. Por tanto, la precision de los resultados obtenidos 
mediante estas simulaciones es limitada, pero dieron las bases para establecer la distribucion 
de las tensiones normales y tangenciales en la zona interfacial hormigon-FRP. Posteriormente, 
los estudios han puesto entasis en el analisis no lineal con elementos finitos de la interfaz, con 
el objetivo de simular el proceso de deslizamiento entre el FRP y el hormigon [ixj. 
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El avance del uso de los elementos finitos en el analisis del comportamiento de la interfaz 
hormigon-epoxi-FRP ha sido relacionado con la simulacion de la propagacion de las fisuras en 
el hormigon y en la zona interfacial. En ese sentido se desarrollaron dos tipos de aproximacion 
es la simulacion de fisuras: 

• Modelos de fisura discreta. En esta formulacion se predeterminan la direccion y posicion 
de las fisuras en el mallado donde, por cada incremento de carga, se requiere que la malla 
de elementos finitos se modifique de acuerdo con la propagacion de cada una de las 
fisuras. Asegura una buena precision de la simulacion del proceso de formacion y 
propagacion de fisuras, pero tiene un elevado coste computacional por el remallado. Este 
modelo ha sido utilizado, entre otros, por Perera et al. [vii], Yang et al. [viii], Niu y Wu [x; xi] 
y Niu et al. [xii] quienes, a excepcion Perera et al., han enfocado sus modelos hacia el 
comportamiento de la interfaz en vigas reforzadas con FRP. 



• Modelos de fisura distribuida. La malla de elementos finitos no se modifica y el proceso 
de fisuracion se determina a traves de las ecuaciones constitutivas no lineales. Las fisuras 
son simuladas como discontinuidades locales que son distribuidas dentro del solido 
aplicando la teorfa de la mecanica del continuo. Entre los estudios que han utilizado este 
tipo de modelos para el analisis de la interfaz hormigon-epoxi-FRP estan Wu y Yin [xiii], Lu 
et al. [ix; xiv], Ebead y Neale [xv], Coronado y Lopez [xvi], 

Algunos investigadores han estudiado el deslizamiento entre el FRP y el hormigon 
considerando que el fallo ocurre por el modo de fractura de deslizamiento en el hormigon, 
debido a que en la mayorfa de los ensayos de adherencia, la superficie del hormigon fisura por 
cortante. Es por ello que simulan la interfaz tomando elementos de espesor nulo o simplemente 
no modelan el adhesivo [viii; x; xi; xii], En otros casos, colocan una capa de elementos de 
interfaz entre el FRP y el hormigon, en la que el deslizamiento es simulado en funcion de la ley 
constitutiva de la zona interfacial como un fallo de estos elementos [xiii; xv], Otra alternativa en 
la simulacion numerica es modelizar el deslizamiento mediante un modelo de fractura en la 
interfaz del hormigon, en los elementos adyacentes a la capa del adhesivo, donde el proceso 
del deslizamiento depende del modelo constitutivo asignado al hormigon [ix; xvi], 

A pesar que, desde el ano 2002, Ueda et al. [xvii] Recomendaron estudiar con mayor 
detalle el deterioro de las propiedades de la interfaz de las zonas cercanas a las fisuras en el 
hormigon, en las simulaciones recientes solo se ha considerado el dano en el hormigon y en 
pocos estudios se ha modelizado el adhesivo. Ademas, tampoco se ha tenido en cuenta que el 
epoxi puede tener un comportamiento no lineal en respuesta al dano del hormigon o presentar 
dano. 

Es importante que el FRP y el hormigon armado se puedan modelizar considerando el 
estado de cada no de los materiales que lo componen; una potente herramienta para ello es la 
teorfa de mezclas de serie/paralelo [xviii]. Asimismo, es necesario que en las simulaciones se 
considere el dano, tanto en el hormigon como en el adhesivo y en la matriz polimerica del FRP, 
bajo el principio de que el dano en los materiales produce la reduccion de las propiedades 
mecanicas de la interfaz, lo cual afecta la adherencia. 

3. TEORIA DE MEZCLAS SERIE/PARALELO PARA MATERIALES COMPUESTOS 

La teorfa de mezclas, como gestor del uso de los modelos constitutivos que describen el 
comportamiento de los materiales simples que componen los materiales compuestos, 
constituye una herramienta que puede ser utilizada dentro de una tecnica de elementos finitos 
para simular apropiadamente el comportamiento de dichos materiales. No obstante, es 
necesario el desarrollo de nuevas herramientas que mejoren la simulacion numerica a nivel 
global y local, de manera que se represente de una forma mas precisa el comportamiento real 
de los compuestos, con la finalidad de obtener informacion que no se puede medir en los 
ensayos de los elementos reforzados con FRP. Una excelente alternativa para la simulacion 
numerica es la teorfa de mezclas serie/paralelo [ii; xix], recientemente validada a traves de 
diferentes pruebas experimentales. Esta teorfa se fundamenta en la teorfa de mezclas clasica, 
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propuesta por Trusdell y Topin (1960), la cual considera que el comportamiento de cada 
componente define el comportamiento global del compuesto. 

La teorta de mezclas clasica surgio como una solucion a la limitacion que inicialmente 
presentaban los modelos constitutivos en el estudio de los materiales compuestos, por tratarlos 
desde un punto de vista macroscopico como un material homogeneo, sin tener en cuenta la 
contribucion de sus materiales componentes [xx], Esta teorta utiliza en la mecanica del 
continuo bajo el principio de interaccion de las sustancias que componen el material, 
suponiendo que en el volumen infinitesimal del compuesto participan en paralelo todos sus 
componentes; asimismo, se considera que cada uno de ellos contribuye al comportamiento del 
compuesto en la misma proporcion que su participacion volumetrica [xxi]. La teorta de mezclas 
clasica presenta una restriccion para su uso en los diferentes tipos de materiales compuestos, 
ya que parte de la hipotesis de que los componentes que coexisten en un punto, estan en 
paralelo y tienen la misma deformacion; por lo tanto, se puede utilizar en el campo lineal y, con 
algunas adaptaciones, en casos no lineales espectficos, siempre que no se consideren 
significativos los movimientos relativos entre los componentes del compuesto [xxii]. 

Debido a las diferencias en el comportamiento de los componentes, a su interaccion y a la 
forma en que estan distribuidos en el compuesto, los materiales compuestos tienen una 
elevada anisotropta y un comportamiento no lineal incluso para bajas tensiones, lo que hace 
que su modelacion numerica sea compleja. No obstante, a traves de la teorta de mezclas 
serie/paralelo complementada con el modelo de mapeo de espacios [xxiii] y a su 
generalizacion para el campo de grandes deformaciones [xxii], se obtiene la no linealidad y la 
anisotropta del compuesto como resultado de la interaccion de cada uno de los componentes 
con su correspondiente ecuacion constitutiva. 

3.1.1 Fundamentos de la teorta serie/paralelo [ii] 

En la teorta serie/paralelo se considera que en una direccion particular los componentes se 
comportan en paralelo y en las otras direcciones en serie. Se utilizan las siguientes hipotesis: 



• En cada volumen infinitesimal del compuesto participan en conjunto todos sus 
componentes, es decir, la distribucion de los componentes es homogenea. 

• Los componentes tienen la misma deformacion en la direccion en paralelo (condicion 
de isodeformacion) y la misma tension en la direccion en serie (condicion de 
isotension). 

• La adherencia entre los componentes es perfecta. 

• La contribucion de los componentes en la respuesta del compuesto es proporcional a la 
participacion volumetrica de cada material componente. 

• El volumen ocupado de cada material componente es inferior al volumen total del 
compuesto. 



3.1.2 Definicion de los componentes en serie y en paralelo de los tensores de 
deformacion y de tension [xxiv] 

El tensor deformacion e se descompone una parte en serie e s y otra en paralelo e p por 
medio de los tensores de proyeccion de cuarto orden (P P ,P S ) en paralelo y en serie, 
respectivamente 

£ — £p + G s (1 ) 

Sp — P^, : £ s s - P s : £ 

El tensor de tensiones a tambien se separa en sus componentes en serie a s y en 
paralelo <j p . 

CT = CJp + (7$ 

<7p — Pp : cj . = P s : o 



( 2 ) 
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donde los tensores de proyeccion de cuarto orden se hallan a traves del tensor de proyeccion 
paralelo de segundo orden N P y del tensor identidad de cuarto orden I 

P p = N p <8> N p P s = I - P p N p = e x ® <?, (3) 

siendo e* el vector unidad que indica la direccion del comportamiento en paralelo, es decir, la 
direccion de la fibra. 

3.1.3 Ecuaciones de equilibrio y de compatibilidad en las capas del compuesto [ xxv] 

Con el fin de minimizar el costo computacional, la implementacion numerica de la teoria fue 
desarrollada para descomponer el compuesto c en un determinado numero de capas neap, tal 
que cada capa j sea conformada por una matriz m y una fibra f. Consecuentemente, con las 
hipotesis enunciadas, en cada capa del compuesto j se plantean las siguientes ecuaciones de 
equilibrio y de compatibilidad de deformaciones: 

• Comportamiento en paralelo 



c _ m ^ f ^ 

j£ P - j&p- j£ P 


c mi m . f 7 f 

j a P= jk- j°r+j k -j°r 


( 4 ) 


Comportamiento en serie 






c m f 

j°s = J°s= j°S 


c c - "‘h . m p +f] c .f F 

jb s - jK ,-fc.s -T jK /«>.V 


( 5 ) 



3.1.4 Ecuacion constitutiva de los materiales de las capas del compuesto 



Como la implementacion de la teoria serie/paralelo se realiza a nivel constitutivo a partir del 
estado de deformacion en un punto de Gauss, una vez que se obtiene la deformacion del 
compuesto c e, se calcula el estado tenso-deformacional de cada componente cumpliendo con 
las ecuaciones de equilibrio y de compatibilidad, para luego hallar el estado de tensiones y de 
deformaciones del compuesto. Dicho estado se determina mediante la ecuacion constitutiva 
que rige el comportamiento de cada componente 



n n . 

7 a_ 



e 



n . 



)(e-e p ) 



' = { m - /} 



( 6 ) 



de tal forma que "C es el tensor constitutivo de la matriz m o de la fibra f de la capa j. 



3.1.5 Algoritmo de la teoria serie/paralelo 

En la figura 1 se muestra el algoritmo de esta teoria donde una vez que se determina el 
tensor deformacion del compuesto en cada capa j, se descompone el tensor de deformacion en 
sus partes en serie y en paralelo con el fin de calcular las correspondientes deformaciones en 
la matriz y en la fibra. Segun la teoria serie/paralelo, la deformacion en paralelo de los 
componentes es la misma, mientras que la deformacion en serie es diferente para cada 
componente. Por ello, en el analisis de la parte en serie se requiere una primera aproximacion 
de la deformacion en serie de alguno de los componentes, por lo que se toma como predictor el 
tensor deformacion en serie de la matriz del compuesto, para luego determinar el tensor de 
deformacion en serie de la fibra mediante la ecuacion (5) [xix]. 

La ecuacion (5) de equilibrio de tensiones se verifica calculando la tension residual de 
la capa y Aas y revisando que esta sea menor que una tolerancia Toler, tal como se muestra en 
la ecuacion (7). Cuando la tension residual es mayor, la prediccion de la deformacion de la 
matriz se corrige en forma iterativa. Una vez se obtiene una tension residual menor que la 
tolerancia, se hace la recomposicion de los tensores de tension y de deformacion de la capa j. 

.Aa s . = "ct s - f jG s < Toler (7) 
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3.1.6 Ecuacion de equilibrio en el compuesto 

Una vez realizado el analisis de cada capa, el tensor de tension obtenido para el compuesto 
c g corresponde a la suma de los tensores de tension de las capas del compuesto neap de 
acuerdo con el porcentaje de participacion volumetrica de cada capa k ■ 



M 7=1 



( 8 ) 



Deformaciones de la capa de compuesto j en un PG 
U C £ = C ,E 
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Figura 1 Esquema de solucion de la teoria de mezclas serie/paralelo en pequenas 
deformaciones para una capa de un compuesto en un punto de Gauss para un paso i+l [xviii] 



4.- ANALISIS DEL COMPORTAMIENTO DE LA ZONA INTERFACIAL HORMIGON- 
EPOXI- FRP 

En este apartado se comparan los resultados obtenidos utilizando la teoria serie/paralelo, 
con los resultados experimentales de un ensayo de doble cortante realizado en el Laboratorio 
de Resistencia de Materiales y Estructuras de la Escuela Tecnica Superior de Ingenieros 
Industrials de Terrassa de la Universidad Politecnica de Cataluna. 
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Son pocas las investigaciones que consideran los modelos tridimensionales, en el estudio 
del comportamiento de la interfaz hormigon-epoxi- FRP, pero en ellos realizan simplificaciones 
que no consideran el comportamiento no lineal la zona interfacial. Entre las mas recientes se 
encuentran: Un modelo tridimensional de un ensayo de doble cortante con placas de GFRP 
[xxvi]. El hormigon y el GFRP fueron representados con elementos solidos suponiendo un 
comportamiento elastico y lineal, mientras que la interface se simulo como dos sistemas de 
resortes de rigidez constante. Una simulacion tridimensional de un ensayo de code simple 
[xxvii]. Aunque se considero el comportamiento elastico lineal en el hormigon, el epoxi y el 
FRP, se evidencia que la distribucion de tensiones es diferente a la obtenida cuando se 
considera un estado piano de tensiones. Dicha diferencia que se debe tanto a la incidencia de 
la relacion entre el ancho del FRP y el del bloque de hormigon (b F Rp/b h ), como a la influencia 
los coeficientes de Poisson de cada uno de los materiales que conforma la interface. 

En este articulo se realiza una simulacion numerica del test de doble cortante en tres 
dimensiones, mediante el programa de elementos finitos PLCDYN [iii], y se analiza el 
comportamiento no lineal de la zona interfacial hormigon-epoxi- FRP, aprovechando las 
ventajas que ofrece la teoria serie/paralelo en determinar el comportamiento tanto de los 
componentes como del compuesto. 

4.1 Geometri'a del ensayo de doble cortante de la interfaz hormigon-epoxi-CFRP 

La geometria del ensayo se indica en la figura 2, donde las bandas de CFRP se encuentran 
adheridas desde el extremo del bloque de hormigon. Una de las laminas tiene menor longitud 
con el proposito de inducir el fallo en ese sector y que, de esta forma, se garantizara que la 
instrumentacion utilizada en el ensayo registrase todo el proceso de carga hasta el fallo. 

Vista lateral (Laminado largo) 

Zona con epoxi 




^ LOO ^ 

Figura 2 Geometria de la probeta del ensayo a doble cortante [xviii] 

Solo se modelizo el bloque al que se le aplica la fuerza de traccion, con una malla 
compuesto por hexaedros de 8 nodos. El modelo tiene cinco materiales simples y sus 
propiedades se muestran en la tabla 1 . 

De acuerdo con la teoria serie/ paralelo, en el hormigon armado y en el CFRP, las 
matrices de estos compuestos son el hormigon y la matriz polimerica, respectivamente, 
mientras que el acero y las fibras de carbono constituyen el refuerzo. Ademas, para considerar 
su correspondiente anisotropia, se tiene en cuenta que las barras de acero y las fibras de 
carbono estan orientadas a 0 grados de la direccion en que se aplica la carga. En la figura 3 se 
indican los cuatro materiales compuestos asignados al modelo y en la tabla 2 se muestran los 
porcentajes volumetricos en los que participan los componentes en cada material compuesto. 
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Tabla 1 Propiedades mecanicas de los materiales componentes [xviii] 



Material 


Criterio de 
limite 


Ex=Ey=Ez 

[MPa] 


V 


CTC 

[MPa] 


at 

[MPa] 


Gc 

[kPa-m] 


Gt 

[kPa-m] 


Hormigon 


Mohr-Coulomb 


3.65-1 0 4 


0.25 


52.0 


3.0 


55.0 


0.5 


Acero 


Von Mises 


2.0CM 0 b 


0.00 


270.0 


270.0 


2000.0 


2000.0 


Epoxi 


Mohr-Coulomb 


1 .25-1 0 4 


0.20 


80.2 


22.9 


36.0 


3.0 


Matriz polimerica 


Mohr-Coulomb 


1 .20-1 0 4 


0.20 


87.5 


29.2 


36.0 


3.0 


Fibra de carbono 


Von Mises 


2.26-1 0 b 


0.00 


2200.0 


2200.0 


14500.0 


14500.0 




Tabla 2 Porcentajes de los componentes en los materiales compuestos [xviii] 



Material 

compuesto 


Matriz de 
Hormigon 


Acero a 0 s 


Matriz de 
Resina 


Matriz 

Polimerica 


Fibra de 
Carbono a 
0 e 


1 


100 










2 


49.8 


50.2 








3 






100 






4 








0.32 


0.68 



4.2 Descripcion de los modelos constitutivos de los materiales 

En las recientes investigaciones de la interfaz hormigon-epoxi-FRP, en las que se han 
modelizado ensayos de adherencia, se han considerado diferentes criterios de dano para el 
hormigon. Sin embargo, se observa que no se ha tenido en cuenta el dano en el adhesivo o en 
la matriz del FRP, efectos que evidentemente, reducen las propiedades mecanicas de los 
materiales que componen la interfaz y, por ende, afectan la adherencia [xviii]. Por ello, en este 
articulo para analizar el comportamiento de la interfaz de una forma integral, se aplico el 
modelo de dano de Kachanov en el hormigon, el epoxi y la matriz polimerica del compuesto, 
mientras que para el acero se considero un modelo elasto-plastico y para la fibra de carbono se 
utilizo como un modelo elastico y lineal. En todos los casos se ha considerado la hipotesis, que 
no hay dano inicial en los materiales. Ademas se supuso un estado inicial de adherencia 
perfecta en la interfaces hormigon-epoxi, epoxi-CFRP y hormigon-acero. 

4.3 Rigidez y resistencia 

En la figura 4 se observa la curva de carga - desplazamiento de la simulacion numerica 
(SIM). En la prueba experimental (EXP), la carga maxima obtenida fue de 59.30kN y, en la 
simulacion, fue de 59.26kN. Por otra parte, tanto en el ensayo como en el modelo numerico, el 
fallo se produjo con poca ductilidad. 
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Sim: simulacion numerica exp: experimental 
Figura 4 P-A con el ensayo de doble cortante 



4.4 Dario 

Una de las ventajas del uso de la teorfa serie/paralelo consiste en que, al asignar un modelo 
constitutivo a cada componente de los compuestos, se puede observar la evolucion de las 
variables internas, tales como el dano o el endurecimiento plastico en el compuesto y en los 
componentes. Por ello, en la simulacion se puede estudiar el dano en el hormigon y en el 
epoxi, lo que permite analizar el mecanismo de fallo de la interfaz de una forma integral. El fallo 
se produjo deslizandose el laminado y desprendiendose una seccion del extremo del bloque de 
hormigon, tal como se aprecia en la figura 5 (a). Al compararla con la figura 5 (b), se observa 
que la parte desprendida del hormigon en la probeta tiene una forma similar a la seccion que se 
deforma en el modelo numerico. 

El mayor dano n el hormigon se produce hacia los bordes superior e inferior de la zona con 
refuerzo, con un maximo del 75% hacia donde termina el FRP, como se puede ver en la figura 
5 (b). Este dano se extiende hasta 0.28m desde el extremo del bloque, es decir, a 0.07m mas 
alia del laminado. Asimismo, en la figura 5 (c) se ve que el epoxi tiene mayor nivel de dano en 
el extremo donde se desprende el hormigon y se va reduciendo hacia el centro de la zona con 
refuerzo. Ademas, en la interfaz con el hormigon el maximo dano es del 22%, mientras que en 
la interfaz con el laminado es del 45%. 
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(c) Dane se al apoxi an al lado coito para 53.2Sk2? 0 0 

Figura 5 Dario en el bloque de hormigon y en el adhesivo 



Como la zona reforzada con mayor dano es la de los bordes superior e inferior, en la 
figura 6 se muestra la evolucion del dano el hormigon y en el epoxi a lo largo del borde 
superior. En lo referente al dano en el hormigon, en la figura 6 (a) se observa que hay una 
distribucion de dano uniforme con un 67% y que el maximo dano ocurre en el extremo donde 
termina el laminado dentro del bloque, con un 75% correspondiente al momento del fallo. 

En lo que respecta al dano en el epoxi mostrado en la figura 6 (b) para una longitud de anclaje 
de 0.21m la distribucion de dano no es uniforme; asimismo, el maximo dano en el epoxi ocurre 
en el extremo del bloque donde se adhiere el refuerzo. La poca uniformidad de la distribucion 
de dano en el adhesivo y su concentracion en el extremo se debe a que, por la corta longitud 
del laminado, hay una mayor concentracion de tensiones. 



Con lo visto en la prueba y en la simulacion se ve la importancia de tener en cuenta el 
dano en el epoxi, lo que muestra que su comportamiento no es elastico lineal como se ha 
supuesto en algunas investigaciones [vii; ix; viii; xxvii]. 
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Figura 6 Dario en hormigon y en el adhesivo en el lado del laminado corto para diferentes 
cargas 

4.5 Tensiones normales y tangenciales en la interfaz 

Otra de las ventajas que tiene la teoria serie/paralelo es que durante todo el proceso de 
carga, permite identificar el estado de tensiones y de deformaciones, tanto en el compuesto, 
como en sus componentes. Debido a que el fallo se produce en la zona reforzada con el 
laminado corto, en este sector se analiza el comportamiento de las tensiones axiales y 
tangenciales del hormigon, del epoxi y del FRP en el borde superior de la zona reforzada, es 
decir, a una altura de 0.15m de la base del bloque. 

En la figura 7 se muestran las tensiones axiales y tangenciales en el hormigon. En lo que 
respecta a las tensiones axiales, en la figura 7 (a) se observa que, hasta una carga aplicada de 
14.48kN, la mayor tension se presenta en el extremo donde se adhiere el CFRP. Luego, como 
serial del deslizamiento, para los otros niveles de carga, la maxima tension, con valores entre 
1.85MPa y 2.0MPa, se produce a lo largo de la zona reforzada hacia el extremo donde termina 
el FRP. Finalmente para la carga de fallo, hay una tension de 2.31 MPa a 0.04m del refuerzo, lo 
que muestra la perdida de capacidad de la zona reforzada. 
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Figura 7 Distribucion de las tensiones en el Hormigon en el lado del laminado corto 

Referente a las tensiones tangenciales indicadas en la figura 7 (b), se observa que para 
los diferentes niveles de carga hasta 49.53kN, las mayores tensiones ocurren en el extremo 
donde se adhiere el FRP, donde se alcanza la maxima tension de 1.68MPa con una carga de 
14.48kN. Luego, para niveles de carga superiores a 49.53kN, las mayores tensiones estan 
entre 0.75MPa y 0.80MPa, y se producen a lo largo de la zona con refuerzo hacia el otro 
extremo laminado. Finalmente, con la carga ultima, donde finaliza el refuerzo, la tension se 
incrementa a 1 .OMPa. 

Al analizar las graficas de la figura 7 se deduce que el comportamiento del hormigon en 
la interfaz corresponde a la accion conjunta de las tensiones tangenciales y axiales. Para los 
primeros niveles de carga (hasta 39.66kN), las tensiones tangenciales alcanzan valores entre 
un 67% y un 89% de las tensiones axiales y luego, para los otros niveles, equivalen entre un 
38% y un 48%. Ademas, se aprecia que las tensiones tangenciales tienen mayor incidencia en 
el dano que ocurre en el extremo del bloque de hormigon que en el lado donde termina el 
CFRP. 



En la figura 8 se observan las tensiones axiales y tangenciales en el epoxi, en la 
interfaz con el FRP donde se produce un mayor indice de dano. En la figura 8 (a) se ve que las 
mayores tensiones axiales ocurren a 0.025m de donde se adhiere el CFRP, hasta una carga de 
49.53kN. Luego las mayores tensiones se dan a 0.05m, donde se alcanza la tension maxima 
de 23.3MPa para la carga ultima de 59.26kN. Por otro lado, al incrementarse el nivel de carga 
en el epoxi, se amplia el area de desarrollo de las tensiones. 
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Figura 8 Distribucion de las tensiones en el epoxi en el lado del laminado corto 

La evolucion de las tensiones tangenciales en el epoxi es aproximadamente parabolico, 
tal y como puede verse en la figura 8 (b). Inicialmente hasta una carga de 27.71 kN, las 
mayores tensiones se presentan en el extremo del bloque. Luego, a medida que se incrementa 
la carga, el sector en el que se desarrollan las mayores tensiones (entre 8.50MPa y 8.80MPa), 
va cambiando hacia el otro extremo donde termina el laminado. 

Las tensiones tangenciales alcanzan valores entre un 34% y un 42% de las tensiones 
axiales, por lo que el comportamiento del epoxi depende de la accion conjunta de las tensiones 
tangenciales y axiales. Sin embargo, el porcentaje de participacion de las tensiones 
tangenciales en el epoxi es menor que en el hormigon. Al comparar la evolucion de las 
tensiones con el dano en el epoxi, se observa que la longitud de adherencia de 0.21m es 
insuficiente para que se desarrolle toda la capacidad del CFRP, ya que el dano se concentra en 
el inicio de la zona. 

En la figura 9 se muestran las tensiones axiales y tangenciales en el CFRP. Las 
tensiones axiales en el laminado observadas en la figura 9 (a), tienen una tendencia casi lineal 
hasta una fuerza de 53.65kN y, despues tienen una distribucion semi-parabolica. Para los 
diferentes niveles de carga, la mayor tension ocurre en el extremo donde se adhiere el 
laminado y la menor se produce donde finaliza el refuerzo. 

Las tensiones tangenciales en el laminado tienen una distribucion de tensiones 
parabolica, tal y como se muestra en la figura 9 (b). Para cargas menores o iguales a 
14.48kN, la mayor tension se presenta en el extremo donde se adhiere el laminado; luego 
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ocurre a 0.05m hasta una carga de 39.66kN. Para los otros niveles, el sector con la mayor 
tension cambia y la longitud de distribucion aumenta, manteniendo un valor de maxima tension 
tangencial entre 1 2.0MPa y 1 2.5MPa. 





|-*-~S3:n -t-'-t-tii-s — -Z--KN jSsSn ;; gw -a-5c 

Figura 9 Distribucion de las tensiones en el laminado del lado del laminado corto 



Al comparar las dos graficas de la figura 9 se puede ver como en la zona donde 
termina el laminado, las tensiones axiales y tangenciales son bajas en el CFRP incluso para la 
carga de fallo, lo que indica que no se alcanza a desarrollar la capacidad del FRP. Por otra 
parte, ya que las tensiones tangenciales son desde un 2% hasta un 4% de las tensiones 
axiales, se evidencia que el CFRP esta solicitado basicamente por acciones axiles. 

Al comparar la distribucion de tensiones y su magnitud, se observa que las tensiones 
tangenciales maximas en el hormigon para cargas mayores a 39.51 kN estan entre un 6% y un 
8% de las tensiones en el laminado, mientras que con respecto a las tensiones axiales la 
relacion es menor a un 0.1%. Por lo tanto, el epoxi, al no ser tan resistente como el laminado y 
al tener mayor resistencia que el hormigon, debe tener la capacidad suficiente para garantizar 
que las tensiones se transmitan y distribuyan en forma adecuada desde el laminado al 
hormigon. 

Coincidiendo con lo encontrado por Chen y Pan [xxvii], en la simulacion realizada en este 
artfculo, al observar el comportamiento del hormigon, del adhesivo y del FRP, se concluye que 
en la simulacion tridimensional, la contribucion de las tensiones axiales y tangenciales es 
diferente a lo indicado en los estudios en los que se ha considerado el estado piano de 
tensiones. En lo que respecta al comportamiento de la interfaz, en el hormigon y en el epoxi, 
tanto las tensiones tangenciales como las tensiones axiales tienen un porcentaje de 
participacion significativo, mientras que en el FRP predominan las tensiones axiales. 
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4.6 Deformaciones 




9M*C CC 


-4- 3M-C 05 


— SII-OlC 


— B- SW-015 


4-SM-C2C 


EX=-C CC 


-*-£> C=-C 05 


— EXPC 10 


15 


-*-E>C=>CX 



Sim: simulacion numerica exp: experimental 

Figura 10 Evolucion de las deformaciones E xx en el compuesto [xviii] 



En la figura 10 se presentan las deformaciones en el compuesto de acuerdo con la ubicacion 
de las galgas en la probeta. Se observa la similitud que existe entre la evolucion de las 
deformaciones medidas en el ensayo y la obtenida en la simulacion a 0.05m, 0.10m, 0.15m y 
0.20m del extremo donde se adhiere el refuerzo. No obstante para 0.0m, con la galga se 
obtuvo una respuesta lineal diferente a lo calculado mediante el modelo numerico, pero la 
pendiente es parecida. Esto probablemente se debe al hecho de que la galga se encuentra en 
la frontera de la zona de adherencia del laminado, mientras que el punto de Gauss analizado 
esta dentro del laminado adherido al hormigon. En este caso, en general, se puede decir que 
en la simulacion y en la prueba experimental, el compuesto tiene una relacion de carga 
deformacion similar. 
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Figura 11 Distribucion de las deformaciones en el compuesto E xx [xviii] 



114 



Maritzabel Molina , Juan Jose Cruz , Sergio Oiler , Alex H. Barbat y Luis Gil 



En la figura 11 se muestra la distribucion de deformaciones obtenida mediante la simulacion 
numerica y el ensayo de laboratorio a lo largo del centro del laminado para los diferentes 
niveles de carga. Se observa que hasta la carga de 56.75kN el comportamiento es similar; a 
partir de este valor, el modelo no Simula el aumento excesivo de deformacion que ocurre a 
0.10m del extremo. Dicho incremento experimentalmente indica la flexion local del laminado a 
causa de una fisuracion localizada. Sin embargo, esta pequena discrepancia no es relevante 
frente a todas las ventajas que tiene el uso de la teoria serie/paralelo en el analisis de 
materiales compuestos. 

5. CONCLUSIONES 

• Debido a que el deslizamiento entre el FRP y el hormigon se produce porque 
sobrepasa la resistencia de la interfaz o del hormigon y, que en ningun caso, se supera 
la resistencia del laminado, el FRP no se esta aprovechando eficientemente como 
refuerzo en elementos de hormigon. La mayorta de las investigaciones estan 
enfocadas a resultados experimentales, que si bien han mostrado aportes importantes, 
no han resuelto incertidumbres esenciales del comportamiento de la interfaz y de la 
perdida de la adherencia. Consecuentemente, es necesario utilizar la simulacion 
numerica como una herramienta complementaria a la experimentacion, para darle al 
proyectista estructural un instrumento que le permita mejorar el calculo de los 
elementos de hormigon armado reforzados con FRP. 

• En este trabajo a traves de la teoria serie/paralelo se estudia el comportamiento de la 
interfaz hormigon-epoxi-CFRP y se comparan los resultados con los obtenidos en un 
ensayo. Igualmente, se modeliza el adhesivo que es susceptible al dano y se considera 
su comportamiento no lineal. Ademas, es importante utilizar modelos tridimensionales, 
ya que permiten calcular el proceso de fallo de una forma integral. No obstante, se 
requiere de una formulacion que permita simular de una forma mas adecuada el 
deslizamiento localizado que puede existir en algunos sectores de la interfaz. 

• La teoria serie/paralelo es una herramienta versatil para analizar el comportamiento de 
materiales compuestos y de sus componentes dentro de una estructura. Por una parte 
permite que cada componente del compuesto se analice mediante la ecuacion 
constitutiva que sea mas conveniente para predecir su comportamiento (elasticidad, 
plasticidad, dano etc.), y por otra, al calcular el comportamiento del material compuesto 
en cada punto de Gauss, lo acopla al resto del mallado, definiendo a su vez, el 
comportamiento global de la estructura. Por tanto, esta herramienta numerica conlleva 
un mejor analisis de estructuras con materiales compuestos sin ocasionar excesivos 
costos computacionales, ya que permite el uso de diferentes tipos de no linealidades y 
utiliza un unico mallado para el analisis de la estructura, de los materiales compuestos 
y de sus componentes. 
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